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Résumé. — On compare la déviation d’un jet par écran déflecteur et par adhérence à une paroi 
convexe, en caractérisant ces écoulements par l’interaction entre la veine quasi-isentropique et les 
zones adjacentes dissipatives (zone de mélange et couche limite le long de la paroi courbe). 

La déviation par écran amplifie le mélange des deux zones ; la déviation par adhérence le mini- 
mise de telle sorte quetle processus gouvernant le phénomène devient le champ de forces centrifuges 
lié à la déviation. Théorie et expériences détaillées. 


Abstract. — Deflection of a jet by a screen is compared to a jet reattaching along a convex wall 
(so-called Younc-Coanp4a Errecr) in terms of a proposed theory for the interaction between the 
nearly isentropic stream and adjacent dissipative flows (mixing zone and curved boundary layer). 

Screen deflection emphasizes the mixing process ; adherence deflection eliminates mixing to such 
an extent that the centrifugal field involved by the deflection appears to become the governing 


process. 


1. Historique. — Le phénomène de déviation 


… d’un jet par adhérence à une paroi convexe placée 
. tangentiellement au jet, est bien connu, depuis 
—. Young {1}, [2] ; il s’agit au début de jets issus de 


PRE NE 


fentes étroites ou d’orifices de faible diamètre 
(fig. 1). 

Chilowski observe à nouveau le phénomène dans 
des expériences sur des projectiles [3], [5]. Étudiant 
l’effet Chilowski, le colonel Lafay arrive à l’idée 
de placer la paroi convexe en prolongement de la 
paroi de la fente d’où sort le jet [3], [5], et utilise 


l'effet propulseur produit par la déviation du jet 
pour réaliser un moulinet. Le renforcement de la 


déviation par grille d’aubes est utilisé sur les freins 
de bouche [4]. 

Les expériences sont reprises à partir de 1932 
par H. Coanda qui propose une paroi convexe à 
facettes planes, et de nombreuses applications 
pratiques de la déviation du jet issu d’une fente 
étroite [5]. 

La déviation par adhérence est observée dans 
différentes structures aéronautiques [6] où elle se 
produit naturellement ; dans ce cas elle est généra- 
lement ressentie comme une gêne, un effet parasite. 


Both theory and experimentation are presented. 


Elle connaît toutefois un regain d’intérêt par son 
application à l’hypersustentation par soufflage au 
bord de fuite combiné avec un volet [7], [8] 
(Shrouded jet flap). 

Une étude mathématique du phénomène a été 
présentée par L. CG. Woods [9] ; nous nous y repor- 
terons en détail dans ce qui suit. 


2. Introduction. — La production simultanée de 
la déviation d’un jet par adhérence à une paroi et 
par un écran intercepteur conjugué présente des 
propriétés remarquables permettant de s’affranchir 
de l’effet d’échelle, de dévier un jet de n'importe 
quelle dimension et d’amplifier les déviations élé- 
mentaires produites séparément par la paroi et par 
l'écran [10] à [15]. 

Au cours de ses recherches sur l’origine de ce 
dernier phénomène [13], l’auteur a été conduit à 
étudier systématiquement la déviation par adhé- 
rence, et la déviation par écran ; il s’est proposé de 
démontrer l'importance considérable du rôle joué 
par diversés interactions entre le jet et le milieu 
ambiant. En particulier, ces phénomènes de dévia- 
tion entrent dans la classe des écoulements carac- 
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FiG. 1. — Déviation d’un jet autour du bout arrondi 
d’une éprouvette, équilibrée par le poids de l’éprouvette. 


térisés par l'interaction entre une veine presque 
isentropique et des zones adjacentes d'écoulement 
dissipatif. Considérant le mélange des deux zones 
par transfert de quantité de mouvement comme le 
processus physique déterminant le gradient maxi- 
mum de pression que l’écoulement peut subir avant 
décollement, L. Crocco et L. Lees [17] ont proposé 
une théorie de cette classe d’écoulements. 

Cette communication a pour but de montrer que 
la déviation par écran amplifie le mélange des deux 
zones ; mais que par contre la déviation par adhé- 
rence le minimise au point qu’il devient négligeable; 
de ce fait, le processus physique gouvernant le 
phénomène devient alors la présence du champ de 
forces centrifuges, lié à la déviation du jet. Une 
théorie semi-empirique de l’action de ce champ est 
présentée ci-dessous. 


3. Étude descriptive. — 3-1. CARACTÈRES GÉNÉ- 
RAUX DE LA DÉVIATION PAR ADHÉRENCE. — La 
condition mécanique essentielle d’une déviation est 
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que la quantité de mouvement déviée soit équi- 
librée par une réaction latérale, qui se présente 
comme une somme : de surpressions quand on 
dévie par écran ; de dépressions quand on dévie par 
adhérence. Mais la différence de base entre ces deux 
modes est que la géométrie de l’écoulement caleu- 
lable en fluide parfait suffit à déterminer l’angle de 
déviation dans le premier mode (par écran) ; dans 
le mode par adhérence, au contraire, l’angle de 
déviation, indéterminé en fluide parfait, est fixé par 
le décollement de la couche limite. Il est donc indis- 
pensable de considérer aussi la zone dissipative le, 
long de la paroi. 


F1G. 2. — Équilibre d’un flux incurvé avec et sans obstacle, 


Étudions l’équilibre du flux incurvé (fig. 2). La 
condition d’équilibre transversal de la portion 
courbe de l’écoulement 


(1 Je) (p[dn) = v?IR (1) 
exprime que les forces centrifuges sont équilibrées 


par un gradient transversal de pression. Supposant 


l’enthalpie totale constante dans tout le jet : 
(dple) + vdv = 0 (2) 


on peut intégrer partiellement (1) le long d’une 
ligne IE normale aux lignes de courant, et on 
trouve la relation : 


= Ya exp f 
Vi E exp |, 


entre la vitesse d'écoulement le long de la paroi V; 
et la vitesse le long de la ligne de jet Vz ; Vz est 
supérieur à Vz. Parallèlement, la pression à la 
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paroi prest inférieure àla pression atmosphérique pz 
Nous supposerons que la ligne IE choisie sur la 


- figure est celle qui passe au point où la pression 


minimum est atteinte. La vitesse à la paroi Vz est 
proportionnelle à V; et croit avec l’intégrale du 
second membre, qui apparaît au premier abord 
comme une fonction croissante de k/R;, en dési- 
gnant par } la hauteur de veine à dévier. 

Cet écoulement du fluide parfait est limité longi- 
tudinalement par l'intervention des forces de visco- 
sité ; la portion déviée de cet écoulement se termine 
à la section du jet qui passe par le point de décol- 
lement et dans laquelle on retrouve partout la 
pression atmosphérique ; le fluide ralenti formant 
la couche limite le long de la paroi est d’abord 
accéléré, mais la survitesse s'accompagne d’une 
dépression qui passe par un maximum, Car elle 
augmente en aval le gradient longitudinal de pres- 
sion que la couche limite doit remonter ensuite et 
facilite ainsi le décollement. Donc la déviation sera 
d'autant plus difficile sur la survitesse nécessaire le 
long de la paroi sera plus élevée (J. Bertin). Consi- 
dérant les conditions aux limites, on remarquera 
que la couche limite est assujettie à remonter globa- 
lement une différence de pression longitu- 
dinale (Pz — Pr) égale à celle qu’on aura réussi à 
établir transversalement. 


3-2. APPLICATIONS. — à) Il en résulte que la dévia- 
tion obtenue par simple rapprochement d’une paroi 
ne sera notable que dans des cas très particuliers, 
soit que V4 soit très faible, soit que la hauteur k à 
dévier soit très faible, même avec une vitesse très 


grande, en prenant À, assez grand (h/{R; — 0,05) : 


il ne suffit pas que A/R, soit faible avec k quel 
conque, car pour dévier de l’angle.x« souhaité, il 
faut que l’écoulement colle sur une longueur À; « 
de paroi ; la couche limite croît avec cette lon- 
gueur À, « ; on ne pourra donc pas dépasser une 
certaine valeur À, max, ni par suite dévier par ce 
seul procédé un jet de hauteur supérieure à un 
certain maximum. 


b) Il est manifeste que ces dernières conditions 
ne se trouvent pas remplies lorsqu'on veut dévier, 
par exemple, le jet d’un réacteur ; la vitesse Vz est 
importante ; la hauteur À imposée de jet à dévier 
est grande. Enfin, le rayon de paroi convexe doit 
être limité de toutes façons à des valeurs telles 
qu’on puisse la caréner convenablement. La dévia- 
tion naturelle est alors très faible, on ne peut pas 
dévier en rapprochant une paroi courbe, et au 
contraire il n’y a pas grand inconvénient à laisser 
cette paroi en place en régime non dévié. Un simple 
examen de cette formule montre que pour het R7 
donnés, il n’y a pas d’autre moyen de diminuer la 
survitesse nécessaire à la déviation qu’en modifiant 
la condition à la limite extérieure Æ de la nor- 
male IE pour y rendre V, le plus petit possible. A 
cet effet, on dispose transversalement à l’écou- 


lement un petit obstacle de hauteur suffisante pour 
créer en amont un champ de pression et un ralentis- 
sement de l’écoulement (fig. 2). Le fait essentiel 
pour les applications est que cet obstacle peut être 
petit et qu’on peut même le remplacer par un jet 
auxiliaire de fluide de faible débit (2 à 3 % “ 
débit principal) transversal au jet principal : 
élimine ainsi totalement l'introduction d’un A 
matériel sur le trajet direct. 


3-3. ACTION CONJUGUÉE D'UN ÉCRAN ET D'UNE 
PAROI DE DÉVIATION. — On peut montrer par des 
expériences simples que la présence d’une paroi 
convexe permet d’amplifier considérablement la 
déviation due à l’obstacle seul, et donc que leffet 
de la combinaison des deux moyens ; obstacle et 
paroi courbe, est un phénomène nouveau remar- 
quable. 


3-3-1. Expérience n° 1. — Pour cette expérience, 
on a utilisé un venturi-tuyère de révolution monté 
à l'envers (fig. 3) : au voisinage du col règne une 


Fi. 3. — Jet dévié par décollement dissymétrique. 


dépression qui peut être très grande (surtout en 
fluide compressible) ; un petit obstacle placé peu 
après le col sur une partie de la circonférence pro- 
voque un déséquilibre de pression et le jet vient 
instantanément se plaquer sur la partie restante en 
éventail : la déviation obtenue est très supérieure à 
celle qu’on aurait obtenue à l’aide du seul bord 
évasé de tuyère ou du seul obstacle. 


3-3-2. Expériences de J. Bertin. — Une autre 
expérience montre également de façon frappante 
comment l’action dynamique du bord de déviation 
peut être transformée en action cinématique : la 
figure 4 montre un venturi rectangulaire dont les 
parois divergentes sont mobiles autour d’axes situés 
au voisinage du col de l’écoulement. Sous l’effet des 
dépressions qui règnent le long de ces parois 


k A 


mobiles, celles-ci tendent à se rapprocher si on les 
laisse libres. Si on n’en libère qu’une, maintenant 
l’autre immobile (cas dela figure), la position d’équi- 
libre correspond à un Jet fortement dévié autour 
de la paroi immobilisée. On notera la fermeture très 


— Jet dévié par le déséquilibre d’un Venturi 
à paroi mobile. 


Fic. &. 


élevée de la section de tuyère obtenue sans blocage 
puisque la paroi aspirée est libre dans son mou- 
vement horizontal. 

Nous ne nous étendrons pas sur l’explication 
détaillée de ces phénomènes, qu’on trouvera dans 
les références [11], [12], [13]. 

3-4-1. Considérons (fig. 5) l'équilibre du jet non 
dévié entre la couche limite intérieure et la zone de 
mélange extérieure et appliquons un raisonnement 
classique ([16], t. 1, p. 85) : lorsqu'on introduit 
l'obstacle, on perturbe cet équilibre par ralentis- 
sement des filets fluides extérieurs dont les parti- 
cules sont au surplus rapprochées de la paroi de 
déviation ; une particule extérieure ainsi déportée 
transitoirement vers une zone où le champ de 
forces centrifuges est plus intense, alors que sa 
vitesse propre est insuffisante, n’est pas en équi- 
libre et, sous l’effet du gradient de pression corres- 
pondant, elle se rapproche encore davantage de la 
paroi de déviation ; au contraire, une particule 
rapide, en aval de l'obstacle, déportée par la turbu- 
lence dans la zone de mélange avec l’atmosphère 


n 


composée de fluide ralenti, s’y trouve soumise à une 
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F1G. 5. — Établissement d’une circulation 
par passage du régime direct au régime dévié par écran. 


force centrifuge non équilibrée par un gradient de 


‘ pression et tend donc à se déporter encore plus loin. 


La première action décrite a pour action de dévier 
davantage les filets fluides, les deux ont pour effet 
d'augmenter la largeur de la zone de mélange exté- 
rieure. De plus, cette zone de mélange est instable 
puisque les échanges de quantité de mouvement 
entre fluide rapide et fluide lent qui lui donnent 
naissance ont tendancé à augmenter. En résumé, le 
coefficient de mélange [17] avec le jet est nota- 
blement augmenté. La déviation finale résulte de la 
position d’équilibre prise à la limite intérieure, où 
les phénomènes sont exactement contraires. Une 
particule rapide déportée transitoirement par la 
déviation vers l’intérieur de la couche limite s’y 
trouve soumise à une force centrifuge non équi- 
librée par un gradient de pression et s’éloigne donc 
à nouveau de la paroi de déviation ; de même, une 
particule ralentie de couche limite déportée dans le 
fluide rapide, ne possédant pas de force centrifuge 
suffisante, « retombe » dans la couche limite. Done 
l'écoulement a une position d'équilibre stable, ainsi 
que sa couche limite, et de plus, l’épaisseur de la 
nouvelle couche limite est plus petite que l’ancienne 
puisque les échanges de quantités de mouvement 
entre fluide rapide et fluide lent sont amoindris par 
la déviation. En conséquence secondaire, le point de 
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décollement est reporté notablement plus loin en 
aval, donnant ainsi une mesure de la déviation 
obtenue, du moins localement. 


3-4-2. Cas d’une déviation par adhérence seule- 
ment.— L’explication physique qui précède permet 
de rendre compte de la même façon de l’effet d’une 
paroi de déviation seule : considérons le cas d’un 
jet de faible hauteur, devant laquelle l'épaisseur 
de la couche limite ne soit pas négligeable. Le fait 
de remplacer d’un côté du jet une zone de mélange 
instable par une couche limite stable et amincie, 
explique un certain déséquilibre qui peut se tra- 
duire par une déviation très importante du jet. En 
outre comme la hauteur du jet est faible, le champ 
des pressions occupe presque toute la zone isentro- 
poque. De ce fait, le zone dissipative du côté 
opposé à Ja paroi se trouve profondément 
modifiée [7]. Le champ des dépressions agit comme 
une répartition de puits le long de la paroi et 
aspire fortement l’air ambiant, dont l’écoulement 
est perpendiculaire à la paroi, et même à contre- 


yes 


courant du jet. Les échanges de quantité de mou- 
vement augmentent fortement pour les raisons 
décrites au $ 3-4-1, parmi lesquelles laspiration 
par le champ de dépressions prend la première 
place. 


3-4-3. Circulation. — Dans les deux cas, le pro- 
cessus décrit a pour effet d'établir autour du dévia- 
teur une circulation finie correspondant à la poussée 
latérale obtenue. 


4. Expériences de déviation combinée par écran 
et paroi. — 4-1. RÉSULTATS EXPÉRIMENTAUX. — 
Les considérations précédentes ont été vérifiées et 
précisées par des expériences sur des Jets issus de 
bouches d’air comprimé. La tuyère se rapprochait 
le plus possible du schéma d’écoulement bidimen- 
sionnel dont on avait calculé la géométrie, en effec- 
tuant les corrections de compressibilité. Comparant 
les répartitions théoriques et expérimentales de 
pressions, on en déduit le rôle joué par les phéno- 
mènes dissipatifs (viscosité et turbulence) (fig. 6). 


——.—— Avec paroi 


mere Sans paroi 


Théorique 
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Pression génératrice : 
17,8 pièzes 
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a) Déviation par écran seulement. — L'efficacité 
du déviateur de jet subit deux pertes d'importance 
comparable : décollement de la couche limite et 
zone de gaz morts côté obstacle ; détente incom- 


‘plète, même à l'extrémité, côté opposé. Alors que le 


régime est subsonique, la pression reste supérieure à 
la pression atmosphérique. 

b) Déviation par écran et paroi convexe adaptée. — 
(C'est-à-dire une paroi circulaire placée en gran- 
deur et position sur le cercle de courbure minimum 
de la ligne de jet.) 

La paroi est notablement en dépression ; de plus 


la pression baisse notablement le long des parois 
amont de la tuyère qui subit une détente complé- 
mentaire. L’anomalie signalée en (a) disparait ; 
le coefficient de striction (facteur de blocage) de la 
veine se rapproche de l’unité, c’est-à-dire que le 
blocage diminue considérablement ; la déviation 
augmente d’une vingtaine de degrès ; on observe un 
important entraînement d’air ambiant (du côté de 
la paroi). En résumé la paroi agit comme une 
trompe sur le jet dévié par l’obstacle ; elle abaisse 
la pression statique à l’extrémité du canal, action 
équivalente à une ouverture de tuyère. 


6 A 


c) Déviation par obstacle et paroi non adaptée. — 
.Si la paroi est placée de telle manière qu’elle morde 
sur la ligne de jet théorique (trop rapprochée, ou 
rayon trop grand) elle se trouve en surpression ; la 
déviation diminue ; le coefficient de striction dimi- 
nue (le blocage augmente) ; il y a donc fermeture de 
tuyère, et remontée de pressions statiques le long 
des paroïs de tuyère. 


4-2. Discussion. — a) Déviation par obstacle seul. 
— Les décollements du côté écran étaient prévi- 
sibles. La remontée de pression du côté opposé peut 
être expliquée comme due à l'interaction de la 
couche limite et de la zone isentropique dans la 
tuyère. Ce phénomène est l’inverse de la déviation 
par adhérence. 

S'il n’y a pas d’obstacle, la surpression est 
nulle : il y a même une petite dépression à l’extré- 
mité de la tuyère. L’obstacle dévie les filets fluides 
isentropiques vers le bas, alors que la couche est 
astreinte à suivre la paroi plane. Il en résulte une 
augmentation du coefficient de mélange de la zone 
isentropique et de la zone dissipative [17], qui a 
pour effet une augmentation de pression statique. 
La pression à la paroi est égale à la pression à la 
frontière des deux zones, sur une même perpendi- 
culaire à la paroi. Cette frontière entre les deux 
zones a été reportée plus haut par suite de l’aug- 
mentation du coefficient de mélange ; donc le 
niveau de pression statique est relevé, eu égard à la 
direction du gradient de pression. 


b) Déviation par obstacle et paroi de déviation 
adaptée. — Le long de la paroi de déviation, le 
phénomène est complètement différent ; la couche 
limite est astreinte à suivre la courbure de la paroi, 
donc elle est soumise au champ des pressions qui 
équilibrent les forces centrifuges (qui n’agissent 
pas sur elle puisqu'elle est composée de fluide ra- 
lenti) : le raisonnement du $ 3-4-1 peut être repro- 
duit ici, à la condition que le champ des pressions 
qui diminuent les échanges de quantité de mouve- 
ment s’établisse réellement : ce qui exige une posi- 
tion correcte de la paroi. 

c) Si au contraire, la paroi mord sur la ligne de 
jet, le raisonnement n’est pas valable, et on est 
ramené au cas (a) précédent : il se produit des 
surpressions et un décollement prématuré de la 
couche limite. 


5. Expériences de déviation par adhérence pure. 
— 5-1.— La déviation obtenue la seule adhérence à 
une paroi a été étudiée systématiquement, en vue 
de déterminer quantitativement ses caractéris- 
tiques et ses limites. C’est un exemple typique 
d'interaction d’un écoulement isentropique et d’un 
écoulement dissipatif adjacent, puisque la paroi ne 
peut agir sur le jet qu’au travers de la couche 
limite. D’après ce qui précède, le rôle essentiel sera 
joué, non par le coefficient de mélange, mais par le 
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champ des pressions, donc d’après l’équation (3) 
par la courbure relative hfr de paroi, paramètre 
dont dépend principalement la différence trans- 
versale de pression (px — pr) égale à la différence 
de pression longitudinale que la couche limite est 
astreinte à remonter. 

La plupart de nos expériences ont porté sur la 
configuration obtenue en prolongeant une paroi de 
la tuyère rectangulaire par une paroi en arc de. 
cercle de 90 degrés sans solution de continuité. 


5-2. DEscriprion pes Essais (fig. 7). Notations. 
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La tuyère expérimentale rectangulaire de hau- 
teur h, largeur Z, comporte un canal droit de 
longueur L' prolongé sur une des parois par un 
quart de cercle QR de rayon r. Le jet est limité par 
des panneaux latéraux sur deux faces ; il est libre à 
partir de O sur la 3€ face, et à partir du point de 
décollement D sur la face qui porte la « paroï de 
déviation ». Prenant O comme origine, on peut 
définir des axes Ox, parallèle à OB, et Oy perpen- 
diculaire. La réaction latérale est dirigée suivant Oy 
et la traînée suivant Ox. 


5-2-1. On mesure la déviation par l’arc QD. La 
déviation maximum est limitée par construction à 
90 degrés. On pourrait facilement obtenir une 
déviation plus grande, mais cette limitation artifi- 
cielle fait apparaître des phénomènes intéressants 


comme on le verra ci-après. On mesure les pressions . 


le long des parois ; elles sont inférieures à la pres- 
sion atmosphérique, la dépression passant par un 
maximum en un point M du bord de déviation. Les 
dépressions mesurées sont très diminuées, ainsi que 
la déviation, si le jet n’est pas limité par des 
panneaux latéraux. L'alimentation par l’atmo- 
sphère de la zone RDA, est donc ici un phénomène 
parasite, qu’on a évité dans cette étude expéri- 
mentale. 

On fait varier toutes les dimensions géomé- 
triques, et, dans une certaine mesure, le nombre de 
Mach du jet. 

La technique des mesures n’appelle aucune 
remarque particulière. On appelle : 


\ 


à 


A 


N°4 


P, la pression totale génératrice ; 
p la pression en un point courant de la paroi de 


. déviation ; 


Pa la pression atmosphérique ; 
w la vitesse au point courant ; 
w, la vitesse à l'infini aval, ou par détente à la 


. pression atmosphérique ; 


D le débit-masse de la tuyère. 


5-2-2. Pour la représentation des résultats, outre 
les arcs de cercle Yy et y», on définit le coefficient 


… de poussée local : 


je DEEE 
__ Pa P Que ERA (4) 


Po — Pa wà 


Presque tous les essais ont été effectuées à 
M = 0,48 ; les corrections de compressibilité ne 
dépassant pas 2 %, nous pouvons interpréter les 
résultats en supposant le fluide incompressible. 

K passe par un maximum X7, en M : il n’est 


- jamais nul, même le long de RD, en raison de 
” l'existence d’un courant de retour RDA, qui pro- 


duit une petite dépression en aval de D. 
La poussée latérale est donnée par l'intégrale des 


pressions Ÿ (résultante verticale). 


Rapportée a la pression dynamique elle permet 
de définir le coefficient latéral moyen : X, ; c’est 
un coefficient de portance. On peut également 
rapporter la poussée latérale à la poussée globale : 


Dw,, définissant un coefficient K, de poussée laté- 


rale. Les deux coefficients it liés en fluide 


 incompressible par la relation : 


(L' + r) Kp — 2ehKy #2hKy 
« (l'effet de striction est négligeable en l’absence 
. d’obstacle). 
5-2-3, — Les variables mesurées ou déduites 


des mesures sont donc : 
Yu; Yp ; K(x+) sur chaque paroi ; Kw ; Kp ; Ky. 


Les paramètres dont on a 1 pRorCRe à déterminer 


; le rôle sont : … 


au été": de 7 


| 
| 
| 


r : 10 60 80 120 160 200.240. 280. 480. 640. 2 000 mm. 
h : 10 20 30 40 60 80 120 160. 200. 240. 280 mm. 


L : 120 mmet 54 mm. 
_L\: 0 100 200 et 300 mm. 


M : de 0 45 à 0.75. 


5-3. THÉORIE SEMI-EMPIRIQUE, POUR UN CRITÈRE 
1 
D'ADHÉRENCE. 


9-3-1. — D’après l’analyse faite au $ 3 (établis- 
sement d’une circulation) on peut donner une des- 
cription élémentaire de l'influence du bord de dévia- 
tion sur le jet en profondeur, qui correspond assez 
bien aux phénomènes observés : l'influence est 
nulle si le jet sort droit parallèlement à Ox ; elle 
serait complète si les lignes de courant devenaient 
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des portions de cercles concentriques à QR. Ima- 
ginons donc que dans le cas intermédiaires les 
lignes de courant puissent être assimilées à un 
faisceau de cercles (fig. 8a), du moins dans la zone 
re EAN 
Q G = 


= 
| N 


Fic. 8a (17€ phase). 


très voisine de la sortie QO, car le point M reste 
très voisin de Q : supposant M et Q confondus, 
nous trouverons par cette théorie locale une expres- 
sion de Ky. 


5-3-2. 1er cas. — Le faisceau va de la droite Ox 
au cercle de diamètre : 


OOR=NT 


en désignant par Àr le rayon de courbure juste à 
l'extérieur de la couche limite. Soit a? la puissance 
du faisceau : 


a = h(h + 2àr) = 0Q.0Q. 


En un point P quelconque de OQ (OP = y), le 
rayon de courbure est tel que : 
2R = (a?[y) — y. 
L'intégrale formule (3) donne alors : 
wo = Wall + (R]2àr)] = w. (5) 


Lorsque À varie de l'infini (influence nulle) à 1 
(couche limite nulle) par exemple par une augmen- 
tation progressive du rayon r, Wy varie de 4, à Wa 
(1 + A/2r) : dans cette première phase, l'effet du 
bord de déviation se réduit à l’accélération de la 
couche limite, mais ne se traduit pas par une dévia- 
tion globale de l'écoulement. 


5-3-3. 2€ cas. — En dessous d’une certaine cour- 
bure relative, la couche limite est suffisamment 
réaccélérée (À voisin de 1) pour qu'une augmen- 
tation nouvelle de r (ou une diminution de h) 
commence à avoir pour effet un recollement pro- 
gressif de l’écoulement sur une certaine longueur 
de paroi (M reste voisin de Q, mais D se déplace 
vers l’aval). Dans ce cas il existe une déviation 
réelle, que nous traduirons par la considération 
d’un faisceau de cercles limité du côté extérieur 
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par une ligne de jet de rayon fini r’, et du côté inté- 
rieur par le cercle QQ' de rayon r. 
Soit À (fig. 8b) le centre du faisceau que nous 


(E 


> 


re 
| \ 
| 
| 


Fic. 8b (2° phase). 


prenons comme origine, en posant : AO = y, ; 
QS = y; AP =Y;%o—Y1—h. 


Au point P le rayon de courbure est encore. 


donné par : 
2R = (a*[y) — y. 
Le calcul donne alors sans difficulté pour la sur- 
vitesse maximum : 
r'2r +h 
rgr —h 


wQ Wa = nf: 

Les cas extrêmes à considérer sont : r° infini : 
on retrouve la formule précédente pour À = 1 ; 
r' = r + h (écoulement giratoire) ; dans ce cas : 


WmM = Wall + hr). (5°) 


5-3-4. — Ces calculs conduisent à mettre sous la 
forme suivante, la relation entre la survitesse 
maximum et la courbure relative : 


(wm — Wa) [wa = œhfr = ac = (1 + Kw)? — 106) 


Le coefficient « caractérise le type d’adhérence 
obtenu : s’il est compris entre 0 et 0,5, l'influence 
du bord de déviation reste locale et intéresse 
d’abord la couche limite ; s’il est compris entre 0,5 
et 1, on peut constater un recollement et une 
déviation globale du jet, l’évolution progressive 
vers l’écoulement giratoire marquant l’apparition 
d’une circulation (cf. $ 3). 


5-3-5. — Soit c, la valeur de hfr pour laquelle 
x — 1/2 ; une expression de « dans le 2° cas, où la 
variation de € n’agit plus sur la couche limite, peut 
être déduite de considérations de similitude. 
On trouve [13], facilement, pour € < €, : 
wM | cél(e + 1) (ce + 2) 
a: 1 VAE actu 
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5-3-6. — Dans le premier cas, la survitesse 
dépend de considérations sur la couche limite ; 
cependant le raisonnement du $ 5-3-2 suppose 
implicitement l’existence locale d’une loi de simi- 
litude, c’est-à-dire que le rapport À des rayons de 
courbure ne dépend pas de À et r séparément mais 
seulement de c ; en outre l’angle y, avant décol- 
lement (en supposant toujours M et Q confondus) 
ne varie pas beaucoup pendant cette phase. Il est 
déterminé par le degré de divergence que la section 
est susceptible de supporter avant décollement, 

. F 
SOI : À ñ 
est égale à : 


(1 — cos yÿs) ; cette variation de section 


(wo — Wa) [Wa 


en désignant par w, la vitesse moyenne dans la 
section OQ. Cette dernière variation est elle-même 
une fraction de la survitesse : (@y [Wa) — 1. 

On peut donc écrire y, étant petit : 


(ar [h) (Y312) = k2h[22r. 


D'où une loi de la forme : À = kc/y) pour c > & 
et par suite : #y —= Will + (yp/k)] ne dépendra 
pas de c dans cette phase. 


5-4, — RÉSULTATS D’ESSAIS. 


9-4-1. — L’allure générale des phénomènes est 
bien celle qui a servi de base au calcul précédent : 
les paramètres les plus intéressants sont k et r. 
La répartition des pressions se déforme en fonc- 
tion de c comme l'indique la figure 9 (r constant, 
h variable): a) Courbure élevée : Forte dépression 
atteignant son maximum en M très près de Q, et 
par suite fort gradient de pression au delà de M, 
donc décollement en D très près de M ; faible varia- 
tion de cette courbure en fonction de k. b) Courbure 
modérée : À partir d’un certain seuil quand À dimi- 
nue, la courbe des pressions s’évase, le point D se 
déplaçant vers l’aval progressivement. Le maxi- 
mum reste élevé, mais la forte dépression se main- 
tient sur un arc croissant. Le point M se déplace 
aussi vers l’aval, mais la dépression atteint presque 
son maximum très près de Q. c) Courbure faible : 
Lorsque la déviation maximum (90 degrés) a été 
atteinte, une nouvelle diminution de À ne fait que 
diminuer les dépressions sans améliorer l’adhérence 
qui est déjà très près de l’optimum (écoulement 
giratoire). 


5-4-2. Variations de la dépression maximum 


(coefficient ÆX3). — Nous nous reportons à la 
définition du critère d’adhérence : 
ac — (1 + Kay) /2 — 1 — survitesse. 
SPA Aer ac = Constante = 05 © 
2 En a < 7 
SL Te PORÉSPS dub: s (7) 


c(e + 2) + ci 
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RE I ë T = 120 mm 

Fe w? = 17,5 pièzes 3 

À 12 
J 

| () 

A an 

2 o2:°? 2 

# ae 

= d È 2 

DT PS 
V4 2 Le 
Lg F7 gr | 

PNY Fa 
= =D ; ES 5 

/ L 

7 
4 ]-6 
KA | 

EN y, 
1-8 
1- &] 
410 


:—.— h = 80 


Points expérimentaux : 


Dépressions maxima 


+ À — 280 
o ” = 160 
v , — 120 
Ab 080 


suivant formules (7) et 
correction (8) 


À - a Courbes théoriques extrêmes 


Asym ptote 


1,0 1,5 2,0 25 C 
Fi. 10. — Dépressions maxima. 
Points expérimentaux : X À — 280 
» O — 160 
» \v —M20 
» A 80: 
Courbes théoriques extrêmes ( ————h} — 280. 
suivant formule (7) et 
correction (8) | ——— .— h 80. 


Cette courbe passe par un maximum pour ce —c;. laquelle : & 


— 0,5. La figure 10 montre que pour 


Pour déterminer, c,, on trace la courbe expéri- € < Co, l'expression théorique donne une repré- 
mentale œ(c) et on cherche la valeur de c pour  sentation très honnête des résultats, si l’on tient 
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compte de la fragilité des hypothèses de similitude ; 
c, a les valeurs indiquées sur cette figure. 

Pour c > c, la survitesse n’est pas constante 
mais décroît lentement, suivant une loi de la 
forme : 


Li RC MNT CINE CEE NUOISE— Re (8) 


Les coefficients X, X”, c, déterminés expérimen- 
talement et liés par la relation ci-dessus sont 
influencés par les différents paramètres : k,r, M, L 
et L'. Seul L' semble avoir une importance appré- 
ciable. 

Dans la mesure où les couches limites jouent un 
rôle secondaire (valeurs raisonnables de L et L'), 
on peut admettre * 


K — 0255 (+ 6 %) 
K'—116 (+ 4%) 
co = 056 (+4 %). 


5-4-3. Variation de la dépression moyenne et de 
la poussée latérale. — La figure 11 montre la varia- 
tion de X, en fonction de c, très analogue à celle 
d’un Cz de profil en fonction de l’incidence, ou 
mieux du Cz d’un cylindre tournant en fonction de 


la survitesse périphérique de rotation (effet . 


Magnus : [18]). 


Hnmnm 
HU 
t 
& 
4Ox0 
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L'effet d’adhérence correspond à la partie 
linéaire, le long de laquelle la déviation est de 
90 degrés. Il passe par un maximum en fonction 
der, ainsi que le pourcentage de poussée latérale X,,: 
cet effet dépend de l’échelle, et diminue quand r 
est trop grand (développement de couche 
limite : $3). 

Il est intéressant de noter la variation de X, 
(fig. 11) en fonction de fr ; le pourcentage de 
poussée latérale relativement faible pour une cour- 
bure élevée, augmente brusquement quand l'effet 
d’adhérence commence à se produire, puis reste à 
peu près constant ; en raison de l’étroïte relation 
entre ce coefficient K,, la déviation y», et la pro- 
duction de circulation, c’est une confirmation 
remarquable de analyse donnée au $ 3 sur laquelle 
est fondée la théorie semi-empiriques $ 5-3. 


6. Théorie de l’effet Young. — Le calcul de la 
déflexion d’un jet autour d’un cylindre circulaire 
a été fait par L. C. Woods [9]. Comme la forme de la 
paroi est imposée, on est conduit à une équation 
intégrale en w ; on calcule une expression appro- 
chée en remplaçant w sous le signe intégrale par sa 
valeur moyenne, notée w, le long de l’arc d’adhé- 


Fre, 42: 


rence D'OD de longueur 27 (fig. 12). On admet en 
outre le long de cet arc que : 


Yo = 0, p—=#wr6; pp = wr. 


en désignant par 
0 l’angle de la vitesse avec Ox ; 


le potentiel au point P tel que OP — r0. 
Si l’on pose : 


Re PE ne (9) 
La vitesse au point P est donnée par : 


w dé o 2 
log — à Arc tg { sh? o— ch? o th? o(0/y) }/2, (10) 


6-1. COMPARAISON DES GRADIENTS DE PRESSION 
LONGITUDINAL ET TRANSVERSAL, — La courbe w/w, 
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ou X en fonction de 0 a des tangentes verticales 
pour 0 — + y (fig. 13) en raison des singularités 
du type (0 + y)? dans le gradient de pression 
Cette propriété très importante serait une source de 
difficultés dans un calcul de couche limite ; mais 
comme elle permet de rendre compte de l’allure de 
Pécoulement, il faut lui attribuer une signification 


15 Y 
HTGAETE: 


physique : les singularités résultent de l’égalité des 
différences de pression entre O et D d’une part, 
O et O’ d’autre part, sans qu’il soit nécessaire de 
faire intervenir la couche limite : dans le plan de 
Phodographe, D et D’ sont des points de rebrous- 
sement. Le rapport des gradients moyens de pres- 
sion correspondants est directement proportionnel 
à © ; d’où la discussion qui suit : 


6-2-. GRADIENT TRANSVERSAL FAIBLE PAR RAP- 
PORT AU GRADIENT LONGITUDINAL ; G PETIT. — 
C’est le cas du $ 5-3-2. Pour fixer les idées calculons 
un maximum de oc d’après la figure 9 : r — 120; 
= 00 #1, Le, 


Y = YÿD — ym = 0 35 6 — 12. 


Pour o < 2, l'expression de la vitesse maximum 
est peu différente de : - ï 


Wo —= Wa €Ÿ. 


Elle est indépendante de À et r, et même de w 
(cf. $ 5-3-6). Les écarts constatés avec les résultats 
d'expériences ($ 5-4-2) sont donc imputables à la 
seule viscosité. La valeur de y n’étant pas négli- 
geable et fixant l'intervalle d’abscisse entre les 
tangentes verticales, la courbe X(6) a bien l’allure 
pointue de la figure 9. 


6-3. GRADIENT LONGITUDINAL FAIBLE PAR RAP- 
PORT AU GRADIENT TRANSVERSAL : G GRAND. — 
C’est le cas du $ 5-3-3 où la circulation apparait, 
ainsi que l'effet Young. Calculons un ordre de 
grandeur minimum d’après la figure 9 : 


O2 0 y = [20 Lo > 10. 
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L'expression de la survitesse maximum est done 
très peu différente de : 


wo 


log 


a | 
& 
a 
nl 
En 
& 
& 
na 


T 
Wa 2 
# est voisin de w,, et w, diminue parce que l’inter- 
valle d’abscisse entre les tangentes verticales 
s’allonge tandis que log w,/w, reste inférieure à Ar. 

La courbe des pressions a donc nécessairement 
l’allure aplatie, à forts gradients aux extrémités, 
notée sur la figure 9. On trouve entre cet le critère 
d’adhérence « la relation (c petit) : 


(1 + ac) log (1 + ac) = c. 


6-6. EXPRESSION DU GRADIENT LONGITUDINAL. 
— Si l’on se fixe l’angle y, l’écoulement corres- 
pondant sera nécessairement décollé en D par le 
gradient local infini. Pour savoir dans quelle 
mesure le développement de la couche limite peut 
intervenir autrement, il faut considérer le gradient 
de pression en amont de D, par exemple au voisi- 
nage de O. Ce gradient a pour expression : 


dæ w2 Y ] 0 
— © — ee DRE 0 
#75 + (5) exp (2 Tu Hip Cote ve (11) 


en posant : 


tedu— (sn: 6 — ch? 6 tg? 5.) 


a) 6, petit : ÀAu voisinage de O, le gradient a 
pour expression approchée : 


b) 6, grand : Au voisinage de O, l'expression 
approchée devient 


0 

th 6, — 

dw Fe __ wà x) ? Ve 
mg = 8 = Sexp(rs sh 62 


Pour une même valeur de £: le rapport de ces 
deux gradients est très grand, de l’ordre de : 
(81182) = Sh 6: exp (2y1 — (T{62) v2). (12) 

Par exemple, pour : 


Met DS PO on 10 


? 


on trouve : 


&g1 = 13 000 go. . 


Le gradient g, est donc négligeable à cause de la 
valeur élevée de 6, : pour la couche limite, tout se 
passe comme si le jet était droit et la paroi plane ; 
les longueurs parcourues avant décollement peuvent 
être grandes ; c’est en cela que consiste le phéno- 
mène d’adhérence. 


6-5. ConcLusion. — L’interaction d’un jet et 
d’une paroi courbe peut se présenter sous plusieurs 


12 A 
aspects, en fonction des facteurs caractéristiques : 
AS 
CR NE ER 
k 2 Wa C 


La couche limite le long de la paroi doit remonter 
un gradient de pression qui se présente princi- 
palement comme une fonction exponentielle 
de (— 6) au point P divisant l’arc OD dans un 
rapport donné. 


a) Écoulement décollé. — Au-dessus d’une pre- 
mière valeur critique © — €,, le gradient de pres- 
sion est élevé et pratiquement indépendant de c. 
L’interaction reste locale, et l’arc d’adhérence reste 
faible et à peu près constant. On peut admettre : 
Y1 # 20 degrés # 0,35 radian. Parallèlement, 6 est 
inférieur à une première valeur critique 6, de 
l’ordre de l’unité : e, est l’homologue du rapport c, 
défini au $ 5-3-2. 

b) Zone intermédiaire. — En dessous de c;, la 
couche limite peut remonter en partie le gradient 
de pression ; y augmentant lorsque c diminue, 
oc augmente rapidement ; le gradient de pression 
diminue en fonction exponentielle de 6, d’où aug- 
mentation accélérée de y. [Il y a donc évolution très 
rapide du régime (a) vers le régime (c) très différent. 

c) Effet Young. — Il est caractérisé par l’annihi- 
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lation du gradient de pression le long d’un grand 
arc de paroi courbe, lorsque c devient inférieur à 
un deuxième rapport critique c,, auquel correspond 
une valeur pseudo-critique 6. Toutefois, l’exis- 
tence de ces valeurs pseudo-critiques est liée à celle 
d’un maximum par construction de l'angle . 
Pour y: — 90 degrés, 6, est de l’ordre de 10. 


d) Effet d'échelle. — Les valeurs numériques 
données ci-dessus correspondent aux conditions de 
vitesse des expériences (M — 0,48). Lorsqu'on aug- 
mente la vitesse, toutes choses égales par ailleurs, le 
point de décollement remonte vers l’amont ; il faut 
diminuer À pour retrouver l’effet Young, les sur- 
vitesses relatives diminuant en conséquence 
($ 5-4-2-2), cependant que 6, augmente. Le sens 
de la variation de 6,est moins certain. L'effet Young 
dépend donc du nombre de Reynolds : w, r ywi, 
suivant une loi qu'il serait intéressant de déter- 
miner. 

L’exploration que nous avons pu faire est incom- 
plète et ne couvre pas les vitesses usuelles à l’échap- 
pement des moteurs à réaction. L'action conjugués 
d’un écran et d’une paroi de déviation ($ 3) peut 
être expliquée très clairement par cette description 
finale [13]. 


Manuscrit reçu le 7 janvier 1958. 
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ÉTUDE DES SELS FONDUS PAR DIFFRACTION DES RAYONS X 
AUX TEMPÉRATURES ÉLEVÉES (Suite) 


II. STRUCTURE A L’ÉTAT LIQUIDE DES CHLORURES LiCl, NaCI, KCI, BaCI; 
ET DU FLUORURE Car.. 


CONSIDÉRATIONS GÉNÉRALES SUR LA STRUCTURE DES HALOGÉNURES FONDUS 


Par JEerzy ZARZYCKI, 


Centre de Recherches des Glaceries de Saint-Gobain (France). 


Résumé. — L’analyse des courbes de distribution radiale déduites des spectres de diffraction X 
des sels fondus au moyen de la transformée radiale de Fourier confirme les résultats déjà obtenus 
pour les fluorures alcalins : à l’état liquide il y a conservation de l’ordre à courte distance du cristal 
mais décroissance simultanée du nombre de coordination principal. Ceci prouve la structure lacu- 
natre des halogénures fondus. On montre que ces lacunes constituent en réalité des « fissures fluc- 
tuantes » entre des amas ioniques ou « îlots » formés en moyenne par 1-2 couches d’ions entourant 
un ion central donné. 


Abstract. — The analysis of the radial distribution curves deduced by Fourier transform from 
the X-ray diffraction spectra of molten salts confirms the results already obtained for alkali fluo- 
rides : in the liquid state the low-range order of the crystal is conserved but the first coordination 
number is much lower. This proves the lacunar structure of molten halides. It is shown that 
the holes are actually ‘‘ fluctuating cracks ” 
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around a central ion. 


I. Introduction. — Dans les articles précé- 
dents [1], [2] nous avons étudié la structure à l’état 
liquide des fluorures alcalins LiF, NaF et KF. 
Comme suite de ces travaux, nous allons présenter 
maintenant l’étude de quelques autres halogénures 
alcalins et alcalinoterreux fondus, notamment celle 
des chlorures LiCl, NaCI, KCI et BaCI, ainsi que 
celle du fluorure CaF,. 

L'étude des chlorures LiCl et KCI liquides a été 
déjà tentée par Lark-Horowitz et Miller [3], [4]. 
Ces travaux se situant avant la mise au point de la 
méthode générale d'interprétation utilisant l’inver- 
sion de Fourier dans le cas des corps amorphes 
hétéroatomiques, les auteurs cités ont dû introduire 

_ dans leurs calculs certaines simplifications impor- 
tantes concernant les facteurs de structure 
(fit = 0 dans le cas de LiCI ; fret = fa— dans le 
cas de KCI) pour se ramener au cas des corps 
monoatomiques ce qui a dû influencer certainement 
Je résultat final recherché. 

Nous avons donc jugé intéressant de reprendre 
dans son ensemble l’étude des chlorures avec les 
moyens plus perfectionnés et les méthodes d’inter- 
prétation plus rigoureuses dont nous disposons 
actuellement. 

En complétant ainsi les résultats déjà obtenus 
concernant les fluorures nous tenterons de donner 
une image plus précise de la structure des halo- 
génures fondus. 


IT. Résultats expérimentaux et analyse des 
spectres. — En ce qui concerne la méthode expéri- 


between ionic clusters grouping 1-2 ionic shells 


mentale utilisée on se reportera aux articles parus 
précédemment [5], [1]. 


200 


100 
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Fic. 1. — Spectre de diffraction de LiCIl fondu à 620 °C. 


Les figures 1 à 5 montrent les spectres de dif- 
fraction de différents sels maintenus à l’état liquide 
aux températures respectives suivantes (à -- 20 0C 

rès) : 
LiC1 à 620 °C, NaCI à 820 °C, KCI à 810 °C, 
BaCI, à 1 000 0C et CaF, à 1 500 °C. 

Les spectres ont été enregistrés en employant la 
radiation MoX%«, strictement monochromatique, 
les temps de pose s’échelonnant de 8 à 10 heures 
selon la nature de échantillon. 
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D’une manière générale on constate que les 
spectres des chlorures sont plus « condensés » que 
les spectres des fluorures correspondants. Cette 
diminution systématique des diamètres des halos 
de diffraction s'explique par l’augmentation géné- 
rale de toutes les distances interatomiques, l’ion F- 
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F1c. 2. — Spectre de diffraction de NaCI fondu à 820 °C. 
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Fi1G. 3. — Spectre de diffraction de KCI fondu à 810 0C. 


étant remplacé par l’ion CI de rayon plus'impor- 
tant. 

Les deux halogénures alcalinoterreux BaCI, et 
CaF, ont été spécialement choisis de manière à ce 
que le rapport des rayons ioniques de l’anion et du 
cation soit sensiblement le même dans les deux c 


ra [rat + = 1,34 PRIT CAT = ae 


On constatera la très grande analogie existant 
entre les deux spectres (fig. 4 et 5), le spectre 
de CaF, étant simplement une version plus 
«dilatée » de celui de BaCI.. 

On notera également le caractère très accusé du 
spectre de LiCI. Comme pour LiF, ceci s'explique 
par le grand écart entre les facteurs de structure du 
cation et de l’anion ; cet effet étant d’ailleurs 
encore plus sensible dans le cas du chlorure. 
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Les spectres de NaCI et de KF enfin présentent 
aussi certaines analogies. 

Les spectres précédentsont été soumis à l’analyse 
radiale de Fourier [6] ; les notations employées sont 
celles des articles précédents. 

Le tableau I rassemble les données numériques 
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Fic. 4. — Spectre de diffraction de BaCI, fondu à 4 000 °C: 
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Fic. 5. — Spectre de diffraction de CaF, fondu à 1 500 °G 
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introduites dans l’analyse de chaque cas particulier, 
notamment : les nombres électroniques effectifs 
moyens X%x et X, du cation et de l’anion ainsi que 
la densité d, du sel fondu à la température £ de 
l'expérience. 

Pour les chlorures ‘alcalins cette densité a été 
calculée d’après les données de Jaeger [7] con- 
firmées par les travaux plus récents de Van Arts- 
dalen et Yaffe [8]. La densité du BaCI, fondu a été : 
empruntée à Drossbach [9], celle du CaF, fondu aux 
travaux de Bääk [10]. 


TABLEAU I 

Kx Ka d °C 
LiCI 2, AT 6 dé20) — 4,49 
NaGlre 410,2 081737 duo — 1,54 
KCI 18,9 0494 déio = 1:50 
BaCl, 58,6 15,7 1000 = 312 
CaF, 24,5 8,25. di5oo — 2:79 


Les figures 6 à 10 montrent les courbes de distri- \ 
bution radiale des halogénures fondus corres- 
pondant aux spectres représentés sur les figures À 
à 5. 
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F1c. 6. — Courbe de distribution radiale de LiCl 
fondu à 620 °C. 


F1G. 7. — Courbe de distribution radiale de NaCl 
< fondu à 820 °C. 


lat Rr'Pm 


Frc. 8. — Courbe de distribution radiale de KCI 
fondu à 810 °C. 


Fic. 9. — Courbe de distribution radiale de BaCl, 
fondu à 1 000 °C. 


F16. 10. — Courbe de distribution radiale de CaF, 
fondu à 1 500 °C. 


III. Interprétation structurale. — En suivant la 
méthode déjà employée dans le cas des fluorures 
alcalins, nous avons repéré sur les courbes de dis- 
tribution radiale les premières distances inter- 
atomiques des composés cristallisés correspondants. 
Ici encore on constate l’accord entre les premières 
distances interatomiques du solide et les distances 
les plus fréquemment rencontrées dans le liquide. 
Ainsi par exemple dans le cas des chlorures alca- 
lins LCI, NaCI, KCI, les courbes comportent deux 
maxima bien résolus ; ils correspondent respecti- 
vement aux deux premières distances inter- 
atomiques : 


X—Cl et X—X (ou Cl— Cl) du cristal. 


La résolution est cependant moins bonne que 
dans le cas des fluorures — en effet il est déjà 
impossible de situer avec précision une troisième 
distance interatomique. 
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En ce qui concerne les composés alcalino- 
terreux BaCI, ou CaF, (fig. 9 et 10), seule la pre- 
mière distance interatomique est résolue ; elle 


JOURNAL DE PHYSIQUE 


N°4 


correspond d’ailleurs à la distance de contact entre 
lP’anion et le cation dans le cristal. 
Dans le tableau IT nous avons établi la compa- 


TABLEAU II 


LiCI 
RONDU CRISTAL 
na(= 74) 4,1 6 
nX(= n4) 9 49 
. BaCl, 


PO 
FONDU CRISTAL (1) 


5,5 8 


(?) Forme « haute ». 


not (— 2n$a) 


raison entre les nombres de coordination relatifs à 
l’état fondu (déduits des courbes de distribution 
radiale) et ceux relatifs à l’état cristallin pour tous 
les composés précédemment étudiés. 

L'étude des chlorures alcalins confirme donc les 
résultats suivants obtenus déjà à propos des fluo- 
rures : 10 Le schéma de coordination du cristal 
persiste à l’état fondu : il y a conservation de 
l’ordre à courte distance. 29 Les nombres de coordi- 
nation à l’état fondu sont sensiblement inférieurs 
aux nombres de coordination à l’état cristallin— il 
y à formation des lacunes. 

Le cas des composés alcalinoterreux est plus 
complexe. Les renseignements sur les distances 
interatomiques à l’état liquide sont par trop incom- 
plets pour qu’on puisse en tirer des conclusions 
certaines. La baisse importante du premier nombre 
de coordination jointe à la diminution de la densité 
à l’état liquide et, en même temps, conservation de 
la distance de contact des ions ne sont pas en 
contradiction formelle avec la formation des lacunes 
à l’état liquide. Toutefois d’autres réarrangements 
structuraux avec modification du nombre de coor- 
dination principal, restent également du domaine 
des possibilités. 

La suite de cet exposé sera exclusivement limitée 
au cas des halogénures alcalins. 


IV. Précisions sur la structure « lacunaire » des 
halogénures alcalins fondus. — En étendant au cas 
des chlorures le mode de raisonnement déjà employé 
pour les fluorures nous pouvons dire que la struc- 
ture des halogénures alcalins fondus est essentiel- 
lement caractérisée par la présence des lacunes ou 


NaCI KCI 
OO RO 
FONDU CRISTAL FONDU CRISTAL 
TA RARE 5,2 6 
9 12 10 12 
CaF, 
FONDU CRISTAL 
n$ (= 2n62) 6,8 8 


vides dans un «réseau » liquide où l’ordre à courte 
distance est le même que celui à l’état solide. Un 
schéma bidimensionnel possible d’une telle struc- 
ture a déjà été indiqué ([1], fig. 7). Les -questions 
qui se posent immédiatement sont celles concer- 
nant la forme, les dimensions et la répartition de ces 
lacunes. S'agit-il de lacunes à l’échelle ionique ou 
bien de véritables cavités de dimensions plus impor- 
tantes ? 

A l’aide de données empruntées à d’autres 
domaines nous essaierons maintenant de déve- 
lopper ce modèle pour le rendre un peu plus quanti- 
tatif. 


19 VOLUME TOTAL DES LACUNES. — [La pro- 


portion (en volume) des lacunes formées dans le. 


sel fondu est donnée par l'expression 
Pres er 
Vr 

qui représente l'accroissement de volume à la 
fusion rapporté au volume total du sel fondu au 
point de fusion. Un travail récent de Schinke et 
Sauerwald [11] renferme une série de mesures 
permettant le calcul de æ pour les différents halo- 
génures précédemment étudiés. Ces valeurs ont été 
rassemblées dans le tableau III. 


Pour les composés en question le part des lacunes 


représente donc 23 % à 15 % environ du volume 
total occupé par le fluide. En ce qui concerne LiF on 
rapprochera le chiffre de 22,7 % obtenu direc- 
tement de celui (21 %) calculé par nous précé- 
demment [1] à partir des mesures du paramètre de 
réseau. L'accord est tout à fait satisfaisant. 


TABLEAU III 


LiF NaF KF LiCI NacCIl KCI 
æ [11] 0,227 0,215 à 147 0,207 0,200 0,147 
L cal/mol 6200 [13] 790013] 6700 [13] 4970 [14] 6560 [15] 6 100 [16] 
tr °C 844 980 856 614 801 768 
o ergs/cm? [17] 255,2 201,6 143,2 440,2 114,1 97,4 
1 ,95 19 1,4 4,55 1,93 


dr g/cem? [7] 1,7 
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20 DiMENSIONS DES LAGUNES. — Pour aborder le 
calcul des dimensions des lacunes il faut disposer 
d’une donnée supplémentaire, par exemple de la 
surface totale des cavités. Dans ce but, on peut être 
tenté d'identifier la chaleur de fusion du composé 
au contenu énergétique dû à la formation des 

cavités, soit : 


L= S5 


formule dans laquelle on désigne par : 

L la chaleur de fusion du cristal ; © l’énergie 
superficielle du liquide au point de fusion ; S la 
_ surface totale des cavités. 

En réalité, l’énergie superficielle effective à 
l’échelle moléculaire doit être inférieure à l'énergie 
. superficielle macroscopique usuelle (Frenkel [12]); 
elle est d’ailleurs impossible à préciser dans l’état 
actuel de nos connaissances. 


Fic. 11. — Schéma de distribution des lacunes dans 
l'hypothèse : 
a) des trous « sphériques » ; 
b) des « fissures fluctuantes ». 


ÉTUDE DES SELS FONDUS PAR DIFFRACTION DES RAYONS X 


AA 


Nous devrons donc nous contenter de la rela- 
tion précédente qui, toutefois, ne doit être regardée 
que comme une première approximation, la 
surface S ainsi obtenue étant une valeur évaluée 
par défaut. 

Pour aller plus loin, il est maintenant indispen- 
sable de faire une hypothèse sur la forme des lacunes. 

-a) Hypothèse des « trous » sphériques. — L’'hypo- 
thèse la plus simple consiste à considérer un ensem- 
ble de P cavités sphériques et de rayon moyen r 


(fig. 11a). 
On dispose alors de trois relations : 


Vr— Vpr = (4/3) rr° P 
SET 
L = So 


r = 8(Vr — Vpr) o/L 


soit, pour un sel fondu de masse molaire M et de 
densité dr, le rayon moyen des cavités est : 


r = 3 x Mo(dr L). 


Pour nous faire une idée du nombre des «trous » 
introduits dans le réseau, nous considérerons 
d’autre part le rapport 


k = PJ2N 


du nombre de trous au nombre total d’ions con- 
tenus dans une molécule-gramme du sel fondu 
(N : nombre d’Avogadro). 

En utilisant les données numériques du 
tableau III, nous avons calculé les valeurs de r et 
de k pour les sels liquides à la température ## de 
fusion. Les résultats ont été rassemblés dans le 
tableau IV. L'examen de ce tableau montre que les 


TABLEAU IV 
LiF NaF KF LiCI NaCI KCI 
r (À) 0,97 0,84 0,69 1,18 0,94 0,82 
k 0,72 1,54 272 0,70 1,80 2,58 
ep (À) 3,22 3,00 3,95 4,35 3,70 4,70 
25 (À) 0,65 0,56 0,46 0,79 0,63 0,55 
m 15,4 8,3 121 19,3 8,6 19,4 


rayons moyens des cavités ainsi définies seraient 
inférieurs aux rayons des ions constituant le réseau. 
(Étant donné les approximations faites, ces chiffres 
ne représentent que des ordres de grandeur.) 
Pour une même série, en allant dans l’ordre 
Li = Na — K on constate cependant que si les 
dimensions des trous ont tendance à diminuer, leur 
nombre, par contre, s’accroît dans de fortes pro- 
portions. 
Pour les sels à cation lourd il y aurait donc pré- 
_ sence d’un très grand nombre de trous de faibles 


dimensions. Ainsi il faudrait admettre qu’il se 
forme en moyenne jusqu’à 2,7 trous par ion consti- 
tutif du « réseau » liquide. Il nous semble difficile 
de concilier une telle proportion de cavités avec les 
résultats fournis par la diffraction des rayons X. 
Il est en effet peu concevable qu’un tel « cortège » 
de cavités associé à chaque ion du réseau n’entraîne 
pas fatalement la disparition de l’ordre à courte 
distance, qui pourtant, comme nous l’avons prouvé 
expérimentalement, est le même que celui dans le 
cristal primitif. 
2 
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Aussi avons-nous été conduits à formuler une 
autre hypothèse dans laquelle le volume des cavités 
est réparti d’une manière différente. 


b) Hypothèse des « fissures fluctuantes ». — Selon 
cette deuxième hypothèse, les cavités ne seraient 
pas uniformément dispersées dans la masse, mais 
au contraire auraient tendance à se grouper en feuil- 
leis ou autrement dit constitueraient des «fissures » 
entre des ensembles ioniques ou « îlots » plus ou 
moins séparés les uns des autres (fig. 11b). L’en- 
semble étant soumis à l'agitation thermique, il y a 
passage des ions d’un îlot à un autre voisin ou, ce 
qui revient au même, déplacement des surfaces de 
séparation — d’où le nom de «fissures fluctuantes » 
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qui nous paraît approprié pour désigner ce phéno- 
mène. Le modèle en question laisse intacts des 
groupements où l’ordre à courte distance caracté- 
ristique du cristal peut être conservé. 

Nous désignerons par e le rayon moyen d’un îlot 
et par à la demi-distance moyenne séparant les 
îlots les uns des autres. En supposant les îlots : 
complétement isolés les uns des autres on aura alors 
les quatre relations suivantes : 


(Vr = (4/3)x(e + 2 
| SH TpaZ 
Léa 

L = S6 


: £ L 
Nina RÉSEAU CONTINU 
s EL A 
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Fic. 12. — Relation entre le nombre de coordination apparent et le degré de morcellement du réseau. 


Z étant le nombre d’ilots par molécule-gramme du 
sel fondu. On en tire : 


Ô = (Vr — Vpr) 5/L 
ou : 
25 — (2J3)r. 


L’écartement entre les îlots est donc du même 
ordre de grandeur que le rayon des «trous » dans 
l'hypothèse a). La valeur de & s’obtient ensuite en 
résolvant l’équation : 


(e + fe — 3 Mo](dr L). 


Une troisième valeur caractéristique sera m, 
nombre moyen d’ions constituant l’ilot : 


m = 2N/]Z. 


Dans le tableau IV nous avons rassemblé les 
valeurs de p, 23 et m calculées à partir des données 
du tableau IIT (ici encore ces nombres ne consti- 
tuent que des ordres de grandeur). 

L'examen du tableau IV montre que les «îlots » 


groupent en moyenne 8-19 ions, c’est-à-dire com- 
portent 1 ou au maximum 2 couches ioniques com- 
plètes entourant un ion. central donné. L’écartement 
entre les différents îlots varie relativement peu avec la 
nature du sel, il est de l’ordre de 0,6 À. On voit que 
les dimensions des îlots sont suffisantes pour assurer 
la conservation de l’ordre à courte distance du 
cristal dans un réseau ainsi disloqué. Quant à 
lorientation des îlots les uns par rapport aux 
autres, elle doit être parfaitement désordonnée car, 
malgré les faibles distances séparant les divers 
groupements (une fraction des rayons ioniques), il 
n’apparaît pas de maxima supplémentaires sur les 
courbes de distribution radiale étudiées. 

Reste à examiner l'influence du morcellement sur 
le nombre de coordination apparent du réseau. Pour 
effectuer ce calcul nous supposerons que les dif- 
férents îlots possèdent la structure du type NaCI, 
qu’ils contiennent tous un même nombre = d’ions 
et qu'ils sont complètement isolés les uns des autres. 

Dans ce cas le nombre de coordination apparent 
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principal de l’ensemble du réseau est identique au 
nombre de coordination apparent principal. N, 
d’un îlot quelconque. N,, est la moyenne pondérée 
de tous les nombres de coordination principaux NV; 
affectés aux m atomes de l’ilot : 
2 ni Ni 
Nn = 
m £ 

Ce calcul a été effectué pour différentes valeurs 
de m,c’est-à-dire pour les îlots comportant un 
nombre croissant de couches concentriques en for- 
mation entourant un ion central donné. La fonction 
(discontinue) N,, — f(m) a été représentée sur la 
figure 12. 

Ceci correspond au cas idéal où tous les îlots 
seraient totalement indépendants. Dans l’hypo- 
thèse où ces îlots seraient partiellement soudés les 
uns aux autres (fig. 11b) le nombre N,, se trouverait 
automatiquement augmenté. La fonction 
Ny = f(m) donne donc la limite inférieure du 
nombre de coordination apparent pour un réseau 
morcelé en îlots identiques mais pouvant par contre 
comporter des « ponts ». Sur le même graphique nous 
avons reporté les points correspondant aux divers 
halogénures étudiés, définis par les coordonnées 
(ra, m). 

On constate que l’essaim de points ainsi définis 
se place bien dans la région du plan délimitée d’une 
part par la droite N = 6 correspondant au réseau 
continu et d’autre part par la fonction N,, — f(m) 
correspondant au cas des «îlots » indépendants. 

Les résultats expérimentaux ne sont donc pas en 
désaccord avec l'hypothèse des « fissures fluc- 
tuantes ». L'examen de la figure 12 permet en outre 
de tirer les conclusions suivantes concernant la 
structure des halogénures fondus : 
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1° Les différents « îlots » constituant le liquide 
sont incomplètement séparés les uns des autres. En 
effet tous les points expérimentaux se trouvent 
nettement au-dessus de la limite N,, — f(m) corres- 
pondant au cas idéal des îlots indépendants. 

(En réalité, pour des raisons techniques les 
nombres #* ont été obtenus pour des températures 
légèrement supérieures aux points de fusion indi- 
qués sur le graphique ; ils représentent sans doute 
des valeurs approchées par défaut). 

29 Les sels de lithium sont constitués par des 
îlots plus gros mais aussi plus indépendants que les 
sels de sodium ou de potassium ; ils se rapprochent 
davantage du cas idéal. 

30 À l’intérieur de chaque groupe caractérisé par 
un même anion la grosseur des îlots est en relation 
avec la température de fusion du composé. 

Ainsi pour les séries : 


LiF (8440) > KF (8560) > NaF (9800) 
et 
LiCI (6140) — KCI (7680) —> NaCI (8010) 


les nombres d’ions « retenus » statistiquement 
autour d’un ion central donné varient respecti- 
vement de 14 à 8 et de 19 à 8, c’est-à-dire que le 
degré de morcellement du réseau augmente avec la 
température de fusion du sel. 


V. Conclusion. — L’étude des halogénures 
fondus par diffraction des rayons X nous a fourni 
un modèle structural qualitatif d’un liquide ionique 
aux températures élevées. Cette image a pu être 
précisée et les résultats rendus plus quantitatifs 
grâce aux données thermodynamiques empruntées 
à d’autres domaines de la physique. 


Manuscrit reçu le 7 janvier 1958. 
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QUELQUES CONSIDÉRATIONS SUR LES AMPLIFICATEURS ÉLECTRONIQUES 
A COEFFICIENT D’AMPLIFICATION NÉGATIF ÉLEVÉ. 
APPLICATIONS AUX APPAREILS DE MESURE. : 


Par G. BLET, 


Résumé. — Nous étudions dans ce travail les amplificateurs électroniques à courant continu 
et à coefficient d'amplification négatif élevé, dans le cas particulier où une boucle de contre-réaction 
relie la sortie à l’entrée. Nous ferons les calculs pour deux montages caractéristiques : à impé- 


dance d’entrée très faible, ou très élevée. 


Dans une seconde partie, nous envisagerons le cas où le circuit d’entrée de l’amplificateur est 
un appareil à spot lumineux, tel que galvanomètre, électromètre, et nous examinerons les réper- 
cussions d’un tel montage sur l’impédance d’entrée, la période propre et l’amortissement. 

Quelques exemples numériques sont joints à chaque cas pour en donner des exemples typiques 


d'utilisation. 


Abstract. — This paper is concerned with D.C.high negative gain feedback electronic amplifiers. 
The calculations are made for two typical cases : those of a very small, and of a very high input 


impedance. 


In a second part, the case is considered of an input consisting in a light-spot instrument, such as 
a galvanometer or an electrometer : the effects of such an input on impedance, period and damping 


are examined. 


Each case is illustrated by some numerical examples. 


“ 


1. Généralités. — Soit un amplificateur à cou- 
rant continu dont le coefficient d'amplification K 
possède une valeur absolue élevée et le signe moins. 
Cela signifie qu’à une tension d’entrée v positive, 
par exemple, correspond une tension de sortie 
V — Kv qui sera négative et | | fois plus grande 
en valeur absolue. Un tel amplificateur doit avoir 
obligatoirement une entrée dissymétrique avec un 
pôle à la masse et de même une sortie dissymétrique 
avec un pôle à la masse. Il existe de nombreux 
schémas pour de tels amplificateurs, de même que 
plusieurs réalisations industrielles. Nous n’insis- 
terons pas sur leur conception qui relève de la 
technique électronique classique et ne présente 
aucune difficulté spéciale. Nous ferons remarquer 
seulement que de tels amplificateurs se comportent 
en courant alternatif comme des déphaseurs : 
en effet, si  — Vo COS ot, la tension de sortie 
V = Ko = Kv, cos wt peut aussi s’écrire : 


V = — |Kl|v, cos ot = |K|v, cos (ot + x). 


Il devient alors indispensable que dans le 
domaine de fréquences utilisé il n’y ait pas de 
rotation de phase. Sifl’on doit travailler à fréquence 
fixe, le problème se simplifie : l’adjonction d’un 
pont déphaseur à la sortie permettra toujours pour 
la fréquence F d’ajuster le déphasage à 1800 exac- 
tement. 


2. Étude de deux types de cireuits. — 2.1. Crr- 
CUIT A IMPÉDANCE D'ENTRÉE RÉDUITE. — (Consi: 
dérons le circuit de la figure 1, où À est un ampli- 


ficateur du type ci-dessus, r la résistance d’entrée 
et À la résistance de contre-réaction. S est la source 
de résistance interne p et de force électromotrice v,. 


F16. 1. — Circuit à impédance d’entrée très faible, 


Soit v la chute de potentiel aux bornes de r. Le 
courant débité par la source S sera : 


i= (va — v)le. (1) 
‘La tension à la sortie de l’amplificateur étant : 
FA EKD (2) 
le courant débité dans la résistance À vaut : 
d = (Y —v)]R (3) 
la chute de tension dans r étant : 
D = it v'}7r, (4) 
De ces équations nous tirons : 
Lo . v + L4 _ Dee 2 
d’où 
2 = [++ au — 0]. (5) 


N°4 


De cette expression rigoureuse, nous pouvons 
tirer une expression approchée exacte, avec assez de 
précision dans le cas où p et r ne sont pas très 
petits devant À, et avec bien entendu |X| > 1, 
hypothèse de base de cette étude. 

Nous aurons : 


(ele) = (vIR) (1 — K). (6) 
Avec la même approximation nous aurons : 
V = —vi{(Rle). (7) 


Remarquons que le quotient v,{ représente le 
courant de court-circuit de la source 1, = Dole ; 
d’où : 

L4 de 1 R. (8) 

Calculons maintenant l’impédance d’entrée appa- 
rente : elle est le quotient de la tension d’entrée v 
par le courant débité 1. 

; = vfi= purs — v) (9) 
SOIL : 

1 


Un + U = AIR ER 


YA = 
Avec les mêmes réserves que précédemment nous 
aurons pour expression approchée : 


ZE ARIK. (11) 


Les deux expressions simplifiées finales [(8) et 
(11)] montrent que la tension de sortie et l’impé- 
dance d’entrée sont proportionnelles à À, résistance 


de contre-réaction. La tension de sortie est indé-. 


pendante du coefficient d'amplification, alors que 
l’impédance d’entrée lui est inversement propor- 
tionnelle. 

En posant Æ' = rJ(R + r), taux de contre- 
réaction inférieur ou égal à l’unité, la relation (10) 


devient, tout en restant rigoureuse : 
D nl) KK'). (12) 


Nous en déduisons également que le courant 1’ 
débité par l’amplificateur est pratiquement égal au 
courant de court-circuit de la source : 


ÿ = in + (AJKK)]. 


2.2. CIRCUIT À IMPÉDANCE D'ENTRÉE ÉLEVÉE. — 
Considérons le circuit de la figure 2, où A est 


Fi. 2. — Circuit à impédance d’entrée très élevée. 


toujours un amplificateur du même type, r la résis- 
tance d'entrée, r’ et À les résistances de contre- 
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réaction. S est la source de résistance interne p et 
de force électromotrice v,. 

Soit v la chute de potentiel dans r et v’ celle 
dans 7’. 

Le courant débité par la source a pour valeur : 


= (D9— 0 + v')fo avec v = ir. (13) 
La tension de sortie V — Xv crée un courant 
D = UK v IR (14) 
et la chute de tension dans r’ vaut : 
= (i+i})r. (15) 
De ces relations on tire : 
; vfr + (Ko — v')[R = v'fr 
SOIL : 
v' = (or fr) (R + Kr)[(R + r'). (16) 


En combinant les équations (13) et (16) on 
obtient : 
Do == o + Lee 
e RO RE TE 
à rapprocher de la relation (5). 
De là on peut déduire : 


(17) 


en (18) 


Cette expression peut se simplifier avec les 
mêmes conditions d’approximation qu’au chapitre 


précédent et donne : 
Fr RIr). (19) 


Nous pouvons calculer l’impédance apparente 
d’entrée par la relation : 


Z = (v — v')]i qui devient 
Æ r'R + Kr 
a 


Sous les mêmes réserves que pour la dernière 
équation, nous obtenons l’équation simplifiée : 


7. 


En posant Æ' = r'J(R + r'), taux de contre- 
réaction, nous aurons si Xr > R : 


(21) 


Z = ri — KK') (22) 

et 
= — v,]K. (23) 
2.3. APPLICATION A UNE PHOTOPILE. — Nous 


allons faire cette application en utilisant ces deux 
circuits pour la mesure (et éventuellement l’enre- 
gistremrent), d’une part du courant de court-circuit, 
d’autre part de la force électromotrice en circuit 
ouvert, d’une photopile au sélénium. 

Soit une photopile de: résistance interne 
e — 10#ohms. Prenons r = 106, r' — 104, R = 104 
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et un amplificateur de POP oCE d'amplification 
K — — 104. 


2.3.1. Mesure du courant de court-circuit. — Dans 
le cas du circuit d’entrée à faible impédance, nous 
trouvons avec les formules simplifiées : 


— —RJK — 1 ohm ... exact à 10% près. 


Le courant 1’ débité par l’amplificateur est égal 
au courant de court-circuit de la cellule : 


FAN 


exact à 10—4 près. 

2.3.2. Mesure de la tension à vide. — Dans le cas 
du circuit à impédance élevée les formules sim- 
plifiées nous donnent : 


Z== r4 — KK) "0,5 102 o0hms. 


La tension de sortie est : 
P= 905. M a2 10 E près: 


2.3.3. Variations de Z avec K et K'. — Il est à 
remarquer que dans le cas de l’entrée à impédance 
très faible, la valeur de celle-ci ne dépend pas de r. 
Seul l’écart entre le courant mesuré 1’ et le courant 
de court-circuit réel à, diminue lorsque r augmente. 
Par contre, dans le montage à impédance d’entrée 
élevée, celle-ci est proportionnelle à r et la tension 
de sortie n’en dépend pas. Il est très intéressant 
qu'avec une résistance d’entrée de 1 mégohm dans 
l’amplificateur on puisse, par une commutation très 
aisée, passer d’une impédance d’entrée de 1 ohm à 
une impédance d’entrée de 5 000 mégohms. L’appa- 
reil de mesure pourra dans les deux cas être un 
même microampèremètre branché en série avec À. 
Par lecture directe il indiquera le courant de court- 
circuit dans un cas, dans l’autre, le produit de son 
indication par À donnera la force électromotrice en 
circuit ouvert. 


Fig. 3. — Variation de l’impédance d’entrée 
en fonction de Æ et K. 


M = Z]r pour impédance élevée 
Mo = Z}r pour impédance faible. 
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Le graphique de la figure 3 montre, à la pré- 
cision des formules simplifiées, comment varie 
l’impédance d’entrée en fonction de X et de K’. 
Nous avons appelé M, et M, les rapports : 


Mc = Zr dans le cas de l’impédance d’entrée élevée. 
M, — Zjr dans le cas de l’impédance d’entrée réduite. 


3. Application au cas où le cireuit d’entrée est 
constitué par un appareil de mesure à spot. — Soit 
un galvanomètre à spot lumineux de sensibilité 
6, radians par ampère. Le montage est celui de la 
figure 4. 


f _D 
PT CT EEE RE c 


FiG. 4. — Schéma optique de l’éclairement du miroir 
et des cellules. 


Un objectif de courte focale f, diaphragmé par . 


une ouverture rectangulaire de surface o = LI forme 
d’une source de surface s une image de surface S à 
une distance D dans le plan où se trouve le miroir 
du galvanomètre de surface s’. Celui-ci forme sur 
une cellule double c, c,, à une distance D', une 
image de surface 6”. 

Soit ® le flux émis par la source. Le flux capté 
par l’ouverture o vaut : 


D = DoJkx f? ; 
il est réparti sur une surface S qui vaut : 
S = sD?]f?. 
Le miroir du galvanomètre en prend la fraction : 
est soit un flux © — Dr LL 
qui se trouve réparti sur l’image 6’ de « en pro- 
duisant un éclairement : 


pp De TRES 
eo D'AMER 2 SUD" 
ou 
= (Ds’Jixs) (1/D'?) 


Pour un courant de : ampères, la rotation du 
cadre sera de 10, radians produisant un dépla- 
cement du spot de 10, D' 
éclairée sera : 


+ 1 0, D’L pour une cellule et 
— i 0, D'L pour l’autre. 


: la variation de surface 


Le à 


No4 


No 4 
La variation différentielle de flux est donc : 
A® = 20, D'LE = ®(s’[2xs) 0, (L[D')i. 


Les valeurs de ZL et D’ ne sont pas arbitraires : 
elles doivent être compatibles avec le diamètre du 
condenseur et la hauteur utile des cellules. 

Nous voyons donc que la sensibilité du montage 
sera d'autant plus grande que : 
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19 la lampe est plus lumineuse ....... ® grand 
20 ]a lampe a un filament ramassé..... s petit 
30 le galvanomètre a un grand miroir. s' grand 
49 le galvanomètre est sensible ....... 6, grand 
5° l’ouverture est une fente allongée .. L grand 
60 la double cellule est proche du galva- 

MOIS PONS NU ee ee D' petit. 


Soient y la sensibilité des cellules en ampères par 
lumen et r, en ohms leur résistance de charge ; un 
courant 2 dans le galvanomètre produira aux 
bornes des cellules une tension différentielle égale à: 


AD = ÿro A0. 
Exemple : Prenons des données tout à fait clas- 
._ siques. 
Pour: la:lampe-......... D — 600 lumens. 
SE rem?, 
Pour le diaphragme... L — 4 cm. 
Pour le galvanomètre... s° — 1 cm?. 
0, = 105 radian/ampère. 
Pour les cellules .:... NN DE=250Cm: 
y —= 3.105 ampère/lumen. 
re = 2.106 ohms. 
Nous trouvons : 
Avr 401 x, 


Soit r, en ohms, la résistance du circuit galvano- 
métrique (y compris celui-ci) ; le coefficient d’am- 
plification de l’ensemble sera : 


KT 10/7. (25) 


Si r — 100 ohms, nous obtenons un gain de 108. 
indépendamment de l’amplification que l’on peut 
mettre en série après les cellules. En général une 
simple double triode à charge cathodique suffira, en 
maintenant le gain à la même valeur, avec une 
impédance de sortie faible. 

Reprenons une expression littérale de X en con- 


servant le montage optique décrit ci-dessus, mais 


SN des | 


en laissant arbitraire la sensibilité du galvano- 
mètre ; nous aurons : 


K = — 1040, fr. (26) 

3.1. MONTAGE A FAIBLE IMPÉDANCE D'ENTRÉE. 
— Dans ce cas, le courant de sortie de l’amplifi- 
cateur est égal au courant de court-circuit de la 
source, avec un écart excessivement faible, et 
l’impédance d’entrée est égale à : 


= — RJK = Rrf104 0,. 
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Rapportée à la résistance du circuit galvano- 
métrique, elle est : 


Zlr = RJ1040, 


avec R =10%éet-0,—="105. 

La résistance apparente n’est que le millionième 
de la résistance vraie. 

La tension disponible V = — 1, R est indépen- 
dante de la sensibilité et des caractéristiques du galva- 
nomètre tant que (v. (5) et (26)) : 


100, > R(T +1). 


Avec 0, — 105 et r — 102, nous voyons qu’une 
source, même de résistance interne très faible, 
comme p — 1 ohm, peut convenir avec, À — 105. 
L’inégalité ci-dessus devient 101 >5>107. La tension 
de sortie est V — — 105, 


Exemple : 

Soit un couple cuivre-constantan de sensibilité 
40 pV/degré et de résistance interne de 4 ohms. 

Pour 10 degré le courant de court-circuit sera : 


do = 40.10 —°/4 = 107$ ampère 
et la tension de sortie 
Pt — 102 volt: 


Cette tension d’un millivolt par millième de 
degré est disponible aux bornes de la résistance de 
son choix ; elle peut être injectée dans un enre- 
gistreur courant : le galvanomètre utilisé est stan- 
dard et robuste. 

Le meilleur galvanomètre utilisé directement 
(un Zernicke, par exemple) ne donnerait qu’une 
déviation de 1 millimètre à 1 mètre et un tel galva- 
nomètre est plutôt fragile. 


3.2. MONTAGE A IMPÉDANCE D'ENTRÉE ÉLEVÉE. 
— Reprenons les équations (22) et (23), l’impé- 
dance d’entrée croît avec la sensibilité du galvans- 
mètre et le taux de contre-réaction. La tension de 
sortie ne dépend que de celui-ci et varie en raison 
inverse. 

Pour X' voisin de l’unité, V est voisin de v, et 
==. 101 


4. Application à un électromètre. — Pour un 
électromètre, le seul montage intéressant est celui 
à impédance élevé. 

La sensibilité d’un électromètre à quadrants 
pouvant. atteindre 10 radians par volt, le montage 
optique étudié précédemment nous donnera un 
coefficient d'amplification À — — 105. 

L’électromètre étant supposé parfaitement isolé, 
son impédance ohmique n’est pas modifié, mais sa 
capacité apparente est réduite, de même que la 
quantité d'électricité demandée à la source pour 
effectuer la mesure, 
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En effet, la tension aux bornes n’étant que 


v = vo(i — KK'), si K' = 5 et K —=.—10, 
p — v,10-4; d’où la charge d'électricité 
demandée : 


7 — C0 CNED 
=. C’ V5. 


La capacité de l’électromètre est divisée par 
10 000. 


5. Dynamique des amplificateurs à gain négatif 
élevé associés à un appareil à spot lumineux. — 
5.1.1. Cas du galvanomètre. — L’équation clas- 
sique du galvanomètre est : 

2 2 

1 (+) += mi (7) 
où l représente le moment d'inertie du cadre ; 
» f le coefficient d'amortissement mécanique ; 
» D, le flux maximum ; 

» R la résistance du circuit ; 

» Cle couple de rappel ; 

» 6] ' l'angle de déviation ; 

»& lecourant traversant l'appareil. 


Quelques relations fondamentales relient ces 


valeurs. 
Si l'amortissement est critique : 


Î +7 2 VIC 


Ts = 2x VIIC 
et 
CO = ®, 1 à l'équilibre. 
Faisons les calculs en u. 6. m. c. g. s. pour le 
galvanomètre standard envisagé précédemment en 
2 


supposant / négligeable devant à Prenons comme 
re 
caractéristiques : À; — 2 000 ohms et 7, = 2,5 


secondes. 
On trouve : 


10$/2 x c. g. ss: 


D, Re Til10:x0, D, 
TREROTÉ I LMO RACE RSA IE T0 L28 RACE 
Cu—="R° TS 1107 /x03 CRM 0T CS 00S; 


Dans le cas du montage à impédance d’entrée 
faible, le circuit de contre-réaction injecte dans le 
galvanomètre un courant : 


TR = KDIR 
Soit 
De KA, 
si le cadre tourne de 0 en posant Æ' = rJ(R +r) 
et, en u. 6. M. c.g.s., 
à — KK’0/100, 
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d’où l’équation représentant les mouvements du 

galvanomètre autour de sa position d’équilibre : 
d?\ dO KK' ®, 

À) & * (c 106% 


d20 
+ (+ )o = 0. (28) 
Do 
RTE 
de 8 vaut : 


— C, donc le nouveau coefficient C” 
CE CA ENR ER) (29) 


Cette nouvelle équation nous montre deux 


résultats essentiels : 


Ne4 


\ 


10 le galvanomètre n’est plus amorti et son, 
régime devient oscillatoire. La nouvelle résistance 


critique devient : 
Re = Rs MCIC- 
20 La période du galvanomètre diminue et 
devient : 
TS ÆNTOV CIC 
La courbe de la figure 5 montre la variation du 


facteur /C[C' en fonction de la valeur de la résis- 
tance À de contre-réaction pour À = — 101 7r. 


EN TRS CNT AB TENTE Te 
Fic. 5. — Variation du couple de rappel en fonction de R 
VCIC’ = 1[V1 —KK. 


Dans le cas du montage à impédance élevée, les » 


résultats sont les mêmes, au coefficient de À près. 
Nous voyons donc que les variations de résis- 
tance critique et de période sont considérables. 


5.1.2. Cas de l’électromètre. — L'équation de 
l’électromètre est : 


d?20 d0 
de + f— + CO = av. 


À dt 


Les équations permettant de déterminer les cons- 
tantes sont les mêmes que précédemment : 


amortissement critique... f — 2 V1IC 
périodes LVL CARRE Ty = 2x VIIC 
SENSIDIN EE NN A REE OJv = x]C 


ou, si 06, est la sensibilité en radians par volt : 
«f[C = 0, 10—$ en radians par u. é m. 
Prenons un bon électromètre Curie : 


Ty, = 20 secondes  N 


EUR 


N° 4 


D'où : 
He A107 JT 
Co "107 
HEAR) LE 


Soit v la tension correspondant à une rotation 0 : 
la tension réinjectée vaut vKK', d’où : 
_d20 ,dô ; 
IS HIS +CU—KK)0—- 0 
19 L’amortissement diminue et le régime devient 
oscillatoire : mais là, il n’y a pas de résistance cri- 
tique permettant de retrouver l'amortissement cri- 
tique. 
20 La période est diminuée. 


VCIC' = 1JV1 — KK’ 


Les variations en fonction de Æ' sont données 
par la courbe de la figure 6, pour À — — 105. 


(même formule que (29)). (31) 


10 10? 105 10 10$ 106 
Fic. 6. — Variation du couple de rappel en fonction de 
AK = (R + r)fr'.VCIC = 1]V1 — AK. 


La période de l’électromètre peut donc être 
réduite considérablement. 


5.2. INTRODUCTION D'UN TERME D'AMORTIS- 
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Si XK' est beaucoup plus grand que l’unité, on 
peut écrire CXK' = C'et négliger D /R,, d’où : 


Ry=2VIIC = TIr (32) 


relation qui indique que l’ensemble Ry a une cons- 
tante de temps de l’ordre de la nouvelle période. 
Si la source a une résistance p beaucoup plus 
faible que À, le terme ®%/o n’est plus négligeable 
et la nouvelle valeur y’ est plus petite que Y. 


9.2.2. Eleciromètre. — Le calcul est semblable : 
| f2 RCKK" y = 9 VTC. 
Avec les mêmes remarques que ci-dessus : 
Te 
RyY = Fe 
On peut noter que dans ces conditions limpé- 
dance du condensateur à la fréquence des oscil- 
lations vaut : 
Yo = RI/2: 
Exemple : 


Soit un électromètre Curie de faible sensibilité 
avec 6 volts seulement sur l'aiguille : 0, — 0,25. 
Prenons 
R—=106 et r» = 5.104. 
Nous avons 


RE LORS TO 2 


SEMENT PAR CONTRE-RÉACTION CAPACITIVE. 
Supposons dans tous les cas que le taux de contre- 
réaction n’est pas très élevé, c’est-à-dire que r et r’ 
sont petits devant À. Plaçons aux bornes de À un 
condensateur de capacité y. Le courant qui le tra- 
| vérse vaut : % 

; dv 

LASER 4 du: 

Le second membre de l’équation de l’appareil de 

mesure devient : 


+ COKK À RCKK'ydOfdt. 
La période reste ce que la contre-réaction 
ohmique l’a faite, mais l’amortissement est aug- 
menté. 
5.2.1. Galvanomètre : 
(D2 Re) — RCKK' y = 2 VIC’. 
Il est possible par le choix de y de conserver 


l'amortissement critique avec la même résis- 
tance R'et la nouvelle période T. 


C 


L 


Tt — 


10 s 
Fic. 7. — Enregistrement de la déviation du spot d’un 
électromètre en fonction du temps : 
a) sans contre-réaction, 


b) avec contre-réaction À, , 
c) avec contre-réaction R et amortissement C. 
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d’où : 


KK'=.-—41,25.102 et VA —KK 11. 

Donc la nouvelle période, au lieu de 20secondes, 
sera de 1,8 seconde. L’électromètre moins amorti 
est oscillant. Pour le mettre à l’amortissement cri- 
tique, il faut mettre aux bornes de R une capacité 
de 0,55 microfarad. 

Les trois courbes de la figure 7 montrent de 
façon très nette les trois stades de fonctionnement. 

Ces trois courbes représentent l’application et la 
suppression d’une tension sur l’électromètre : 


So T; — 20 secondes 
T — 1,8 seconde 
régime oscillaioire 


a) sans contre-réaction 
b) avec contre-réaction R ... 


c) avec contre-réaction À et 
amortissement C....... T —1,8 secondes 


régime amortr. 
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N. B. — Les oscillations subsistantes sont dues 
au fait que l’électromètre en essai n’était malheu- 
reusement pas soustrait aux vibrations du labo- 
ratoire. | 


* 


Conclusions, 


Il ressort de cette étude que l’emploi des ampli- 
ficateurs à coefficient d'amplification négatif élevé 
et à contre-réaction apporte une très grande sou- 
plesse à la résolution des problèmes d'amplification. 

Dans le cas d’une résistance d’entrée, on peut 
agir à volonté sur l’impédance d’entrée. 

Dans le cas d’un appareil à spot lumineux, on a 
en outre le choix de la période et de l’amortis- 
sement. Il faut noter de plus que dans ce cas l’appa- 
reil de mesure à spot fonctionne pratiquement en 
appareil de zéro et ses performances n’influent pas 
sur la précision des mesures. 

Manuscrit reçu le 16 décembre 1957, 
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LE CYCLOTRON DE L'ÉCOLE SUPÉRIEURE DE PHYSIQUE ET DE CHIMIE DE PARIS 


Par M. G. DUPUY, 


Chef de Travaux. 


Résumé. — L’auteur décrit en détail un petit cyclotron dont le projet et la réalisation sont dus 
à un groupe d'élèves de l’École Supérieure de Physique et de Chimie de Paris. Cet instrument est 


utilisé à des fins pédagogiques. 


Abstract. — The author describes in details a small cyclotron which has been designed and 


constructed by some students of the Ecole de Physique et Chimie in Paris. 


used for educational purposes. 


. Nous avons entrepris en 1953 la construction 
l'un petit cyclotron pour les besoins de l’ensei- 
mement au laboratoire de Physique atomique de 
"École. 

L'étude et la mise au point de l’appareil ont été 
ffectuées par des groupes d'élèves de quatrième 
innée des promotions successives dont le temps de 
séjour au laboratoire est de vingt semaines par 
année scolaire. Ils ont eu à mettre en œuvre leurs 
‘onnaissances dans les divers domaines de la Phy- 
ique : magnétisme, haute fréquence, technique du 
ide, physique nucléaire, etc. 

Les dimensions de l’appareil sont restreintes 
urtout pour des raisons budgétaires. M. R. Lucas, 
Directeur de l’École, a bien voulu nous aider et 
ous encourager ; qu’il veuille bien trouver ici nos 
emerciements pour l’appui et l’intérêt qu’il n’a 
essés de manifester au cours de cette réalisation. 

La construction a été presque totalement effectuée 
vec'les moyens de l’École grâce à M. J. Luc, colla 
orateur technique, que nous nous voulons remer- 
ier tout particulièrement pour sa contribution et 
es suggestions. 

Les performances sont évidemment en rapport 
vec les dimensions du cyclotron. Actuellement, de 
açon courante on obtient un flux de protons dont 
énergie est de 600 000 électron-volts et dont 
intensité est de 6 microampères ; Ce qui permet 
effectuer un très grand nombre de réactions 
ucléaires et d’investigations concernant l’émission 
es particules. Par la modification de la ligne haute 
‘équence, qui intentionnellement a été réglée pri- 
itivement sur une fréquence de 12 mégacycles, 

sera possible d’élever l’énergie des protons 
isqu’à 1 million d’électron-volts environ. 


Principe du fonctionnement. — Je rappellerai 
ès brièvement le principe de fonctionnement de 
t accélérateur dû à Lawrence. 

Une chambre étanche (fig. 1), dans laquelle on 
it régner un vide élevé, contient deux électrodes 
‘euses, semi-cylindriques, appelées dés. Elles sont 


The instrument is 


disposées symétriquement et forment approxi- 
mativement un cylindre, les bords rectilignes des 
dés étant séparés par un petit intervalle. La 
chambre est placée entre deux pôles d’un électro- 
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FiG. 1. — Schéma de la chambre du cyclotron. 


amant dont l’axe de révolution passe par le centre 
de symétrie des dés. Le champ magnétique d’in- 
duction, supposé uniforme, courbera la trajectoire 
des particules chargées issues d’une source d’ions, 
localisée entre les dés, au voisinage immédiat du 
centre de symétrie. Les dés sont soumis à une diffé- 
rence de potentiel alternative de fréquence conve- 
nable : soit v = V, sin 2r ft. 

Les ions, dont la charge positive est un multiple 
entier Z de la charge élémentaire e, sont accélérés, à 
partir du lieu de leur production, vers le dé qui est 
négatif à l'instant considéré ; en même temps le 
champ magnétique d’induction B courbe les tra- 
jectoires. Pour simplifier, nous supposerons que le 
champ électrique alternatif a toutes ses lignes de 
force parallèles au plan de symétrie des électrodes 
et perpendiculaires aux bords rectilignes des dés. 
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Dans cette hypothèse, les ions ont en partant de la 
source, une trajectoire qui se confond avec une des 
lignes de champ électrique en l'absence du 
champ B. Dès qu’une particule pénètre à l’intérieur 
d’un dé, celui-ci soustrait l’ion à toute action élec- 
trique, si on néglige les forces provenant des autres 
ions présents dans l’électrode, car il sert alors de 
cage de Faraday. 

A l’intérieur, sous l’influence de B, la trajectoire 
d’une particule de masse FM est circulaire et par- 
courue à vitesse constante v telle que 


Mv?fr = ZevB 
on voit que la vitesse angulaire de rotation 
wo = vfr = ZeB]M (4) 


est invariable à condition toutefois que B, M, Z 
restent constants. 

Si, pendant le temps que la particule met pour 
décrire une demi-circonférence, la différence de 
potentiel alternative a changé de sens, l’ion sera de 
nouveau accéléré à son passage entre les dés. Dési- 
. gnons par V la valeur moyenne de la tension inter- 
électrodes relative à la durée de passage ; l’énergie 
de la particule s'accroît de ZeV pendant la tra- 
versée et sa vitesse d’une quantité Av telle que 


ZeV = (1/2) M A v?. 


Le rayon de courbure de la trajectoire sous 
l’autre électrode augmentera à son tour pour que 
l'équation (1) reste vérifiée. Le processus conti- 
nuera tant que les ions se présenteront entre Les dés 
à des instants favorables à l’accélération. Schéma- 
tiquement, la trajectoire sera constituée à l’inté- 


rieur des électrodes par un ensemble de demi-. 


circonférences dont les rayons croissent réguliè- 
rement et quisont reliées tangentiellement à travers 
l’espace qui sépare les dés. Il y aura synchronisme, 
c’est-à-dire possibilité d’accélérations successives 
d’un ion, lorsque la demi-période 1 /2f, de la diffé- 
rence de potentiel appliquée entre les électrodes, 
sera égale au temps £ qui s’écoule entre l’entrée et 
la sortie d’une particule dans un même dé ; or, 
comme nous avons 


t est indépendant du rayon ; c’est un invariant 
avec les hypothèses faites. 
La condition de synchronisme, dite encore de 


résonance, est par suite 


Es ZeB 
7 9x M° 


Quand elle est réalisée la particule quitte 
toujours un dé au moment où la tension accélé- 
ratrice reprend la même valeur V. Après un certain 
nombre de révolutions la trajectoire a un rayon r, 
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les ions ont subi n accélérations successives et pos- 
sèdent une énergie cinétique 


DA ZAe Tr 


1 
= — D Es 
ste M 


= 2n°r°f2 M (2) 


ou encore puisque 
(1/2) Mv? = nZeV 
r = V2MV/]ZeBi.Vn. 


L'énergie des particules dépend donc seulement 
des valeurs de B? et de r°. La loi de variation des 
rayons des trajectoires s’obtient à partir de (2) 


2Arfr = AEJE. 


Le Ar relatif à deux spires voisines à l’intérieur 
de la même électrode est, puisque dans ces ondes 
tions AE = 2 ZeV, 


Ar = r.ZeV JE: ; (3) 


Lorsque les particules accélérées sont des pro- 
tons la fréquence de résonance exprimée en méga- 
cycles par seconde est 


fn = 1,52 B 


B étant le champ d’induction en kilogauss. 
On aurait avec les mêmes unités dans le cas des 
deutons 


fa = 0,76 B. 


L'énergie des ions en MeV est calculable faei- 
lement pour les protons et les deutons au moyen 
des formules 


En = 0,28 B?r?10—1 
Eg = 0,24 B?r?10—4 


B étant toujours en kilogauss et r, rayon de la 
trajectoire finale, en centimètres. | 


Domaine relativiste. — Les choses deviennent 
moins simples quand les vitesses atteintes par les 
ions arrivent dans le domaine de correction relati- 
viste ; la masse m des particules n’est plus cons: 
tante, elle croît à chaque accélération. La condition 
de synchronisme, le champ d’induction étant 
toujours supposé constant, n’est plus satisfaite. 
Sur une trajectoire de rayon r la vitesse est infé- 
rieure à celle qu’on aurait à la résonance ; la parti- 
cule arrive avec un retard de phase dans l’espace 
accélérateur ; les retards  s’accumulent et 
s’amplifient avec l'énergie acquise jusqu’au 
moment où lion arrive entre les dés quand la diffé: 
rence de potentiel est nulle. Cette cause relativiste 
de déphasage limite l’énergie que l’on peut commu 
niquer aux particules par ce procédé. 

Pour compenser ces retards de phase il semble 
qu’il suffise, en donnant une forme convenable aux 
pièces polaires, de s’arranger pour que B croisse en 
même temps que r afin de maintenir constant le 
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ZeB 
2x M° 
configuration radiale du champ va à l’encontre 
d’une condition impérative, celle d’obtenir .une 
focalisation du faisceau d’ions dans le plan médian 
des dés pour ne pas en perdre une fraction trop 
importante. Pour cela il convient que B soit au 
contraire une fonction décroissante du rayon. 

Ces deux conditions contradictoires font qu’en 
pratique l’énergie est limitée, pour des protons, 
à 30 MeV environ, bien qu’on ait cherché à res- 
treindre le déphasage total en diminuant le plus 
possible le nombre des accélérations par lappli- 
cation entre les dés de tensions très élevées tout en 
évitant la formation de décharges disruptives. 


rapport f — Malheureusement une telle 


Focalisation électrique et magnétique. — Les 
études théoriques dé Rose et Wilson ont montré 
que les trajectoires des ions étaient compliquées par 
des oscillations verticales. Nous avions supposé 
initialement que les lignes de force du champ élec- 
trique étaient rectilignes et parallèles au plan de 
symétrie horizontal des électrodes. Elles possèdent, 
en réalité, une certaine courbure, comme l’indique 
la figure 2. La composante verticale du champ agit 


Jer Dé Lignes de force 20 


Trajectoire de la particule 


F1G. 2. — Focalisation électrique. 


alors sur les particules dont la vitesse n’est pas dans 
le plan médian, qui est un plan de symétrie pour le 
champ. Un grand nombre de particules ont une 
vitesse oblique par rapport à ce plan ; aux faibles 
énergies par exemple, les chocs contre les molécules 
du gaz résiduel modifient la direction des vitesses. 
De toutes façons la composante verticale de la 
force électrique les entraîne plus ou moins vite vers 
les parois des dés où elles sont neutralisées. L’inter- 
valle accélérateur fonctionne comme une lentille 
ionique. L’examen des lignes de champ électrique 
montre qu’une particule circulant à l’intérieur d’un 
dé, en s’écartant du plan médian, sera d’abord 
ramenée vers lui, puis s’en écartera. Ces deux effets 
ne se compensent pas exactement parce que 
l'énergie de la particule augmente et que le champ 
varie pendant la durée de passage. 


Qualitativement on peut dire que le temps de. 


séjour dans la partie divergente du champ est 
inférieur à celui passé dans la partie convergente 
puisque la vitesse de l’ion augmente, de sorte que 
des deux déplacements verticaux, le premier 
l'emporte ; il y a focalisation. L’effet est encore 
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plus important, comme le montre le calcul, lorsque 
la particule franchit l'intervalle pendant la décrois- 
sance de la tension accélératrice. 

La focalisation par le champ électrique n’a lieu 
qu’au voisinage du centre, dans la région où les 
ions n’ont pas encore une grande énergie ; elle peut 
être améliorée à plus grande distance grâce au 
champ magnétique d’induction. Dès que l’on 
s’approche de la périphérie des pièces polaires, les 
lignes de force ont une certaine courbure, la conca- 
vité est tournée du côté du centre de symétrie. La 
force de Lorentz d’origine magnétique agissant 
sur un 10n en mouvement en dehors du plan médian 
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Fic. 3. — Focalisation magnétique. F force de Lorentz 
agissant sur l’ion I dont la vitesse est perpendiculaire au 
plan de la figure. (B champ d’induction magnétique en I.) 


est dirigée vers ce plan (fig. 3). La décroissance 
radiale du champ, qui donne naissance à cet effet 
de focalisation, peut être provoquée ou amoindrie 
en modifiant l’entrefer. La mise en forme du champ 
se fait de façon empirique au moyen de couronne 
ou de disques de fer placés de part et d’autre de la 


‘chambre du cyclotron. La diminution du champ 


doit être faible afin de rester au voisinage de la 
résonance, l’ordre de grandeur étant de 1 à 2 % 
pour les raisons suivantes. La focalisation électrique 
étant bientot inefficace à une certaine distance du 
centre, il est avantageux de disposer l’orifice de la 
source d’ions au voisinage immédiat du plan 
médian pour avoir le plus grand nombre possible 
d’ions dont la vitesse initiale soit dans ce plan; 
de plus on aura une intensité du faisceau d’ions 
plus grande en aidant la focalisation électrique par 
une focalisation magnétique obtenue en faisant 
décroître le champ d’induction, dès le centre, d’une 
valeur supérieure à celle de la résonance pour 
arriver finalement à une valeur inférieure. 

La décroissance imposée au champ, le long d’un 
rayon, à pour conséquence une variation de la 
phase de la tension accélératrice correspondant à 
l’arrivée des particules entre les dés. 

Puisque le champ au centre est supérieur à celui 
fixé par la condition de résonance, les ions arrivent 
en avance dans l’intervalle accélérateur ; le con- 
traire a lieu, plus tard, quand les ions arrivent dans 
la région où le champ d’induction magnétique est 
devenu inférieur au champ de résonance. 

Comme la phase initiale, ainsi qu’il a été dit, 
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doit être telle que le champ électrique aille en 
décroissant pendant la traversée de l'intervalle 
entre les électrodes, les particules dès les premières 
accélérations, arriveront entre les dés avec une 
avance croissante avec le nombre des passages. Il 
ne faut pas que cette accumulation d’avances les 
fasse arriver à l’instant où la tension appliquée 
entre les électrodes devient nulle ou négative. Il 
convient donc que les ions pénètrent dans la zone 
où le champ est inférieur à celui de la résonance 
avant cette éventualité. À partir de ce moment ils 
prendront régulièrement du retard à chaque pas- 
sage. Le faisceau d’ions doit être utilisé, quand par 
suite des retards successifs, ils passent entre les 
dés alors que la tension passe par un maximum ; 
c’est dans ces conditions que la variation du rayon 
de la trajectoire entre deux traversées est la plus 
grande. Ainsi une fraction minime du faisceau ren- 
contre la paroi verticale du dé quand il se présente 
à la fenêtre de sortie. 


Description des divers éléments 


L'’électro-aimant. — Pour des raisons d'économie 
nous avons utilisé la carcasse d’un alternateur 
formée de deux demi-anneaux séparables de 104cm 
de diamètre intérieur, de 41 cm de largeur et 
de 10 cm d’épaisseur, dont le poids est de 1 200 kg 
environ. Ce sont ces dimensions qui ont été impé- 
ratives dans la suite du projet et qui limitent le 
champ d’induction qu’on peut obtenir dans l’entre- 
fer. Le flux de fuites est en effet plus important dans 
un circuit magnétique de cette forme par rapport à 
celui qu’on aurait dans un circuit à section rectan- 
gulaire où les bobines d’excitation peuvent occuper 
une grande partie du volume disponible. 

Les deux pièces polaires sont des cylindres 
d’acier doux de 50 em de longueur et de 30,2 cm de 
diamètre, l’entrefer a été amené à la valeur fixée en 
plaçant entre les demi-anneaux des cales d’épais- 
seur convenable ; elles ont été retouchées jusqu’à 
obtenir un entrefer de hauteur constante à 1/10 
de mm près. Nous avons choisi, pour avoir une 
bonne marge de sécurité quant à la tenue à l’élé- 
vation de température des bobines d’excitation 
(9 par pôles), du fil de cuivre de 1,3 mm de diamètre 
entouré par de la soie de verre imprégnée d’araldite 
qui garde ses qualités isolantes jusque vers 1800. Les 
faces de chaque bobine sont en contact avec des 
plateaux de cuivre par l’intermédiaire de disques 
minces de mécanite, le même matériau les isole des 
pièces polaires ; on a soudé sur la périphérie des 
plateaux des tubes de cuivre parcourus par un 
courant d’eau qui assure un refroidissement efficace. 

La génératrice à courant continu qui alimente les 
bobines de champ débite au plus 32 ampères sous 
220 volts. Le courant d’excitation de cette dynamo 
est fourni par un groupe indépendant. La mise en 
service de l’ensemble est effectuée à partir d’une 
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cabine comportant les relais temporisés habituels. 
La chambre d’accélération a une hauteur intérieure 
de 70mm. Le champ d’induction est constant dans 
l’entrefer à l’intérieur d’un cercle de 22 cm de dia- 
mètre, il décroît ensuite de 8 % quand on se 
trouve à 14,5 cm de l’axe des pôles, le champ au 
centre étant de 8 000 gauss. La décroissance radiale 
a été ramenée à 1 % environ par correction au 
moyen de plusieurs couronnes de tôle de diamètres 
intérieurs différents empilées dans l’entrefer. 

Dans l’état actuel l’assujettissement du champ à 
la condition de résonance se fait manuellement. 
L'élévation de température des bobines d’exci- 
tation entraîne une variation de leur résistance ; 
l'équilibre de température est atteint en deux 
heures. Le champ est alors stable ; dans la période 
transitoire on corrige la dérive du champ en agis- 
sant sur le courant d’excitation de la génératrice 
au moyen d’un rhéostat. 


Mesure du champ magnétique. — A côté des 
procédés classiques une des meilleures méthodes 
pour mesurer un champ magnétique d’induction à 
mieux que 107 4 près consiste à utiliser la résonance 
magnétique nucléaire. En principe les moments 
magnétiques des atomes s’orientent dans le champ 
magnétique, on provoque des transitions entre les 
différents états quantiques au moyen de champs 
magnétiques variables en haute fréquence de phase 
et de polarisation convenables. La condition de 
résonance magnétique nucléaire est 


ky = x 


u étant le moment magnétique du noyau, I le spin 
du noyau, v la fréquence de résonance, A le champ 
d’induction magnétique. 

Pour les protons 1 — 1/2 et 

u = (1,4100 + 0,0002)107?3 dyne/cem/gauss. 
La fréquence de résonance est alors 42,6 méga 
cycles dans un champ de 10 000 gauss. 

Quand on veut mesurer des champs plus élevés 
on utilise une solution de chlorure de lithium afim 
d’avoir des fréquences de résonance plus petites. 
(lu = 3/2: ur = 1,6443 107? dyne/em/gauss) 

Une condition essentielle est que l’nhomogénéité, 
du champ ne soit pas trop grande. La variation du 
champ dans le volume du corps utilisé ne doit pas 
en effet être supérieure à 10 gauss au maximum 
La sonde est constituée par un petite tube conte” 
nant environ un centimètre cube d’eau tenant 
solution à faible concentration un sel paramagnés 
tique : chlorure ferrique ou sel de manganèse. 
tube est placé dans une petite self dont l’axe est 
perpendiculaire au champ d’induction que Po 
désire mesurer. Elle est alimentée par le cour 
haute fréquence d’un générateur étalonné. ‘4 

L'effet Hall peut servir également à mesurer 1 


+ 
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champ magnétique en utilisant un cristal de ger- 
manium généralement du type n ou mieux un 
cristal d’antimoniure d’indium. Taillé en forme de 
parallélépipède rectangle de 5 mm de longueur 
environ et placé perpendiculairement aux lignes du 
champ, il permet la mesure du champ local à 1 % 
près. Le coefficient de l'effet Hall est sous la dépen- 
dance de la température. Une variation de 109 
entraîne une erreur de 4 % avec le germanium ; de 
même la résistivité du germanium est fonction du 
champ magnétique. 


La chambre étanche et les dés (fig. 1). — La 
chambre est un anneau cylindrique en bronze 
de 360 mm de diamètre intérieur et de 460 mm de 
diamètre extérieur. La hauteur disponible à l’inté- 
rieur de l’enceinte est de 70 mm. Elle est fermée par 
deux couvercles d’acier doux dont l’épaisseur est 
de 20 mm jusqu'aux points situés à 180 mm du 
centre ; l’épaisseur varie brusquement alors et 
passe à 18 mm. C’est dans cette deuxième partie 
des deux couvercles que sont logés les joints 
toriques d’étanchéité en caoutchouc synthétique. 
Elle est munie d’un certain nombre d'ouvertures 
circulaires permettant la liaison avec les pompes à 
vide, le passage des conducteurs haute fréquence 
portant les dés et celui de l’alimentation de la 
source d’ions, de l’électrode déflectrice, etc. Toutes 
ces ouvertures sont rendues étanches par l’emploi 
de joints toriques. Les dés, portés à l'extrémité des 
conducteurs haute fréquence, peuvent être déplacés 
de l'extérieur au cours du fonctionnement du cyclo- 
tron. Dans ce but chacun des conducteurs hori- 
zontaux est lié à une tige isolante verticale en 
stéatite dont l’autre extrémité est solidaire d’une 
membrane en tombac compressible ou extensible 
de façon dissymétrique au moyen de trois vis de 
réglage. Les variations de position des dés influent 
sur l'intensité du faisceau d’ions reçu par une cible, 
qui peut être amené ainsi à sa valeur optima. Les dés 
sont refroidis par contact avec la ligne haute fré- 
quence qui les supporte et qui, elle, est parcourue 
par un courant d’eau. 


L’installation à vide. — Le volume à vider est de 
l’ordre de 20 litres. L'installation comporte un 
groupe de pompage : pompe à palettes, pompe à 
diffusion de vapeur d’huile, dont le débit avec 
robinet écran est de 150 litres seconde à 1074 mm 
de mercure. Par un jeu de robinets, on peut faire le 
vide dans l'installation avec la pompe à palettes 
seule. Il a été fait usage pour réaliser les canali- 
sations de tubes de cuivre qui servaient d’anodes 
pour des lampes de puissance et qui par consé- 
quent ne présentaient pas de porosité. Tous les 
joints sont toriques en caoutchouc synthétique, ils 
sont écrasés dans des gorges par serrage. Le degré 
de vide est mesuré avec une jauge à ionisation, la 
pression obtenue est inférieure à 1075 mm de 
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mercure. En fonctionnement lorsque de l’hydro- 
gène est envoyé dans la source d’ions la pression est 
de l’ordre de 1074 mm. Un tube à décharge est 
également utilisé pour apprécier le vide obtenu 
avant de mettre dans le circuit la pompe à diffusion. 


La source d’ions (fig. 4). — Après divers essais 
nous avons utilisé le type de source à décharge à 
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basse tension. Elle est constituée par une boîte 
plate, en cuivre, refroidie sur sa périphérie par de 
l’eau circulant dans un petit tube de cuivre. La 
face supérieure de la boîte porte un petit tube de 
cuivre amovible disposé parallèlement aux lignes de 
champ magnétique ; son extrémité est taillée obli- 
quement. Il contient un morceau de graphite, percé 
suivant l’axe. Il occupe une fraction de la hauteur. 
Cette source est quasi ponctuelle ; le diamètre du 
canal est de 1,6 mm. A l’intérieur se trouve un 
filament de tungstène chauffé par un courant de 
haute fréquence. Le gaz (hydrogène, deutérium ou 
hélium) arrive dans la boîte par un petit tube de 
cuivre. Le débit est réglé par une fuite variable à 
volonté qui est un robinet à pointeau d’un modèle 
commercial. Les électrons sont accélérés par une 
différence de potentiel continue de 150 volts envi- 
ron appliquée entre le filament et la boîte ; ils 
ionisent le gaz par choc. 

Le courant atteint plusieurs centaines de milli- 
ampères quand la décharge d’arc est amorcée. Le 
champ magnétique agit sur le faisceau d’électrons 
et d’ions pour canaliser verticalement les trajec- 
toires suivant les lignes de force. La décharge tra- 
verse un canal pratiqué dans le graphite ; à l’inté- 
rieur de la source la pression qui règne est de l’ordre 
de 1/100 mm de mercure. Le gaz s’écoule par le 
canal ; on fixe ses dimensions, qui dépendent de la 
vitesse de pompage, de sorte que la pression rési- 
duelle dans la chambre du cyclotron soit seulement 
de 1 à 2.107 4 mm de mercure ; l’idéal étant d’avoir 
une Pression élevée dans la source et inférieure 
à 10—4 mm dans la chambre. La tension haute fré- 
quence de chauffage du filament, provenant de 
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deux triodes oscillatrices montées symétriquement, 
est amenée depuis la cabine par une ligne coaxiale 
dont l’impédance est adaptée. La puissance dis- 
ponible est variable à volonté en agissant sur le pri- 
maire du transformateur qui alimente les lampes 
redresseuses. É 

En marche normale la puissance fournie au fila- 
ment est de 70 watts. À la sortie de la source on 
recueille un courant d’ions dont l'intensité est de 
l’ordre de quelques milliampères ; il est formé de 
protons H+ et de molécules une fois ionisées Hf. 
L’intensité du courant reçu sur la cible, placée sur 
les trajectoires des ions qui ont franchi l’ouverture 
de sortie pratiquée dans le dé, atteint 6 micro- 
ampères. Cette valeur n’est pas la valeur limite ; 
elle a été augmentée jusqu’à 10 microampères en 
chauffant davantage le filament. En maintenant 
intensité à 6 microampères le filament a déjà 
fonctionné plus de 300 heures sans incident. 
L’orifice de sortie de la source se trouve approxi- 
mativement dans le plan de symétrie horizontal des 
dés, les particules susceptibles de subir le cycle 
complet des accélérations étant celles qui pénètrent 
dans les cavités des électrodes au voisinage de ce 
plan. Des plaquettes de cuivre, fixées obliquement 
par rapport aux faces horizontales des électrodes, 
viennent près de l’orifice et augmentent l’émission 
utile de la source. La source d’ions est déplaçable 
de l'extérieur latéralement et verticalement pen- 
dant le fonctionnement ce qui permet de rechercher 


la position donnant le courant le plus intense sur la 
cible. 


Générateur à haute fréquence (fig. 5). — Le 
montage choisi est du type oscillateur à circuit 


Fic. 5. — Schéma de l’oscillateur haute fréquence. 


grille accordé ; les dés et les deux tubes de cuivre 
qui les supportent à une extrémité et qui sont 
réunis à l’autre par une self de quelques spires, 
constituent le cercuit anodique ; il est réglé pour 
qu’il fonctionne en quart d’onde. 

: Le générateur comprend deux triodes montées 
symétriquement ; ce sont des 3.T.4 000 refroidies 
chacune par un ventilateur placé sur le support. 
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Elles sont susceptibles de dissiper une puissance 
maximum de 4 kilowatts. La ligne haute fréquence 
est ouverte ; les pertes par rayonnement sont plus 
importantes qu'avec une ligne fermée, mais le mon- 
tage est plus simple, bien qu’il soit nécessaire 
d'isoler les deux lignes par rapport à la chambre du 
cyclotron, ce qui ne présente pas de difficultés 
particulières : les tensions utilisées ne dépassantpas 
quelques dizaines de milliers de volts. 

L’isolement a été réalisé à partir de 4 lampes 
oscillatrices de puissance hors d’usage dont les 
anodes extérieures en cuivre faisaient corps avec 
une couronne de ferro-chrome soudée elle-même au 
ballon de verre. É | 

Chacun des ballons à été divisé en deux pour 
éliminer le culot, puis, après avoir coupé une partie 
des anodes, les demi-ballons ont été associés par 
paire et resoudés. On a obtenu ainsi un ballon sur 
lequel sont soudés. 2 tubes aux extrémités d’un dia- 
mètre. Sur les portions d’anodes conservées oh a 
brasé des couronnes de laiton assez épaisses, com- 
portant un évidement pour loger un joint torique. 
La jonction avec la ligne haute fréquence est alors 
aisément effectuée. 

L'alimentation anodique du type classique est 
filtrée. Lorsque la puissance demandée par la ligne 
devient trop grande, par exemple lors d’un court- 
circuit entre les dés et la masse au cours de réglage, 
un disjoncteur entre en action et coupe la tension 
anodique ; des milliampèremètres et voltmètres 
indiquent la tension anodique, le courant anodique 
moyen et le courant grille. 


L'accord du circuit grille se fait par condensateur M 


variable. La fréquence du circuit anodique peut 
être modifiée en faisant varier la longueur de la 
ligne par coulissement des tubes télescopiques ou 
en changeant le nombre de spires de la self placée en 
tête. 


Initialement le circuit a été réglé sur une lon-« 


gueur d’onde de 24,9 mètres soit 12,04 mégacyles: 
Le champ magnétiqué' de résonance correspondant 
est de 7 973 gauss. 


Étude du montage. — L'étude de l’oscillateur. 


4 


EL 


peut être conduite de la manière suivante : On 


trace au moyen d’un émetteur étalonné la courbe 
de résonance du circuit de façon à pouvoir modifier 


les éléments pour se rapprocher de la fréquence 


choisie. En même temps, on détermine la qualité 


du circuit. Il reste à mesurer l’amplitude de la 
tension entre chacun des dés et la masse du cyclo=« 


tron. Il est commode d’utiliser pour cela un kéno= 
tron dont la plaque est réunie à l’électrode dont om 
mesure le potentiel et la cathode à la masse par 


l'intermédiaire d’un condensateur shunté par une 
résistance importante de l’ordre de 108 ohms en 


série avec un microampèremètre. La capacité doit 
être, et de beaucoup, supérieure à la capacité 
interne du kénotron. Elle se charge quand la diode 


No x 


est conductrice la différence de potentiel à ses 
. bornes devient égale à la tension de crête, car l’in- 


tensité du courant traversant la résistance ne 
modifie pas de façon sensible la charge de la capa- 
cité puisque avec les valeurs choisies la constante 
de temps du circuit est très grande par rapport à la 
période des oscillations. Si on veut connaître pen- 
dant le fonctionnement du cyclotron la tension 
entre les dés on utilisera une petite sonde circulaire 
isolée de la chambre voisine d’un dé. Elle est reliée à 
une armature d’un condensateur dont l’autre arma- 
ture est à la masse, on forme ainsi un diviseur capa- 
citif. On mesure la tension aux bornes de ce conden- 
sateur avec un voltmètre électronique dont les indi- 
cations sont reliées aux résultats obtenus par la 
méthode précédente. 

Le cyclotron qui a été construit fonctionne avec 
une amplitude de tension entre les dés qui peut 
atteindre 16 000 volts. 


Déflecteur. — Pour éviter les perturbations 
apportées par le champ magnétique dans le fonc- 
tionnement des montages électroniques détecteurs 
de particules on est amené à diriger le faisceau 
d'ions sur une cible située en dehors du champ ; en 
même temps la cible est en dehors du domaine 
d'action de ia haute fréquence. La déviation est 
obtenue en général au moyen d’une différence de 
potentiel continue appliquée entre la masse et une 
électrode déflectrice : lame de cuivre parallèle à la 
paroi verticale du dé et qui est reliée au pôle négatif 
de la tension. 

La force d’origine électrique agissant sur les ions 
est radiale, elle s’ajoute à la force centrifuge 
{ = Mo?}fr ; de cette égalité on déduit en prenant 
les différentielles logarithmiques des deux nombres 


dfff = — drfr 


en remarquant que Mv?reste constant. 
La variation Af de f provient du champ électro- 
statique h agissant sur l’ion de charge ze c’est hze 


_ on a donc très approximativement 


Are he 
ne MDE 7 
ou, en désignant par V la différence de potentiel qui 


pourrait communiquer aux ions pris initialement 
au repos, leur énergie finale 


— Arfr = hr]2V puisque F) Mv? — ZeY. 


Cette relation permet de prédéterminer h quand 
on s’est fixé le rapport Ar /r nécessaire pour amener 
les ions à passer en un point déterminé. 

Le calcul n’est qu’approximatif car le champ 
entre le dé et le déflecteur qui agit sur les ions n’est 
pas constant. Les ions du faisceau pénètrent dans 
ce champ à peu près au milieu du dé où se trouve la 
fenêtre de sortie, c’est-à-dire à des instants pour 
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lesquels le potentiel du dé est voisin de zéro et de 
toute façon va en croissant ; le champ augmente 
donc ce qui est favorable à la déflection ; il faut 
remarquer encore que les 1ons sont conduits vers la 
périphérie des pièces polaires où le champ magné- 
tique diminue rapidement ce qui est une cause 
d'augmentation supplémentaire de rayon de cour- 
bure des trajectoires. 

Après avoir ainsi calculé une valeur approchée du 
champ, on cherche à recevoir le faisceau sur la 
sonde, placée dans un cylindre fixé à la chambre et 
assez long pour que la cible soit en dehors des 
lignes de force du champ magnétique. En faisant 
varier le potentiel continu appliqué sur le déflec- 
teur, le courant sur la cible passe par un maximum 
assez net. On peut ensuite agir sur la position de 
l’électrode déflectrice commandée à l’extérieur par 
la déformation d’une membrane. Il est commode 
d’ailleurs d'observer les dimensions et la position 
du faisceau en le recevant sur une sonde recouverte 
de sulfure de zinc fluorescent qui est vue à travers 
une glace obturant une ouverture de la chambre : 

L’intensité du faisceau reçu par la cible est 
réduite de façon très importante d'environ 90 % 
après déflection. Les ions n’ayant pas tous été 
accélérés au même instant à la sortie de la source 
n’ont pas la même énergie, le faisceau a une cer- 
taine largeur par l’action du champ électrostatique. 
Un certain nombre d’ions vient rencontrer les parois 
et est perdu. 

La tension continue alimentant le déflecteur est 
fournie par un montage redresseur classique ; elle 
peut varier à volonté de 2 000 à 30 000 volts. Une 
self de choc est intercalée entre l’électrode déflec- 
trice et le pôle négatif de la haute tension ; un 
condensateur est placé en dérivation entre cette 
électrode et la masse de façon à ce que les charges 
variables en haute fréquence du déflecteur puissent 
s’écouler au sol sans traverser le montage four- 
nissant la tension. L’électrode déflectrice est portée 
à l’extrémité d’une tige qui traverse un bloc de 
plexiglass servant d’isolant, la membrane et enfin 
une ouverture pratiquée dans la chambre du cyclo- 
tron. L’étanchéité est assurée à la manière habi- 
tuelle. 


Recherche des conditions de fonctionnement. — 
La mise en évidence du fonctionnement du cyelo- 
tron s’effectue au début au moyen d’une sonde : 
lame de cuivre ou de graphite portée à l’extrémité 
d’une tige rectifiée en acier inoxydable non magné- 
tique. La tige peut coulisser dans un tube de 
plexiglas qui traverse une plaque de laiton obturant 
uné ouverture circulaire pratiquée dans la chambre 
du cyclotron. La sonde est ainsi isolée électri- 
quement de la masse ; l’étanchéité au vide est 
assurée par joints toriques pour la plaque et le 
tube, un petit joint, dont le diamètre intérieur est 
égal à celui de la tige, est écrasé à la manière 
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habituelle ; la sonde est facilement déplaçable sous 
vide entre les dés. 

La sonde doit être placée à l’intérieur d’une boîte 
à parois de cuivre mince évidemment ouverte sur 
une face pour que les ions du faisceau puissent y 
pénétrer. Elle est reliée à la masse. La sonde est 
presque en contact avec la paroi verticale de la 
boîte la plus voisine du centre de la chambre. Ce 
blindage de la sonde est indispensable sinon le 
courant reçu par la cible est la somme du courant 
d’ions qui ont subi à partir de la source les accélé- 
rations normales successives et du courant pro- 
venant d'ions vagabonds beaucoup plus nombreux 
formés dans la chambre par le champ de haute fré- 
quence appliqué sur les dés. Ces courants parasites 
empêchent l’observation de la résonance quand on 
fait varier le champ magnétique. La courbe du 
courant reçu par la sonde en fonction du champ 
doit passer par un maximum aigu alors que sans 
blindage on obtient des courbes présentant plu- 
sieurs maxima extrêmement peu prononcés et 
aplatis. 

La largeur de la courbe de résonance permet 
d'apprécier le fonctionnement. Les courants sur la 
cible sont mesurés avec un galvanomètre de table à 
plusieurs sensibilités ; il est inséré entre la sonde 
et la masse et shunté par un condensateur servant 
de court-circuit pour les courants de haute fré- 
quence qui pourraient traverser l’appareil de 
mesure. 

Dès qu’on a obtenu une courbe de résonance il 
est facile de parfaire tous les réglages et de voir 
l'influence des anneaux de correction et de mise en 
forme du champ de façon à obtenir le maximum du 
courant reçu par la cible. 


Utilisation du faisceau. Les émulsions 
nucléaires permettent dans certains cas l’étude des 
réactions nucléaires dans la chambre même du 
cyclotron ainsi que la détermination de l’énergie 
finale des ions du faisceau. Elles sont commodes 
d'emploi par leur faible encombrement et leur 
sensibilité continue ; elles évitent d’avoir à ‘dévier 
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le faisceau par un champ électrostatique pour 


l’amener sur la cible dans une région où le champ 


magnétique a une influence négligeable sur les 
appareils électroniques de détection et où le champ 
de haute fréquence est sans action sur la cible. 
Pour un petit cyclotron on évite ainsi une dimi- 
nution considérable de l'intensité du faisceau qui 
est normalement petite car les ions du faisceau 
n'étant pas monocinétiques, leur déviation par un 
champ électrique au xiliaire conduit à un étalement 
important des trajectoires qui viennent rencontrer 
les parois. Le morceau de plaque nucléaire est placé 


derrière une petite ouverture pratiquée dans une , 


boîte métallique qui le contient et protège l’émul- 
sion de l’action des'ions vagabonds et de la lumière 
provenant de la source d’ions. Le plan de l’émulsion 
est incliné de façon à ce que les particules la 
pénètrent sous une incidence voisine de l’incidence 
rasante. L’examen des trajectoires et la mesure de 
la mesure de leur longueur sont ainsi facilités. 

Par cette méthode l’étude des désintégrations des 


éléments légers par les protons et les deutons est . 


relativement facile ; on peut déterminer le spectre 
d'énergie des particules émises au cours de la 
réaction nucléaire et leur distribution angulaire 
ainsi que la section efficace. 

Les cristaux scintillateurs concourent au même 
but, de plus certains d’entre eux sont sensibles aux 
photons y, ils permettent l’étude des captures 
radiatives, des diffusions inélastique, etc. 

Le cristal choisi suivant le rayonnement à détec- 
ter, est placé au voisinage de la cible il doit être 
recouvert d’un film mince opaque à la lumière 
toujours présente dans la boîte du cyclotron. Il est 
en contact optique avec un guide lumière constitué 
par un cylindre de quartz ou de plexiglas enfermé 
dans un tube d’un métal non magnétique. Les scin- 
tillations sont transmises par une suite de réflexions 
totales jusqu’à un photomultiplicateur ainsi éloigné 
du champ magnétique. Le tube comporte un épau- 
lement qui s'appuie Sur la boîte du cyclotron, 
l'étanchéité étant assurée toujours par un joint 
torique. Manuscrit reçu le 14 janvier 1958. 
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DOSAGE RADIOMÉTRIQUE RAPIDE DE L'URANIUM ET DU THORIUM 
DANS LES MINERAIS COMPLEXES 


Par:G JURAIN et J. P. MAILLOT, 


Centre de Recherches Radiogéologiques, Nancy. 


Résumé. — Exposé d’une méthode radiométrique rapide de dosage de l’uranium et du thorium 
dans des minerais complexes. 

Elle utilise fondamentalement deux mesures d’activité gamma à l’aide d’un ensemble à scin- 
tillation Nal(Tl) en se plaçant à deux valeurs différentes du gain du préamplificateur. On obtient 
ainsi deux résultats qui, associés aux résultats obtenus au cours d’étalonnage, permettent le calcul 
de la teneur en U et Th. Cette méthode de dosage s'applique uniquement à des minerais en équi- 
libre radioactif. : 


. Abstract. — À rapid radiometric method for the determination of uranium and thorium in 
complex ores. 

We use fundamentally two gamma-activity-measurements, with the aid of an assembly of 
Nal(Tl) scintillation detectors by fixing two different values of gain of the preamplificator. Thus 
we obtain two results which associated with those obtained by standardisation, allow the calcu- 
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lation of the percentage of U and Th. 


This method is only applicable for ores in radioactive equilibrium. 


La méthode décrite ici permet un dosage rapide 
de l’uranium et du thorium dans des minerais com- 
plexes à teneur moyenne, à l’aide d’un ensemble 

- de détection à cristal d’iodure de sodium activé au 
 thallium. Une pesée suivie de deux mesures d’acti- 
« vité y permettent par un calcul simple, et après 
… étalonnage du dispositif, d'obtenir les deux teneurs. 
“ Introduction. — Le dosage de l’uranium et du 
… thorium dans les minerais se fait classiquement 
par voie chimique (oxydo-réduction de Uf à US 
+ pour l'uranium ; précipitation de l’oxalate, grillage 
» et gravimétrie de ThO, pour le thorium [1]). Radio- 
« métriquement on peut étudier les rayonnements « 
… des familles U et Th à l’aide des plaques pour usages 
… nucléaires [2], ou par scintillation sur écran de 
… sulfure de zinc activé à l'argent [3]. Les rayon- 
1 nements y permettent également un tel dosage, 
… soit par méthode au compteur de Geiger-Müller [4], 
» soit par scintillation sur cristal de Nal activé au 
 thallium. 
- Il nous fallait une méthode simple, rapide, de 
- précision moyenne et permettant avec un personnel 
non entraîné, et un matériel simplifié et robuste, de 
- communiquer en très peu de temps les résultats des 
teneurs en U et Th, aux différents stades du trai- 
- tement physique d’enrichissement d’un minerai 
complexe uranifère et thorifère (uranothorianite) ; 
ces conditions éliminaient le dosage chimique, long 
et pénible, le dosage « par plaques nucléaires 
(longues manipulations pour le développement et le 
‘dépouillement des plaques impressionnées), et le 
dosage « par scintillation., délicat et perturbé par 
le radon. 


Méthode employée. — Il ne nous restait donc que 
le comptage y par scintillation. Mais nous avons 


d'emblée éliminé le dosage par spectrométrie [5] 


du fait de l’importance et de la complexité de 
l’appareillage. 


Principe. — Le nombre d’impulsions est pro- 
portionnel à la quantité de corps radioactif et la 
hauteur d’impulsion, dans ce comptage, est fonc- 
tion de l’énergie d'émission, énergie qui, dans les 
familles U et Th, s'étend de quelques keV à 
2,62 MeV. En traçant la courbe représentant en 
ordonnées la variation d’intensité du rayonnement 
et en abscisses une discrimination d’énergie 
(hauteur des impulsions admises dans le système 
numérateur), on observera une différence suivant 
la nature des émetteurs de la famille en équilibre 
avec ses produits de filiation. 

Si, pour une énergie donnée, 1 gramme d’ura- 
nium donne une ordonnée y, et 1 gramme de 
thorium une ordonnée y,, un mélange de X gramme 
d'U et de ! gramme de Th donnera une‘ courbe 
dont les ordonnées seront : 


y = Kys + Lys 
Or, sans tracer la courbe, si l’on se place à deux 
abscisses d’énergie différentes, m et 7, on aura 
pour les intensités 
EN M Yn = K.Yim + l'Yom 
ERA Un = KeUan À L-Uon 
donc deux équations permettant le calcul de k et L. 


Dans ce système, Yn et Yyn représentent respec- 
tivement le nombre d’impulsions enregistrées aux 
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deux abscisses d'énergie m et n. Par un étalonnage 
judicieux, on pourra appliquer la méthodes à des 
échantillons de teneurs inconnues (fig. 1). 


INTENSITE A 


MeV 


m n 


DISCRIMINATION D'ENERGIE 
HrG: 1; 


APPAREILLAGE. — L’appareillage utilisé est cons- 
titué : d’un détecteur à scintillation équipé d’un 
cristal de Nal activé au thallium fabriqué par 
Harshaw, de dimensions 25 mm X 25 mm, couplé 
optiquement avec un photomultiplicateur « Radio- 
technique » type AVP 50 (1), celui-ci étant relié à 
un préamplificateur DCS.I Gabon construit sous 
licence C. E. A. par L. I. E. (?) ; d’un bloc de comp- 
tage RPN3 « Mesco » () muni d’une haute tension 
et permettant, soit un comptage prédéterminé, soit 
une durée prédéterminée de ce comptage. (Ce bloc 
peut aussi admettre l’utilisation d’un passeur 
d'échantillons automatique, avec enregistrement 
des résultats) ; le bloc et le scintillateur sont com- 
plétés par un support de sources qui donne une 
géométrie bien définie, semblable pour toutes les 
expériences. Les sources sont constituées par des 
coupelles d'aluminium d’environ 29 mm de dia- 
mètre, contenant 2 grammes de l’échantillon réduit 
en poudre. 
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Le préamplificateur du détecteur a été modifié 
pour permettre l’entrée des impulsions dans le 
bloc RPN3. En effet, la hauteur maxima des impul- 
sions admises dans ce bloc doit être de 5 volts au 
lieu de 12 volts pour les autres blocs de comptage. 
Il a donc fallu supprimer un étage de préampli- 


(3) La Radiotechnique, 
Paris (11°). 

(2) L. I. E. Laboratoire Industriel d’Électricité, 41, rue 
Émile-Zola, Montreuil. 

(5) Mesco, 14, rue René-Coche, Vanves (Seine). 


130, avenue Ledru-Rallin, 


JOURNAL DE PHYSIQUE 


N° 4 


fication. Cette modification n’entraîne pas de 
variation de la valeur des gains de ce préampli- 
ficateur, qui demeurent respectivement 0,15, 1,5 
et 15. 

Tout l’appareillage est placé dans une pièce où le 
mouvement propre est relativement faible (le 
château de plomb ayant été éliminé pour per- 
mettre un transport aisé de l’ensemble). 


Mope oPÉRATOIRE. — Deux échantillons ayant 


été l’objet d’un grand nombre de dosages chi- 
miques servent d’étalons. Il s’agit d’uranotho- 
rianites contenant respectivement : 
a) | 20,0 % d'U 
b) 5,1 % d'U 


54,5 9, Th 
78,7 % Th. 


Deux grammes de poudre sont placés dans la 
nacelle et soumis au comptage dans la géométrie 
déjà définie. Les manipulations consistent en : 


19 Étalonnage. — A) Mesure du mouvement 
propre aux deux gains 0,15 et 1,5, les deux gains 
correspondent aux deux énergies m et n. Ces 
mesures sont exécutées pendant un temps sufli- 
samment long pour intégrer les fluctuations. 
B) Mesure du nombre d’impulsions aux deux gains 


pour l’étalon (a). Après déduction du mouvement w 


propre, on obtient la participation propre de cet 
étalon. C) Mesure du nombre TU aux 
deux gains pour l’étalon (b). 


Soient Ym1 €t Yny les comptages que l’on aimait | 


respectivement observés aux deux gains pour 1 % 
d'uranium. 
De même, soient ymo et Yo les comptages que l’on 
aurait relevés aux deux gains pour 1 % de thorium. 
Les participations propres observées dans les 
expériences B) et C) seront : 


B) Yma = Ka-Ym1 + la-Ymo (1) 
Yna = Ka-Yn1 + la-Yn2 (2) 

C) Ymb = Kb-Ym1 + lb-Yme (3) . 
Ynd = Kb-Yni + lb-Yne- (4) 


k et l représentant respectivement les pourcen- 
tages en Ü et Th de l'échantillon (a), 4; et l, ceux 
de l’échantillon (b). 

Les équations (1), (2), (3) et (4) permettent par 
conséquent de calculer les quatre constantes y»1, 
Ymoy Yn1s Une: 

20 Mesures. — On procède au comptage du 
nombre d’impulsions aux deux gains pour l’échan- 
tillon de teneurs inconnues U et Th. 

Après déduction du mouvement propre, on 
obtient deux valeurs 


(5) 
(6) 


Yme — U.Ym1 + Th.Yn2 
Yne = U.Yn1 + Th.Yne 


” e 
M bn AMAR ES 2 + 


& 


avec : 
PE Do -Yma — la-Ymvd 
PAU Ka — da Kb 
Ka: Ynib — Kp.Yma 
AS lb ka RÉ la ky 
Éus nm Una — lo: Ynb 
Yn1 el. 
ANNEE Ka-Ynb — Kp-Una. 


PET, kà 
Les valeurs de U et de Th peuvent s’écrire : 


__ Ynæ:Ym2 = Yma-Yn2 


cn (7) 
Yma-Yna Ym1-Yn2 


Th — Yme:Yns TT Yna:Yms 
Ym2-Yn1 = Ymi-Yne 


(8) 


30 Exemple de manipulation. — Nous avons dis- 


. posé ci-dessous un tableau montrant les opérations 


successives. (Les résultats sont ramenés au même 
temps de comptage 1 minute.) 


Gain 0,15 Gain 1,5 
Mouvement propre 3 51° 7 
Étalon (a) 1 046 (1 033) 49 833 (49 316) 
Étalon (b) 1467 (1154) 40 005 (39 488) 
Inconnu 884 (871) 41 327 (40 810) 


Les valeurs entre parenthèses représentent la 
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participation propre de la source (c’est-à-dire la 
valeur enregistrée dont on déduit le mouvement 
propre). 

Dans ce tableau, par application des formules (7) 
et (8), on obtient les pourcentages U et Th assez 
rapidement. 


Généralisation et limites de la précision du pro- 
cédé. — Ne disposant malheureusement pas d’un 
troisième étalon suffisamment bien dosé pour 
vérifier la validité générale du procédé, nous 
n’avons eu qu’une confirmation indirecte de cette 
validité. Nous avons modifié très légèrement la. 
valeur d’une résistance du préampli, faisant varier 
du même coup la valeurs des gains (de façon très 
faible). Les coefficients d'étalonnage s’en sont 
trouvés totalement changés. 


AVANT MODIFICATION 


Ym1 = 11,68 
Uma = 14,66 
Uni) 
Yna = 448 


APRÈS MODIFICATION 


Ym1 = 58,9 
Ym2 = 29,8 
Yn1 — 1 672 
Yna = 920 


Nous résumons dans le tableau suivant les 
valeurs obtenues dans les deux cas. 


TABLEAU 
ÉCHANTILLONS AVANT MODIFICATION APRÈS MODIFICATION CHIMIE 

U Th U Th Le] Th 
Uranothorianite (201 n° 7) 16,294 6%2 1:07 15 SL ANS HOT 16/0102 SNA 
Monazite n° 1 (BL0/R0: 8007 OS ES RHONE NCA 
Monazite RA 36 ODA TS EST 0929241300 DSL M GROS 
Monazite RA 43 DEN EMRE RO7 (SONO RUEETERCE — _ 
Monazite n° 2 (SR RES 107 02H07 MERS ROSE — -- 
Thorianite RA 57 SO TD DOS COCA A RANE 8,9:07%72 6007 
Thorianite n° 7 SANS PET TES EMA GORE — — 


La méthode voit sa précision limitée par diffé- 
rents facteurs dont les principaux sont l’erreur sur 
la pesée et la fluctuation dans le comptage. Ajou- 
tons encore le défaut d’homogénéité des poudres 


» échantillons. L’erreur sur la pesée est de 5 mg envi- 


ron, soit 1/400 sur les deux grammes que contient 
la coupelle. L'erreur sur le comptage sera beaucoup 
plus importante, car, si l’on veut une méthode 
rapide, on limite volontairement la durée des 
comptages. Or la fluctuation est distribuée suivant 
une loi de Poisson. Elle diminue lorsque le nombre 
d’impulsions comptées augmente. On a plus préci- 
sément une probabilité de 0,68 de faire une erreur 
inférieure ou égale à 1 4/n, où nest le comptage enre- 
gistré. L’erreur quadratique moyenne ainsi définie 
sera pour un comptage de 10 000 impulsions, 


1/100 en valeur relative. En pratique, il faudra 
faire des comptages de cet ordre. 

La limite inférieure de la méthode, dans les 
conditions où elle est appliquée, se situe aux envi- 
rons de 0,5 à 1 % de teneur globale des deux 
éléments U et Th. Cette limite peut se reculer par 
utilisation : d’un château de plomb qui diminue le 
mouvement propre, mais alourdit considérablement 
l’appareillage ; d’un cristal d’iodure de sodium à 
puits, qui permet d'utiliser une plus grande quan- 
tié de poudre échantillon (aux environs de 5 à 
10 grammes). 


Conclusion. — La méthode décrite permet de 
doser l’uranium et le thorium en effectuant deux 
mesures d'activité y à l’aide d’un détecteur à scin- 


38 À 
tillation, sur un échantillon finement broyé d’un 
poids de deux grammes. 

Cette méthode s'applique à tous les minerais 
n'ayant subi que des traitements mécaniques, à 
l'exception des traitements chimiques susceptibles 
de modifier l’équilibre entre les différents produits 
de filiation des deux familles. Elle peut être con- 
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duite par un personnel non spécialisé avec un appa- 


reillage couramment utilisé et très maniable. Sa : 


précision est satisfaisante, lorsqu'on a à étudier des 
phases d’enrichissement physique et que les teneurs 
ont atteint des valeurs notables. 


Manuscrit reçu le 28 janvier 1958. 
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ÉTUDE DU COMPORTEMENT THERMODYNAMIQUE DES RÉACTEURS NUCLÉAIRES (2) 


Par G. FOURNET, 


Société Alsäcienne de Constructions Mécaniques, Paris. 


Résumé. — Nous avons établi — au moyen d’une approximation mathématique — des expres- 
sions monômes des différents paramètres régissant le comportement thermodynamique d’un type 
très général de réacteur nucléaire. Nous pouvons ainsi prévoir immédiatement les conséquences 
d’une modification quelconque des données de base du réacteur. 


Abstract. — By means of a mathematical approximation we obtain monomial expressions of 
different parameters governing the thermodynamical behaviour of a very general type of nuclear 


reactor. 
reactor data. 


1. Introduction. — C’est grâce à la circulation 
d’un fluide de refroidissement que l’on peut en 
général extraire une certaine puissance d’un 
réacteur nucléaire ; l'étude du comportement ther- 
modynamique d’un réacteur est donc particu- 
lièrement importante. 

Le but de notre présent travail est d'établir une 
_ méthode de calcul de ce comportement et de déter- 
miner ainsi, en régime permanent, les valeurs des 
différents paramètres (débit, températures, pertes 
de charge) définissant l’état de ce fluide de refroi- 
dissement. 

La structure mathématique des équations défi- 
nissant les échanges de chaleur est telle qu’il 
paraît impossible d’obtenir des expressions 
maniables de ces différents paramètres ; la consé- 
quence d’une modification d’une donnée de base 
du réacteur ne peut être connue qu'après de longs 
calculs numériques. Nous nous sommes donc atta- 
chés à essayer de trouver au moyen d’approxi- 
mations mathématiques convenables une expression 
monôme des différents paramètres thermo- 
dynamiques et de pouvoir ainsi prévoir immédia- 
tement les conséquences d’une modification quel- 
conque des données. 


2. Modèle de réacteur. — Nous supposons que 
la partie active du réacteur de forme générale 
cylindrique (rayon À, longueur 2/4) présente un 
certain nombre de canaux régulièrement répartis 
- et tous parallèles à l’axe du cylindre Oz. Nous 
faisons lhypothèse que la répartition macro- 
scopique du flux ® de neutrons thermiques est 
régie par une loi du type 


Dr, z) = ®, J'o(ar) cos bz (1) 
où r désigne la distance d’un point quelconque à 


(1) Le présent travail a été développé à partir d’études 
entreprises par la S. A. C. M."pour le compte du C. E.'A. ; 
ces études avaient donné Henau rapports SACM LT 463 
(juin 1953) et LT 469 (juillet 1953) et plus particulièrement 
au rapport LT 509 d’avril 1954, 


Thus we can foresee immediately the consequences of any modification of fondamental 


l’axe Oz de la partie active et z la distance de ce 
point au plan médian de la pile perpendiculaire à 
l’axe Oz. ; 

Dans chacun des canaux ménagés dans le modé- 
rateur, se trouve disposée une substance fissible. 
Toutes les sections de tous les canaux sont iden- 
tiques ; la substance fissible au centre, est entourée 
d’une gaine qui l’isole du fluide ; cette gaine peut 
comporter des dispositifs augmentant les échanges 
de chaleur et nous désignons par / le « périmètre 
efficace » d’une section de la gaine. La section 
offerte au passage du fluide est S. 

Nous supposons que toute la production de 
chaleur a lieu dans la substance fissible et que la 
circulation du fluide doit être telle qu’en aucun cas, 
la température superficielle 4, des gaines ne puisse 
dépasser une valeur 7, (cf. chapitre 8). 

Nous supposerons également que les écoulements 
ont lieu en régime turbulent et que la valeur du 
coefficient d’échange À est fournie par la relation 


RD IX = 0,023 RS P+U8 2) 


où D est le diamètre hydraulique, À la conduc- 
tibilité thermique du fluide, & et & les nombres de 
Reynolds et Prandtl. 


3. Établissement des équations générales. — La 
répartition macroscopique du flux de neutrons 
thermiques étant donnée par (1), la puissance 
dissipée dans une tranche de longueur dz de sub- 
stance fissible est 


AJ (ar) cos bz dz (3) 


quand r et z sont les coordonnées moyennes de 
cette tranche. 

Si, suivant les notations habituelles [1], on désigne 
par z, le «rayon équivalent » de la cellule que Pon 
attribue à chaque canal, le nombre des canaux du 


réacteur sera 
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et la longueur totale des barreaux de substance 
fissible est ainsi 


R? 
72H. 
1 
La valeur moyenne À de la densité linéaire de 
puissance dissipée dans un barreau est donc 


quand P désigne la puissance neutronique totale 
de la pile et que l’on suppose que toute la pro- 
duction de chaleur a lieu dans la substance fissible. 

Si on désigne par « (variable r) et 8 (variable 2) 
les rapports entre la valeur maximum de cette 
densité de puissance et sa valeur moyenne, on 
obtient dans le sens radial 


4 
1 [ Jotar) 2x r d 
— olar) 2x r dr 


X —= 


rx R? aR 
AE R WOJAER 
TE f J(u) u du el 


et dans le sens longitudinal 


1  dH 
RNA sin bH 
2J cos bz dz 


Le paramètre À de l’expression (3) est alors 
défini par 


B 


Pr? Pa aR bH 


PE AE SR D ST an) ul 


Pour obtenir les équations de base régissant le 
comportement thermodynamique du réacteur, 
nous considérons les phénomènes thermiques dont 
un canal, situé à une distance r du canal central, 
est le siège ; les températures de la gaine et du 
fluide dans une tranche dz dont les coordonnées 
moyennes sont r et z sont respectivement notées 
ta (r, z) et O(r, 2), x 

La première équation est obtenue en écrivant 
l'égalité entre d’une part, la puissance transmise 
de la gaine au fluide qui circule dans le canal — soit 
pour une tranche de canal de longueur dz 


hllte(z, r) — O(z, r)] dz 


où À désigne le coefficient d’échange, / le péri- 
mètre efficace de la gaine — et d’autre part, la 
puissance dissipée dans la substance fissible 


AJ, (a, r) cos bz dz. 
La première équation est ainsi 
hlta(z, r) — O(z, r)] — AJ,(ar) cos bz. (4) 


La deuxième équation peut être obtenue en 
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écrivant que la puissance fournie au fluide sert à 
augmenter sa température d’une quantité définie 
par 


hllta!z, r) — 0(z, r)] ds — GSC d0 


(Cest la chaleur spécifique du fluide, G le débit 
massique unitaire) d’où 
dO(z, r) 4 


= Gt 50" AVole, r) cos bz. (5) 


4. Pile à fente. — 4.1. DÉTERMINATION DU 
DÉBIT MASSIQUE UNITAIRE DE CHAQUE CANAL. — 
Nous allons maintenant appliquer les équations 
générales que nous venons d’indiquer au cas d’une 
pile à fente ; l’entrée du gaz (à la température 6,) 
dans chaque demi-canal (compris entre la fente 


l 
Entrée du gez. par la fente 


Entrée du gaz per la fente 
F1G. 1. — Schéma d’une pile à fente. 


et une extrémité du réacteur (cf. fig. 1)) se fait 
dans le plan z —0 de sorte que l’équation (5) 
peut s’intégrer sous la forme 


O(z, r) — 6, = ser [cos bz dz 


Gr) SC 
SOIt 
GRR er _ AJ(ar) sin bz 
CES (6) 


La température de la gaine {4(z, r) est donc 
fournie par 
cos bz sin bz (7 
hi Gr) SCb (° ) 


Dans le modèle de réacteur que nous consi- 
dérons, nous savons que la température des gaines 
doit être inférieure à une température limite T, ; 
la valeur maximum de (7) devra donc être pour 
chaque canal (repéré par r) inférieure à Ty. Le 
calcul de 


D AT (arts ee 
dz 


1a(z, r) — 05 = AJ,(ar) 


sin bz cos bz 


hi Gr, 5Cb |. 
fournit pour valeur z, de z correspondant au 
maximum de (7) 
sin bzm _ hl 


cos by  Gtr) SCb? 
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nous obtenons donc pour la valeur maximum de (7) 


} 1 1 
ta(r, Z)max — 00 = AJ (ar) VAE PB G2(r) S2 C? b? (8) 


Le débit massique unitaire du canal se trouvant 
à la distance r de l’axe de la pile est par consé- 
quent défini par : 


D HA de, 1. 
AJ, {ar) Vs DE F [b G(r) SCT? i 


La relation fournissant la valeur du coefficient 
d'échange À étant mise sous la forme (cf. (2)) 


RAD 02 COR) 2216, 


soit tous calculs faits 
h — 0,023 Cu+0.2 B—238.D—0.2, 08, 


nous pouvons poser 
hl = s. GS (10) 
avec 


s — 0,023. Cu? —28,D—02 J, (11) 


De même, nous écrirons pour simplifier des 
expressions ultérieures 


bSCG=t.G; (12) 
l'équation (9) devient ainsi 
To — 0% V 1 1 
Aile)  V 5 Gi Tac a) 


La méthode de calcul développée dans l’annexe I 
permet, après avoir posé 


1 1 


SOI 
k== GI — ce ...), (14) 
d'obtenir (cf. (68)) 
it paft}rn 
AJlar) Re? nn) A 
k?+1 
soit 
a k? 1 1+0,8%? 
Ta— 6  Vk+1 ALLEZ y P+i c + 
PAT en) 0 É ; 
+1 
k 
d’où 
G{r) (15) 
1+%? 1 1+%? 
= (42 +1) STFO8 [pe sk 5] TOR pie] 10,88. 
CT 9 
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Tous calculs faits, nous pouvons écrire : 
Fe 0,5(1+%2 1 
Gr) = Ë 1 A PS | 
(0,023 k)F° 
DIRES 1+%? 2 1+7%2 
1+0,8X%? 1+0,8X4? ,8k2 
b À (PAS ON) a SE 
1+%? 0,2%? 5e 2k? 
6 1+0,8%° 1+0,8%? æ 3(1+ 0,8%? 
“ 0‘2%? k? 1 1+%2 
1+0,8%? 7 1+0,8%? 7 1+0,8k 1+0,8%2? 
LT) ET der .[J,{ar)] *, (6) 
4.2. CHOIX DU COEFFICIENT #. — Nous choi- 


sirons la valeur de k (cf. (14)) telle que pour un 
certain canal moyen repéré par r, (auquel corres- 
pond G;,) on ait exactement 


1 1 
Le 
s GY tGo 
soit 
et os 
k=-6@ (17) 
À partir de (13) on obtient alors 
To— 09 (Uk? 1) 
AJ (ar) 1@ Me) 


l'élimination de G, entre (17) et (18) nous fournit 


k AJ ,{(aro)\°:2 

Re (re o) 
soit, après avoir remplacé {et s par leur expression 
(cf. (11) et (12)), 


(19) 


k A 
[C2 11 7 0,02 


2 

3 b0,8. C—9,2 dan P+ à 
.D°,? FA 
Après avoir calculé la valeur du,second membre 


de (20) la courbe représentative des variations de 


2 k 
RUES 


permet de déterminer la valeur de k convenable. 


y(X) (21) 


4.3. PERTE DE CHARGE. — Nous ne calculerons 
que la perte de charge relative aux canaux de la 
pile. Pour obtenir la perte de charge totale, il 
convient donc d’ajouter ce qui concerne l’extérieur 
de la partie active du réacteur. L'expression 

2G? fx 
Ap = Die (22) 
permet de calculer la perte de charge corres- 
pondant à un canal de longueur x. Cette expression 
suppose que la valeur du rapport f/o n’évolue pas 
le long du canal ; quand cette valeur varie il faut 
prendre une valeur moyenne. 
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En remplaçant f par son expression 


—0,2 
f = 0,046 R—02 = 0,046 (23) 


on obtient 
Np=10,002 7 Gr 6m Le DRE 


L'expression précédente de @G nous permet 
d’obtenir 


(24) 


2 1)0,5@9+ 1) 1,8 LE 1,8 
Apr) =; 0,092 res | 1+0,8%? .b 1+ 0,8%? H 
1023 
1,8+1,8X%? 1,8+1,8k* 
1,0+ 0,8%? 7 1,0+0,8%? 
A Tao 0 
1+# 0,2—0,2%2 1,21 
LA T8TroSR ! T'arose 2 140,86 
pere 2: 
1,2 + 0,6%? 1,8% 1,8 À ;gare 
7 1,0+0,8%? ,  1,0+0,8%? 7 1,0+0,8%? ”"1+0,8%? 
D TOFO8R jLOHOB +08 Per d 
(25) 
4.4. DÉBIT TOTAL. — Le débit volumétrique 


d’un demi-canal situé à une distance r, de l’axe du 
réacteur est égal à 


(rs) de G{r;) 5 
q À ‘ e . 

Pour obtenir le débit volumétrique total, il suffit 
maintenant d’effectuer une somme de termes ana- 
logues à (26) pour toutes les valeurs de r; réalisées 
dans le réacteur ; pour simplifier nous remplaçons 
cette somme par une intégrale où le débit (26) 
correspondant au canal r est multiplié par la proba- 
bilité. 


(26) 


5 DATE E 
TT 
de trouver un centre de demi-canal dans une cou- 
ronne circulaire définie par les rayons r et r + dr. 
Nous obtenons ainsi pour le débit volumétrique 
total 


Heure Gr) ST nds ES ii 
el à = af Gore 


T a? 
2 14 | 
= fr Lo] x de. (27) 


L’exposant de J,(x) peut s’écrire sous la forme, 


Le UNE VOS 
1 + 0,8 2 1 + 0,8 °° 


1 + 


cet exposant est en général, légèrement supérieur 
à l’unité. La valeur de l'intégrale 


aR ITR 
[ (Jo(x)] 1 F08E > dx 


peut s’écrire Z(u,(k)) en posant 


Ilu) =f" Ja) + x dr (28) 
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et 
O2 


1 + 0,8 k2 


Nous pouvons donc maintenant écrire l’expres- 
sion du débit volumétrique total Q sous la forme : 


ua (k) (29) 


à z 1 
ga [e [SE tt) 
1 1+% __1+# 
ET TF0  ;1+0.86 (Te —0}) 1-+0,8% 
+R. 0,2%° 2 
a ç 1+08r FE F. +34 +0.845) 
con 1 Es I TE - 
PRO ol PE SET 


Le débit massique s'obtient au moyen de 
M = Qe. 


4.5. TEMPÉRATURE DU FLUIDE A LA SORTIE DE 
CHAQUE CANAL. — La température de sortie d’un 
canal quelconque est définie par 


AJ (ar) sin bA 
OUT, r) — 05 = Oslr) — 00 ns 


soit, tous calculs effectués 


1 
Jose sin bA 


. " Fe TA 1)9,54 + #°) 
S Te = x € 
s k (0,023 k) 
0,84? 0,2k2 1+ k? 
7 1+0,84? 7 1+0,8k? 1+0,8%°? 
b A (Ta — 0,) È 
0,2k? 0,2k? Hd 0 
TT L+0,8% ; Ti+086 æ 34 +0,84) 
0,2K2? # x 0,84? _0,2k4 
1+0,84? 1+0,8%? 7 1+0,8k? 7 1+0,8k%? 
D L S .J{ar) HO (EN 


Pour obtenir la température moyenne 0 de sortie | 


des gaz on peut calculer la moyenne de 04(r) ou 


. bien se servir de la relation 


(6— 0,) MC = P (32) 
ce qui fournit 

ain LB AE 

0% MC OC eR 


d’où en nous servant des valeurs des rapports œ et f. 


(33) 


6— 6, 
Tao 
: 1 0,8k? 
É sa Se on-l0) _,—1r086 
(0,023 k) I(u:1(x)) 
0,284 0,2#? 
ed 1 £ ke3 
A +08 FE og +08 : 
0,2k? 0,2k° CL 
F1+08R  T1+0,8% ,  31+0.8#) 
(Es LL £ À 
0,2% LE 0,8k? 
140,84 IFO8K Ç 1+0.84 
D OitO8R IF Ç 1+088, (34) 


N° 4 


4.6. PuISSANCE DE POMPAGE. — Nous ne cal- 


- Culons ici que la puissance qu’il convient de fournir 


au fluide pour qu'il puisse traverser la partie 


_ active du réacteur (cf. 4.3). 


Si chaque demi-canal du réacteur est alimenté 
séparément, la puissance de pompage est égale à 


2 Ap(r).glr) 


tous les demi-canaux 
dont la valeur en première approximation est 
fournie par 


2x r dr 


R 
2 F Aplr).alr) (35) 
Tous calculs faits nous trouvons ainsi 
2 LS 1)0.5œ+07 28 
W, = 0,368 ie] 0 
Me (0,023 k)% : 
T0: 
Juh)at a Hb 
1+%? 1+%? 
2er A ar En PET 
* 1+0,8k? ” 1+0,8%? 
A (TG — 0) à 
1+%% 0,2—0,4%2 2,8k2 
ie ç *1+08# 1088 p150+0.80) 
__1:2+0,8k° 28 k? 1,8—0,8k4° 
ï 2 T0 0 872 TS 2 
D 1,2+0,4% 1 1+0,8%' 1+0,8X% (36) 
_ avec 
aR 
Ju) = p Jofx)280 +0 x de. (37) 


WT, — 0,368 [ 


D PERTE, 


Si tous les canaux sont alimentés à la fois (dis- 
positif classique) le débit massique g(x) corres- 
pondant à chaque demi-canal devra toujours être 
associé à la perte de charge maximum Ap(0) de 
sorte que dans ce cas, la puissance de pompage est 
donnée par 

W = Q.Ap(0) (38) 
soit 


2,8 
1+0,8k7 


(&2 + De —] 
(0,023 k)* 
2,8 


7 1+0,8k? 


‘ j I(u;(k)) a? tm lo H 


0,2—0,4%2 
1+0,8%? 


2,8k? 


LS ANER 7 1,5(1+0,8%2 
Ps 0 173 : + 


1,24 0,48? 
7 1+0,8k? 


2,8k? 
at 1+0,8%? 


1,8—0,842 
7 1+0,8k? 


S (39) 
.. La seule différence entre W, et W, réside dans 

- l'emploi soit de J(w,(k)) soit de Z(u,(k)), évidem- 
ment la première intégrale possède une valeur plus 
petite que la deuxième intégrale. 
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9. Réacteur sans fente. — 5.1. DÉTERMINATION 
DU DÉBIT MASSIQUE UNITAIRE. — Nous allons 
maintenant montrer rapidement comment une 
méthode de calcul similaire peut être appliquée au 
cas des piles sans fente (cf. fig. 2). 


Entree du gez 


Fic. 2, — Schéma d’une pile sans fente, 


L'entrée du gaz (à la température 0,) s'effectue 


dans le plan z = — Æ de sorte que l’équation (5) 
peut s'intégrer sous la forme 
__ AJ (ar) 
O(z, r) — 0, — Gr) SC CG bz dz 
SOIt 
O(z, r) — 0, = ue [sin 6: + sin bH]. (40) 


La température {,(z, r) de la gaine peut donc 
s’exprimer sous la forme : 


cos bz 
hl 


sin bz + sin bH 
Gr) SCb 


ta(z, r) — 05 = AJ (ar) 


} tr) 
et un calcul analogue à celui exécuté dans le 
paragraphe 4.1 fournit 

ta(3, T)max — 00 — AJo(a, r) 


sin bH HAS ENCRES 
Gi sce * V y Ÿ [On sc 


. (29) 


Le débit massique unitaire du canal se trouvant 
à la distance r de l’axe de la pile est par conséquent 
défini par : 
sin bH 1 


Car) $ + 


Pai==0 


AJ{ar) 


(RL)? 


1 

+Hoan ss l‘) 

L'exemple montre que sin bA est très voisin de 
Punité ; nous pourrons donc employer en première 
approximation l’expression (69) que nous avons 
établie dans l’annexe. En utilisant les mêmes nota- 
tions que celles définies au $ 4.1 nous obtenons 
donc : 


Tae0: 
AJ {ar) 

1 
A VE Ti 1 VE 1 \ve 
MAN ne) PR (EE 


&k À 


ou nous rappelons que # est défini par 
à il 
SG Can [ 
La fonction e(r) étant toujours petite devant 
l'unité. 
La relation (44) fournit pour valeur de G(r) 


4 + efr) + ...] (45) 


VIE+1 VR+1—1 
Greta VE En) NPA rose RER Le 
RENTE LL Vi 
y CVE FI +02. pret] 08/HF1+02. (46) 
(Fee) 


Tous calculs faits nous pouvons écrire 


Gr É + cr : 


0,8VR+1+0,2 b O0,8VÆ+1 +0,2 


A ER. 


Vri+1 Vri+1 
a EVE 
ADSVRF1 + 0,2 (TEE 0 8 0,8VK2+1 +0 


Via+1 0,2(VR2+1—1) 2. VRt+1—1 
C O8VH+I+0,2 OSVR+I+02 æ S08V#+1+02 
0,2(VH+1—1) VIA +1 —1 _ 5 


0,8V43+ 1 +0,2 


DOS8VE + 1 +0,2 l 0,8Vk?+1 +0,2 S 


Vr+1 


Jo{ar) SV +1 O2 0 ET) 


5.2. CHOIX DE LA VALEUR DU COEFFICIENT #. — 
Nous choisirons la valeur de k de telle façon, que 
pour un certain canal repéré par r, (auquel corres- 
pond G,) on ait exactement 


HONG DEEE 
Dern ii (48) 
À partir de (43) on obtient alors 
Te—00 _sinbH+VR +1 1+Vk+1 (49) 
AJ (ar) Go 1 


l'élimination de G, entre (48) et (49) nous fournit 
k LEO) RE, 
(1 + Vr au)? TG— 805 


soit encore, après avoir remplacé t et s par leur 
expression : 


= 8 1 ( (50) 


k 1 RARE En 
(Ag Dep on DE UNIES 
D+92 J—1 $08. Feel" À (51) 


Après avoir calculé la valeur du second membre 
de (51) pour un modèle de pile, la courbe repré- 
sentative des variations de 


k 


[1 + Vr Fil? 62) 


Yalk) = 
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permet de déterminer la valeur de k convenable 
pour la pile considérée. 


5.3. PERTE DE CHARGE. — Nous pouvons main- 
tenant calculer la perte de charge correspondant 
à la partie active du canal en appliquant l’expres- 
sion (22) avec x — 2H ; le résultat final est 


1,8 
Aplr) = 0,18 [ie =. EVE ETES 
(0,023 k)VB +1 —1 


RE MEN LL LSVR+1 
0,8/k2+1 +0,2 
b H ADSVEE +1 +02 (Te — 


—0,2Vk?+1 +0,4 
ur 08VRi+ 1+0,2 \æ 


D. 1,8Vk?+1 
0,) OSVE FI +02, 
VE +1—1 
7 0,8VETF1 +0,2, 


1,8VR2+1 
0,8Vk2+ 1 +0,2 


pue i0 


_06/#+1 + 0,6 me 
D 98VE+1+0,2 


VI +1—1 LETTRE P 
O8VR+1+02 ç SV T1H02, 
L8VI+1 
Jolar) 0,8Vk?+1 +02. (53) 


5.4. DÉBIT TOTAL. — Le débit volumétrique 
total peut s’obtenir au moyen de 


DE Gr) Sir Ares 
of il = af 4 G{r) r dr 


FU: 
2SG(0) far En 

ne en af [J'o(x)] 0808 +1 +02 x dr. 
L’exposant de J,(x) 


(54) 


peut s’écrire 


(VIE ETES) 


VE? + = —— 
à VIS TDR 


R2 + 1/0,8 VK? + 1 + 0,2 = 1 + 
cet exposant est en général légèrement supérieur à 
l'unité. La valeur de l'intégrale peut s’écrire 
I(u,(%)) où Z(u) est définie par (28) et où w,(4) 
est égal à 

VR? + 1—1 


k ES yen Ale ie 
si EVE +1+1 


(55) 
Nous pouvons maintenant obtenir l’expression 
générale du débit volumétrique 


: 1 

O = (a + VRAI om 02 Tla(#) 

(0,023 AVFFI—1 "3 
DAMES RES 
0,8V#?+ 1 +0,2 


a? Gb 
VE +1 JAP VE 
A MS8VE+1 +0,2 (TEE 0,8VR+1+0,2 
TS VeAI __02(VR +1 —1) re VEÆ+F1—1 
p—1C OSVEFI+02 LU O8VE+1+02 0,8VR?+1 +0,2, 
0,2(VE+1—1) VR+1—1 VE +1—1 


DOSVEFI+02 y OSVAFI+0,2 s "oser +02, (56) 


Le débit massique s’obtient au moyen de 
M = @e 


No | 


A di 


N°4 


nous ne donnerons pas l’expression de M très voi- 
sine de (56). 


5.5. TEMPÉRATURE DU FLUIDE A LA SORTIE DE 
CHAQUE CANAL. — La température de sortie d’un 
canal quelconque est définie par 


_ AJi(a r) 


nn eo st bH (57) 


soit tous calculs effectués 


ee, SITLO PT 
6, (r)—-0,— er | 0,8VX?+ 1 +0:2 


(0,023 4JVF +1—1 


__ 0.8V#?+1—0,8 O,2VR2+ 1 —0,2 AVE 
O,SVÆ FI +0,2 7 o8Vr ! BUT? 0,2 
b + A OSVE+1+02 (pe 6) CUP e : 


0,2VK2+ 1 —0,2 0,2(Vx2+ 1 —1) 2 VRI +1 —1 
COSVA +1 + 0,2 u 0,8Vk+1+0,2  æ 0,8VR2+ 1 +0,2 


VR+1—1 


__O2(VR+1—1) 
1 O8VE+1+0,2 


D OSVE+1+0,2 


_ 08(Vr?+1—1) 
S 08VE?+1+0,2. 


NUE 
Jo{ar) 0,8VR+1+0,2 (58) 


La température moyenne @ de sortie peut être 
définie par 
(6—6,) MC = P 


ce qui fournit 


. P A.97r R°H 
0 = 59 
OeCrt af 591 


D we 


En nous servant des valeurs des rapports « et B 
nous obtenons 


1 


Ne 65 LEVRTE A1 | osVerrr0 à 
| rer 10e Le | is 
- 5 VA+1 +1 ve VRei+1—1 
I(0) pb  O8VREI +02 ( A ] * O8/HEI +0,2 
I{u3(k)) Pa" 68 
VIF 1—1 VRE+1—1 VRi+1—1 


0,2 


C OSVA+I+0,2 


"OSVR+I+02 0,824 1 +0,2 


ms 


VE I—1 ÿ VRE+1—1 a VR+1—1 
D O.SVR+1+0,2 7 OSVE+1+02 0,8VET+1+0,2 (60) 


4.6. PUISSANCE DE POMPAGE. — Si chaque canal 
du réacteur est alimenté séparément, la puissance 
de pompage est en première approximation déter- 
minée par 


277 dr 


[avg 
0 


Tous calculs effectués nous trouvons ainsi 


2 
Tia 


VIA 2,8 
(1 + VE + nn e +0,2 
(0,023 k)VE+1—1 


0,368 | 


COMPORTEMENT THERMODYNAMIQUE DES RÉACTEURS 


45 A 
F5 RARE ER EM 1 2,8VI+1 
0.8VIP+1 + 0,2 —2 #4 (2 _|osvmri+02. 
J(uA(k)) b a? a? H ( 


28 VR+1—1 
œ150,8VR+1 +02, 


+ 0,6—0,4V%2+ 1 
u 0,8V#2+1 +0,2 


2,8VRk?+1 
0,8Vk?+ 1 + 0,2 


Cane 


0,8Vk?+ 1 —2,6 
0,8V#2+ 1 +0,2 


VR?+1—1 
OSVEEF1+0,2 


—0,4(Vx?+1 +2) 
D 8V#+1+0,2 l 


2,8 


2,8VE?+1 
Jolar) SV +1 +0,2. (61) 
Si tous les canaux sont alimentés à la fois, la 
puissance de pompage sera donnée par 


W = Q.Aprax 


le calcul montre que l’on obtient une expression 
analogue à celles que nous venons d’indiquer à cela 
près qu’il faut remplacer l'intégrale J(u,(k)) par 
I(u,(k)). 


(62) 


6. Utilisation des différentes expressions. — La 
première opération consiste à déterminer la valeur 
de Æ suivant, soit l'expression (20), soit l’expres- 
sion (51). Les membres de droite de ces expressions 
sont en 

Jo (arc) ; 


l’exposant 0,2 montre qu'une variation relative 
importante de r, n’entraîne qu’une faible variation 
relative de 


k ou fi ke 
(k? + 1)05 


GEVE TA)? 

et de Æ par conséquent. Le choix de r, n’est donc 
pas critique et l’on peut prendre RÀ/2 ou encore la 
valeur telle que J?2(a,x) soit la moyenne arithmé- 
tique de 1 et Jt?(aR) ou tout autre valeur com- 
mode, c’est-à-dire correspondante à une valeur de k 
facile à manier. Signalons à ce sujet que pour un 
certain nombre de réacteurs on peut prendre # — 1; 
notre, premier travail dans ce domaine était jus- 
tement limité à ce cas (Rapports SACM LT 469 
et LT 509). 

Avant de développer la méthode générale de 
calcul nous voulons indiquer que dans toutes les 
expressions que nous avons données nous avons 
toujours groupé séparément en disposant des 
points les facteurs suivants : 

les facteurs faisant intervenir la géométrie du 
réacteur : forme et dimensions de la partie ac- 
tive (À, ÆH, a, b), distance entre canaux (par 
l'intermédiaire de x) ; 

les facteurs faisant intervenir la « macrothermo- 
dynamique » du réacteur : puissance neutronique 
(par l'intermédiaire de A), différence T4; — 0, de 
température entre le point le plus chaud des gaines 
et le fluide à l’entrée ; 

les facteurs définissant le fluide de refroidisse- 
ment : masse volumique p, chaleur spécifique à 
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pression constante C, viscocosité u, nombre de 
Prandtl € (nous avons toujours éliminé la conduc- 
tbilité À, en faisant intervenir la définition du 
nombre de Prandtl) ; 

les facteurs définissant la géométrie d’une section 
de canal : diamètre hydraulique D, périmètre 
efficace / et aire de la section de passage S. 

Dans le calcul des membres de droite des équa- 
tions (20) et (51) définissant la valeur de k les deux 
premiers groupes de facteurs ne donnant pas de 
difficulté. 

Les facteurs du troisième groupe ont des valeurs 
variables avec la température. Le premier stade du 
travail consiste à effectuer une estimation des 
conditions moyennes (température et pression) du 
fluide au sein du réacteur ; on obtient ensuite à 
partir de (20) et (51) la valeur de k correspondant 
à ces conditions moyennes. On peut ensuite cal- 
culer, Æ étant connu, les valeurs des différents para- 
mètres et en particulier les répartitions de la tempé- 
rature et de la pression au sein du fluide. Si les 
valeurs moyennes ainsi obtenues sont par trop 
différentes de celles considérées a priori on recom- 
mence le travail en tenant compte de ces nouvelles 
valeurs. Ces approximations successives doivent 
être également effectuées quand la méthode de 
calcul que nous préconisons n’est pas utilisée et ne 
constitue donc pas une complication liée à cette 
méthode. 

Le seul facteur du quatrième groupe qui peut 
présenter une difficulté est le périmètre efficace l ; 
ce périmètre efficace dépend en général du coel- 
ficient d'échange k et du débit massique unitaire G 
par conséquent. Une première solution, calquée 
sur celle que nous venons d'indiquer dans le para- 
graphe précédent, consiste à déterminer la valeur 
de k à partir d’une évaluation /, de  ; la valeur de 4 
ainsi obtenue permet de calculer les différents 
paramètres et la fonction G(r) en particulier ; à 
partir d’une valeur moyenne Gbien choisie on peut 
ensuite évaluer /. La comparaison entre cette nou- 
velle valeur Æ et /, montre si l’on doit ou non pro- 
céder à un nouveau calcul. 

Une deuxième solution, qui constitue un nouvel 
avantage de notre méthode, peut être employée. 
On commence par mettre la relation 


1= 1U(G) 
sous forme d’une expression monôme 
11 Gi 
ce qui permet de poser 
hl = s' G8+5; 


un calcul similaire à celui que nous avons donné 
dans le cas où l’exposant de G est 0,8 peut se 
poursuivre. Nous avons effectivement employé 
cette méthode pour résoudre un problème qui nous 
était soumis. 
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7. Résultats. — Cette méthode de calcul, bien 
que basée sur une approximation, nous a fourni 
des valeurs en très bon accord avec des évaluations 
effectuées au moyen de longs calculs. Le principal 
avantage de cette méthode est de permettre de 
voir à chaque instant ce que l’on fait et de «sentir » 
ainsi l’ensemble du comportement thermo- 
dynamique du réacteur. 

Nos expressions peuvent être également d’un 
grand secours lors de la conception d’un modèle de 
réacteur. Si l’on cherche par exemple, toutes 
choses restant égales d’ailleurs, quel est le fluide 
qui conduit aux plus hautes températures de sortie 
dans une pile à fente il suffit de comparer les 
valeurs de (cf. (31) 


0,2k? 
1+ 0,842? 


0,2? et 
C PE Mo æ 3(1+ 0,8%?) 
pour les fluides envisagés, la valeur de k est déter- 
minée pour chaque fluide à partir de l’expres- 
sion (20) ; on peut également ne considérer qu'une 
valeur de k correspondant à un fluide «moyen ». 
Nous avons souvent remarqué dans des calculs 
analogues que cette dernière méthode ne conduit 
qu’à des erreurs faibles et qu’elle peut être effecti- 
vement utilisée. 


8. Conclusions, — Notre méthode de calcul peut 
être facilement adaptée au cas où le coefficient 
d’échange est fourni par une relation autre que la 
relation (2). Il est également possible de prendre en 
compte d’autres hypothèses ; signalons seulement 
à ce sujet que nous avons mis au point un calcul 
analogue pour étudier le comportement thermo- 
dynamique d’un réacteur où l’on désirait que la 
température de la substance fissible ne dépasse en 
aucun cas une valeur fixée Ty. 

L’idée de base de notre travail est la transfor- 
mation d’une expression quelconque (cf. l’annexe 
et l’expression 64 par exemple) en une expression 
monôme approchée. Nous pensons que de tels pro- 
cédés de «factorisation » sont susceptibles d’être 
employés de façon générale pour résoudre certains 
problèmes technologiques ; l'avantage de ces pro- 
cédés est de fournir des expressions monômes des 
différents paramètres et de permettre par 1 même 
l’étude simple et rapide de l’influence des données 
sur ces paremètres. 

Nous sommes heureux de remercier MM. R. Julia 
et P. Herreng, Directeur Général et Directeur de 
la S. A. C. M. pour nous avoir autorisé à publier ce 
travail. 

Manuscrit reçu le 24 janvier 1958. 


Annexe 


Nous allons chercher la meilleure expression 
monôême fournissant une valeur approchée de 


vy + Va? + y (63) 


No 4 
quand le rapport de + à y est voisin d’un nombre 
quelconque # ; nous posons 

æ—= ky(1 +) (64) 


où € est petit devant l’unité. En ne conservant 
que les termes du premier ordre nous obtenons : 


Va y VR21 E 2e E ...) LA 
y ——— e k2? 
=yve rtf haie a k 


Nous cherchons maintenant à déterminer E, «, 8 


COMPORTEMENT THERMODYNAMIQUE DES RÉACTEURS 


47 À 


Y+VR FA 
E = 1. (67) 


kRE+1+YVR+1 


Dans le cas où y est nul nous pouvons donc 
écrire 


ARC Et x 1 
» 2 LL 
Me ie Fr Ses (68) 
RH 


de même dans le cas où y est égal à l’unité 


es HE VEr E y? 
ee rai Pa 
Ext y8 — Vk2 LA lee RON z k 1HVR TI 
p=yr+VE + Doit AVR nr rives à (80) 
par identification nous obtenons RE +1+ VE 
2 . 
1e k (65) soit encore 
RÉ NM PET V 1 1 
A ALTe PE MIE ARE © k2 1 lon es = 
8 1+yvRe +1 2 De PNEU Ne Ver 
nn FRA 
RARES kR VEHI 
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ÉTALONNAGE D'UN CHRONOTRON 
. A L'AIDE D'UN GÉNÉRATEUR D’IMPULSIONS DU TYPE ROSENBLUM 
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Résumé. — En vue de réaliser un spectromètre à neutrons rapides par temps de vol, un chro- 
notron a été construit en utilisant une distribution de lampes du type « 6 BN6 » à laquelle a été 
ajouté un système d’amplificateurs différentiels, permettant d'augmenter leur temps de résolution. 
Pour étalonner cet ensemble, un générateur d’impulsions de temps de montée 2.10—° s, de lar- 
geur 1085 et d'amplitude variant entre 1 et 45 volts, a été construit à partir d’un compteur Rosen- 
blum. Le temps de résolution du chronotron déterminé à l’aide de ce générateur d’impulsions a été 
trouvé égal à 5.10—10s, 


Abstract. — À fast neutron spectrometer using the time of flight method is described. The 
coincidence system is a distribution device of ‘‘ 6 BN 6 tubes. Resolving time has been impro- 
ved by differential amplifiers. To test this device, a pulse generator ofrise time 2-X 10 °sec. and 
with an amplitude from 1 to45 V has been constructed from a Rosenblum Counter. The chronotron 
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resolving time determined with this pulse generator has been found equal to 5 X 1029 sec. 


s 


Un spectromètre à neutrons rapides par temps 
de vol nécessite un système de coïncidence dont le 
temps de résolution doit être d’autant plus petit 
que le parcours des neutrons est plus faible et que 
leur énergie est plus grande. Lorsque pour des 
raisons d’angle solide, le parcours ne doit pas 
excéder 1 à 2 mètres, le temps de résolution devra 
être au moins égal à 10° $ si l’on veut obtenir des 
résolutions en énergie de 10 % env. entre 10 et 
20 MeV. 


1. Le chronotron. — En vue de réaliser un tel 
spectromètre, des mesures de temps de résolution 
ont été faites sur un chronotron dont voici le 
principe : vingt lampes détectrices de coïncidences 
du type « 6 BN6 » sont montées suivant la tech- 
nique des amplificateurs distribués. Les capacités 
parasites des deux grilles de commande de chaque 
lampe sont utilisées pour construire deux lignes à 
retard à constantes localisées. Les impulsions 
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issues des détecteurs sont envoyées aux deux extré- 
mités opposées de ces deux lignes et se propagent 
en sens inverse, le.temps de parcours d’une lampe 
à la suivante étant de 6.10-10 s. La lampe qui 
détecte la coïncidence fournit une indication sur le 
temps qui séparait les deux impulsions. On réalise 
ainsi un système dont les différents canaux sont 
distribués sur l'échelle des temps de 12.105 $ 
en 12.10-10 $. O’Neil (1955)[1] a déjà réalisé un 
système semblable mais avec un temps de réso- 
lution de 4.10—°s. 

La difficulté, lorsqu'on veut diminuer ce temps 
de résolution, vient du fait que l'impulsion de 
coïncidence apparaît simultanément sur plusieurs 
canaux. La courbe de l’amplitude de ces impulsions 
en fonction de leur position sur l’échelle des temps 
présente un maximum aplati à l'endroit précis de 
la coïncidence (fig. 2a). Un ensemble de 2 ampli- 
ficateurs différentiels et d’un redresseur (fig. 1) 
permet de transformer la courbe enveloppe précé- 
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F1G. 1. — Schéma du chronotron. 
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des différents canaux. 
a) A la sorite des lampes « 6 BN6 ». 
b) A la sortie du système différentiel. 


dente en une courbe plus aiguë (fig. 2b). Ce mon- 
tage est analogue au système différentiel de 
Bay (1956)[2]. Le temps de résolution de l'appareil a 
été étudié en analysant les différents canaux succes- 
sivement. Pour faire cette étude, il est nécessaire de 
disposer d’un générateur d’impulsions très brèves 
dont la largeur n’excède pas 105$ et dont l’ampli- 
tude varie de quelques volts à quelques dizaines de 
volts. 


2. Le générateur d’impulsions Rosenblum. — 
Aussi avons-nous pensé utiliser un compteur Rosen- 
blum (1955) [3], Sahe, 1957 [4] comme source d’im- 
pulsions de référence et avons réalisé un générateur 
dont les caractéristiques sont les suivantes : 

1. Temps de montée : 2,5 105. 

2. Largeur à mi-hauteur : 7 105. 

3. Impédance de charge : peut être réduite à 2 Q. 

4. Amplitude de l'impulsion comprise entre 15 
et 45 V pour une impédance de charge supérieure 
ou égale à 75 Q ; et comprise entre 1 et 10 V pour 
une impédance de 2 Q. 

9. Fréquence de répétition déterminée par l’in- 
tensité de la source de polonium. Dans le montage 
réalisé, la source de 20 u. e. s. (= 15 uC) donne 
environ 400 impulsions par seconde. 

6. Polarité positive. 

Le compteur Rosenblum utilisé est un compteur 
classique comprenant 8 fils de tungstène, tendus 
à 1,2 mm d’une plaque horizontale en acier inoxy- 
dable. En l’absence de rayonnement «, la décharge 
« couronne » obtenue est telle que le courant d’ioni- 
sation collecté atteint plusieurs centaines de micro- 
ampères. L’impulsion, ayant un temps de montée 
de l’ordre de quelques millimicrosecondes, est trop 
rapide pour que son amplitude dépende sensi- 
blement de la valeur de la résistance À,. On réalise 
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alors un diviseur de capacité (fig. 3). Si l’on sup- 
pose € > C5, l’impédance de sortie est déterminée 
par la capacité C, et la résistance À, en parallèle, 
cette dernière étant nécessaire pour la répartition 
des potentiels. Le condensateur C, variable, dont la 
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Sortie 
déclenchement 
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Fic. 3. — Schéma de principe 
du générateur d’impulsions Rosenblum. 


capacité résiduelle doit être très faible et égale à 
une fraction de picofarad, est réalisé par deux 
plaques parallèles, à écartement variable, connec- 
tées aux bornes « plaque » et «sortie » par des fils 
très courts. 


ÉTUDE DE L'IMPULSION. — Le signal est appliqué 
par l'intermédiaire de 20 m de câble coaxial (impé- 
dance 75 Q) sur les plaques de déviation verticales 
d’un oscallographe cathodique « Tektronix » 
type 517, pour compenser le retard (environ 10-75) 
introduit par le circuit de balayage déclenché. La 
bande passante de l’ensemble constitué par la ligne 
de 75 Q et le tube à rayons cathodiques, est telle 


Fic. 4. — Impulsions du compteur Rosenblum. 
Balayage : 140$ s/carreau. 

Sensibilité : 145 V/carreau. 

Impédance de charge : 75 (2. 

C; = 40 pF. 

Co = 1,9 pF. 
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qu’elle ne permet pas de déterminer avec précision 
le temps de montée réel de l’impulsion. 

On obtient une mesure plus précise du temps de 
montée en mesurant l’amplitude à la sortie du 
compteur Rosenblum, lorsqu'une ligne coa- 
xiale (75 Q) court-circuitée, de longueur variable, 
est mise en parallèle. On trouve par cette méthode 
un temps de montée égal à (2,5 + 0,2) 10 $ 
(fig. 4). 


3. Temps de résolution du chronotron. — Le 
chronotron décrit précédemment a été étalonné 
avec le générateur d’impulsions Rosenblum et la 
largeur à mi-hauteur de la courbe de résolution 
d’un canal a été trouvée égale en moyenne 
à 5.107 8 (f216a) 

Un essai a ensuite été fait en détectant à l’aide 
de deux photomultiplicateurs RCA 6810 (cristal 
plastique polystyrène “+ Popop) les rayon- 
nements 8 et y émis (1) par lSAu. 

La largeur à mi-hauteur de la courbe de réso- 
lution a été trouvée égale à 2 10° s (fig. 5b). 
L'écart avec le résultat précédent est dû d’une 
part au fait que les impulsions des détecteurs ont 
une largeur de 2 105$ au lieu de 107$ s avec le 
générateur Rosenblum, et d’autre part aux fluc- 
tuations statistiques dans les cristaux et les photo- 
multiplicateurs. 

Remarquons que ces résultats ont été obtenus 
sans amplificateurs distribués entre les phototubes 
et le chronotron, l’amplitude des impulsions étant 
dans nos expériences supérieure au seuil de 5 volts 
du chronotron. Avec des photomultiplicateurs 
radiotechnique 51 AVP et des rayonnements moins 
énergiques (surtout en formant les impulsions 
à 107$ s), il est indispensable d'utiliser un. ampli- 
ficateur distribué (CRC AMP 95) de gain 8 entre 
détecteur et chronotron. Dans ce cas, on a réussi à 
diminuer le temps de résolution à 1,5 107°s. 

Ce travail a été effectué au Laboratoire de Phy- 
sique Nucléaire de l’Université de Grenoble sous la 
direction du PT Bouchez qui nous a proposé ce 


(1) La période de transition de 1$Au est © 4.10—115 
qu'on néglige par rapport au temps de résolution du chro- 
notron. 
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FiG. 5. — Courbes de résolution. 
a) En utilisant le générateur d’impulsions Rosenblum. 
b) En utilisant deux photomultiplicateurs : RCA 6810. 


sujet et que nous remercions pour ses encoura- 
gements et ses nombreuses suggestions. J. Duclos 
est reconnaissant au C. N. R. S. de lui permettre 
d'effectuer ces recherches et remercie le PT Benoît 
de l’avoir guidé dans le domaine de l'électronique. 
Enfin, J. Fleury est reconnaissant au C. E. À. de 
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SÉLECTEUR MÉCANIQUE POUR NEUTRONS LENTS 


Par G. GOBERT et B. JACROT, 


Section des Réactions Nucléaires à Basse Énergie, C. E. N., Saclay. 


Résumé. — On décrit un sélecteur mécanique de neutrons lents. Ce sélecteur est constitué de 
deux rotors distants de 1 mètre, et faisant 4 mètre de diamètre, tournant à 6 000 tours/minute. 
La largeur de la raie obtenue est, dans les meiïlleures conditions, de 0,2 À et la durée de l’impulsion 
des neutrons de 50 us. Les problèmes mécaniques rencontrés lors de la construction sont analysés. 
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Abstract. — À mechanical selector of slow neutrons is described. 


(diameter 4 m) separated by 1 meter distance. 


the construction are analyzed. 


Depuis quelques années, plusieurs groupes de 
physiciens [1] se sont attachés à l’étude de la dif- 
fusion inélastique des neutrons lents par les solides. 

Le principe des expériences est d'envoyer, sur 
un échantillon de substance à étudier, des neutrons 
lents monocinétiques et de mesurer l’énergie des 
neutrons diffusés à un angle donné. Pratiquement, 
les neutrons incidents monocinétiques sont rem- 
placés dans beaucoup d’expériences [2] par la 
partie d’un spectre maxwellien qui a traversé un 
filtre cristallin [3], de béryllium par exemple. Avec 
cette technique, on a un faisceau composé de 
neutrons de plus grande longueur d’onde qu’une 
certaine longueur d'onde caractéristique du filtre 
(3,95 À pour le béryllium). Comme un spectre 
maxwellien décroît très vite avec la longueur 
d’onde, cette technique peut être utilisée. Cepen- 
dant, pour avoir des expériences plus précises, il 
faut avoir un spectre beaucoup mieux défini en 
énergie. Ceci peut être fait à l’aide d’un spectro- 
mètre à cristal ; c’est la technique qu’utilise Brock- 
house [le]. On peut aussi utiliser un sélecteur 
mécanique [1c]. Le principe général est le suivant : 
un premier dispositif laisse passer un paquet de 
neutrons assez court en temps. Puis, à une certaine 
distance, un deuxième dispositif ne laisse passer 
que les neutrons de ce paquet ayant mis un temps 
donné pour parvenir au deuxième dispositif. On 
réalise ainsi un monochromateur travaillant par 

impulsions, qu’on peut donc associer à un système 
d’analyse par temps-de-vol. Nous avons construit 
un tel appareil en nous imposant les conditions 
suivantes : 1) largeur à mi-hauteur de la raie pro- 
duite : 0,2 dans les meilleures conditions ; 
2) grande facilité de changer la longueur d’onde 
transmise sans changer la largeur de la raie ; 
3) durée de l’impulsion 50 us à mi-hauteur dans les 
meilleures conditions. 

Pour ceci, nous avons adopté un dispositif cons- 
titué par deux roues absorbantes tournant autour 
d’un axe horizontal et portant des fentes radiales. 


Itis constructed of two rotors 


The rotation speed is 6 000 r. p. m. The width 
of the line obtained is — in the best conditions —50 us. 


Mechanical problems encountered during 


Les roues ont 1 m de diamètre et sont situées à 1 m 
de distance. Une boîte de vitesse de rapport 1/1 
connecte les 2 roues et permet la variation de la 
longueur d’onde transmise. 


Déphaseur 
* mécanique 


, Moteur 
électrique 


j L \ 
Couronne de cadmium ne 
Fic. 1. — Schéma de principe du sélecteur mécanique. 


La figure 1 donne le schéma de principe de cet 
appareil. 

Nous décrivons 
parties du sélecteur. 


successivement les diverses 


I. Rotors. — a) DEscriprion. — Chaque disque 
est en duralumin monté sur un noyau en acier 
constituant l’arbre. Le poids de chaque disque est 
de 135 kg. Entre les diamètres 720 et 940 mm, se 
trouve le logement d’une couronne de cadmium 
d'épaisseur 1,1 mm, rendue prisonnière par une cou- 
ronne de même métal que le disque. Cette couronne 
de cadmium est percée de 8 fentes radiales 
(50 ‘X 16 mm) dont le centre est au rayon de 
420 mm. D'autre part, au rayon 250, et dans le 
plan radial défini par la première arête de chaque 
fente à neutrons, est percé un petit trou 
(3 X 0,16 mm) destiné à laisser passer un signal 
lumineux synchronisé avec le début de l’impulsion 
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des neutrons. La position de ce trou par rapport à 
la fente à neutrons est définie à quelques microns 
près. 

Les signaux recueillis par un photomultiplicateur 
(EMI 6099) ont un temps de montée à 
6 000 tours/mn de l’ordre de 1 micro-seconde. 

Le tableau 1 donne les caractéristiques méca- 
niques des divers matériaux utilisés. 
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_— Coupe d’un rotor monté dans son enceinte étanche. 
On voit le montage des joints tournants (a) et des roulements à billes (b). 


Ces rotors tournent dans deux carters à vide 
pour réduire la puissance consommée. 


b) Résisrance. — Le disque ayant des vitesses 
périphériques appréciables, nous avons déterminé 
un profil tel que le métal ait une répartition de 
contraintes acceptable. 

Recherchant une forme facilement usinable, 
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F16. 3. — Schéma de montage des fentes optiques et des fentes à neutrons. 
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TABLEAU I 
Poins MopuLe LIMITE LIMITE COEFFICIENT 
SPÉCIFIQUE  kg/mm? ÉLASTIQUE DE RUPTURE DE Poisson 

kg/dm* kg/mm? kg/mm? 
Rotors AU; G 1 vieilli 254 7 200 29 42 0,34 
Moyeux acier XC 38 f 8 21 000 50 85 0,3 
Visserie acier 30 NC 8 21 000 70 105 0,3 


Fic. 4, — Photographie d’un rotor. 


nous nous sommes inspirés du profil hyperbolique 
étudié par Stodola [4] pour les disques de turbines. 
Ceci étant, nous avons calculé la répartition des 
contraintes par la méthode graphique de Donath 
qui consiste à décomposer le disque en tranches 
successives d'épaisseur identique, et à étudier le 
comportement du métal aux lisières de chacune 
d’elles. 

Les résultats ont été portés sur la figure 5 où 
sont représentées, pour la vitesse de 6 000 tours/mn 
et à différents rayons, les contraintes radiales o,, et 


tangentielles oc. 


À titre comparatif, nous indiquons par les 
courbes Gro €Ù Sto les valeurs des contraintes dans le 


cas où le disque serait homogène, ne comportant 
ni cadmium ni couronne de serrage, ainsi que les 
valeurs 6, pour le cas où le disque serait également 
homogène mais non profilé (disque plan). L’épais- 
seur (25 mm) étant faible par rapport au rayon, 
nous pouvons négliger les variations que peuvent 
présenter les tensions radiales et tangentielles sur 
l'épaisseur du disque. 

Sur la même figure, nous avons tracé les SR 


Gr, €t ©, des contraintes radiales et tangentielles 


C 


Contraintes 
dans 
Ja couronne 


ere dans le cas du 
disque à faces parallèles 
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F1G. 5. — Courbes théoriques de répartition des con- 
traintes à 6 000 t/minute : 
a) disque plan homogène, courbes ot, ; 


b) disque profilé homogène, courbes 6t,, 67, à. 
c) disque profilé avec effet d’une couronne de 
cadmium, courbes Oty CÙ Cry 


Ce dernier cas correspond aux rotors effectivement 
utilisés. 


propres à la couronne. Comme nous avons admis 
des valeurs relativement élevées dans notre déter- 
mination des contraintes initiales, nous avons 
ajouté l’effet de cette couronne à celui produit par 
le cadmium. Pour cette détermination des valeurs 
initiales, nous avons considéré deux cas : 

19 Valeur de la pression spécifique produite par 
le cadmium dans son logement, plus pression pro- 
duite par la couronne supposée déformée. La valeur 
de la contrainte initiale qui s’ensuit a été supposée 
appliquée à la périphérie pour la détermination du 
profil d’un disque sans cavité. 

20 Détermination des contraintes initiales au 
rayon 360 mm (rayon intérieur de la cavité) en 
tenant compte d’un effet de jante sur la section de 
maintien avec tous les effets produits par l’en- 
semble constitué par la jante, la couronne et le 
cadmium. Les valeurs de contraintes à ce rayon 
dans les deux cas ne diffèrent que de 5 %,. L’usinage 
de ce disque a été exécuté par la Société RATIER 
à Montrouge suivant les mêmes méthodes que 
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celles utilisées pour les hélices d'avions : on com- 
mence par un forgeage du disque avec son profil, 
suivi d’une ébauche sur machine. Avant d'aborder 
les autres phases de cet usinage, les disques ont 
subi un vieillissement artificiel par variation de 
température dans des ambiances différentes. Sans 
changer les caractéristiques du métal, cette stabi- 
lsation diminue considérablement les tensions 
internes et évite toutes les déformations dues aux 
usinages et aux variations de température de 
Patmosphère auxquelles peuvent être exposés par 
la suite les rotors. 

Cette méthode nous a permis, après montage 
sur la machine, de constater à la périphérie des 
voiles n’excédant pas + 2/100 de mm. 


c) MESURES PRATIQUES. — Il était essentiel de 
vérifier que les contraintes expérimentales ne 
dépassaient pas celles calculées car nous avons 
adopté un coefficient de sécurité de 2. Ces mesures 
de contraintes ont été effectuées au Centre d’Essai 
des Moteurs et Hélices de Saclay (C. E. M. H.) 
par MM. Le Meur et Darricaud. 

La mesure des contraintes statiques produites 
par les effets précédemment cités a été exécutée au 
moyen de jJauges extensométriques du type 
AOIP 600 Q (18 b 25) fixées avec de la colle 


F1G. 6. — Schéma de montage des jauges de contraintes. 
À chaque jauge de mesure (a) est associée une jauge de 
compensation (b). 


Durofix. Les éléments résistants ont été placés aux 
rayons 170 et 303 mm. L'orientation des jauges au 
collage permettait la mesure des contraintes tan- 
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gentielles et radiales ; elles ont été jJumelées de 
façon à pouvoir obtenir des recoupements pendant 
les essais. 

Pour chacune de ces jauges, le facteur de sensi- 
bilité est À — _ ee +2; 1’où _ et . sont 
les variations relatives de la résistance et de la 
longueur. Les difficultés d’une telle mesure sont les 
suivantes : 

1) Elle doit s'effectuer en régime dynamique, 
c’est-à-dire pendant la rotation de la pièce expéri- 
mentée. Ceci oblige à faire passer le circuit élec- 
trique de mesure par un collecteur, organe délicat 
qui doit être de très bonne qualité et où le frot- 
tement des balais et les variations des résistances 
de contact peuvent avoir des conséquences impor- 
tantes. 

2) L'influence de l'humidité et de la température 
sur la résistance de la jauge doit être compensée en 
plaçant sur le même rayon que les jauges de mesure 
des Jauges témoins montées sur des plaquettes 
intermédiaires et soumises à aucun effort. Ces der- 
nières se trouvent ainsi dans les mêmes conditions 
(vitesse périphérique, frottement sur l’air résiduel, 
température) que les jauges de travail. 

Les résultats pratiques ont été reportés sur la 
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F1G. 7. — Résultats expérimentaux des mesures de con- 
traintes (décembre 1956). On voit que les contraintes 
effectives sont toujours inférieures à celles calculées. 


figure 7 et traduisent un niveau de contraintes 
effectives inférieur de 20 % environ à celui obtenu 
per le calcul. 


d) ÉQUILIBRAGE. VITESSES CRITIQUES. — La 
vitesse critique la plus basse de chaque rotor a été 
déterminée à partir du calcul de la flèche statique. 
Cette flèche est de 3,75 microns. La vitesse critique 
correspondante est de 15 400 tours/mn, ce qui est 
bien au-dessus de la vitesse de rotation utile. 

Par ailleurs, chaque disque a été équilibré stati- 
quement, puis dynamiquement. L’équilibrage 
dynamique a été fait au C. E. M. H. sur une 
machine Biring à une vitesse de 400 tours/mn. La 
précision obtenue est de 5 g/cem. Les corrections 
ont été faites par meulage du flan du disque. 

Les fentes à neutrons peuvent être obturées par 
paires de fentes opposées. Ceci a imposé un con- 
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trôle très strict de l'identité en poids des pièces 
d’obturation. 


IT. Accouplement. Boîte de vitesse. — La 
liaison des rotors doit être assurée de façon rigide, 
afin d’assurer une constance de déphasage entre les 
deux rotors. Pour cela, nous avons solidarisé les 
roues par un système à engrenages de résistance 
mécanique élevée et de haute précision. Nous avons 
accouplé les roues à ce transmetteur à engrenage 
par des accouplements du type « Comelor-Fouge- 
rolles » constitués de plateaux à denture extérieure 
montés sur les arbres et de coquilles à denture inté- 
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les vecteurs de rotation continue (fig. 8). Ce mou- 
vement auxiliaire est oktenu par un moteur à 
courant continu, entraînant un réducteur. Le 
réducteur est une application particulière du train 
de Peckeur [2]. Ce moteur peut également être 
alimenté par un stabilisateur électronique de phase 
construit par M. Eisenberg, qui a pour effet de 
maintenir la phase constante même si la vitesse des 
rotors varie. Ceci est utile car cette vitesse fluctue 
de + 1 % avec le temps. 


III. Bâti. — 4) Le bâti de cette machine est 
constitué d’une plate-forme en tôle soudée croi- 
sillonnée et dont les longerons principaux sont en 
profilés I de 140 mm de hauteur. Cette plate-forme 
de 920 kg est montée par l’intermédiaire de vérins 
hydrauliques sur des supports roulants indé- 
pendants. 


Les enceintes des rotors sont faites par un assem- 


blage de fonds emboutis usinés de 2 cm d’épaisseur 
et de 1,20 m de diamètre, la partie convexe est 
extérieure. Ces enceintes sont liées à la plate-forme 
par des vérins à vis, ceci afin de permettre l’ali- 
gnement des axes. 


IV. Puissance d’entraînement. Système de vide. 
— Cette puissance doit vaincre : 1) les frottements 
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rieure assurant la transmission du mouvement. Ces 
accouplements ont été équilibrés avec une machine 
€ GISHOLT » du C. E. M. H. avec une précision 
de 1 g/em. A l’arrêt, l’ensemble accouplé présente 
un écart angulaire maximum de { degré dû aux 
jeux. En rotation, ces jeux n’ont presque aucun 
effet et le déphasage reste stable à environ { us 
près si la vitesse de rotation de moteur reste cons- 
tante dans le domaine d’utilisation de la machine. 

Le transmetteur à engrenage (construction Maag 
à Courbevoie, Seine) permet de faire varier à 
n'importe quelle vitesse l’angle de phase, en compo- 
sant un vecteur lent, à un planétaire intérieur, avec 


Axe de rotation de la machine 


CA 
78 d 


RTE 


Moteur 
de déphasage 


Rapport 1/500685 
Fic. 8. — Schéma de principe du système d’accouplement des deux rotors. 


purement mécaniques ; 2) les frottements du disque 
sur l’air résiduel ; 3) les forces d’inertie et assurer un 
démarrage progressif. 

Cette puissance doit être réglable et permettre 
d'obtenir plusieurs vitesses continues ; pour cela, . 


6000 


1000 


F1G. 9. — Puissance absorbée par les divers éléments de 
la machine en fonction de la vitesse de rotation 
(novembre 1956). ° 
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nous avons choisi un ensemble Ward-Leonard. Le 
moteur d'entraînement des rotors de 6 kW assure 
un démarrage en 1 heure avec une grande résis- 
tance aux échauffements. 

Pour chaque rotor, l’inertie est sensiblement 
1,25 kgm/s?. Un couple de 0,4 m.kg et une puis- 
sance de 3,5 ch sont donc nécessaires. Pour vérifier 
nos prévisions énergétiques et connaître aussi par- 
faitement que possible les puissances de frottement, 
nous ayons procédé à une série d’essais dont les 
résultats sont portés figure 9. Pour ce qui concerne 
les frottements sur l’air, nous avons porté figure 10 


1000 


100 


100 1000 10 000 


Fc. 10. — Puissance consommée par un rotor tournant 
à l’air libre, en fonction de la vitesse de rotation. 


les résultats de nos essais de juillet 1956, effectués à 
- l'air libre, de 500 à 3 500 tours/mn. Ces essais 
concernent un seul rotor. Afin d'éliminer ce frot- 
tement, nous avons monté chaque rotor dans une 
enceinte de 250 dm* de capacité où le vide dyna- 
. mique est obtenu par une pompe à palettes 
de 38 m° de débit aux pressions ordinaires. 
En fonction de la pression, nous avons tracé, 
. pour deux vitesses, les courbes de puissance sur la 
figure 11. Ces résultats montrent qu’un vide de 
quelques millimètres de mercure est suffisant. 
Dans ce cas, l'étanchéité est assurée par un sys- 
tème de joints tournants du type PACIFIC, com- 
posé d’un grain en acier stellité rectifié, tournant 
avec l’arbre et frottant sur deux grains fixes en 
Carbone (vicarb). L’épaisseur de la partie frot- 
. tante est de 4 mm. La fibre moyenne de cette 
. partie frottante a une vitesse de 15 m/s. 


SÉLECTEUR MÉCANIQUE POUR NEUTRONS LENTS 


57 A 


La lubrification de ce système est assurée par une 
circulation d’huile Shell turbo 3 à 40 Engler avec un 
débit de 1 l/mn ; cette huile entretient le film et 
évacue les calories en passant dans un bac de 
15 litres refroidi par l’eau eirculant dans un ser- 
pentin. La température des paliers est ainsi stabi- 
lisée aux environs de 600, 


2830 tours/mn 
V=148 m/s 


1550 tours /mn 
= 81m/s 


1000 


1 10 100 


mm de Hg 


F1G. 11. — Variation de la puissance consommée en fonc- 
tion de la pression dans les enceintes où se trouvent les 
rotors. La distance entre un rotor et les parois de son 
enceinte sont comprises entre 50 et 100 mm. 


À l'extérieur de cet ensemble, se trouvent les 
roulements à billes en contact avec l’air et lubrifiés 
par un dispositif à vapeur d'huile SKF fonction- 
nant avec une pression d’air comprimé réglée 
de 300 à 500 grammes. 


V. Circuits électriques. — L’appareillage élec- 
trique est composé principalement des dispositifs 
nécessaires à l'alimentation automatique du moteur 
d'entraînement. 

Le groupe Ward-Leonard exécute les ordres émis 
d’un pupitre de commande, par l'intermédiaire d’un 
jeu de relais. 

Le groupe principal comprend une excitatrice, 
alimentant, sous une tension continue de 120 volts, 
lies inducteurs du moteur, de la génératrice et du 
moteur de déphasage. Le démarrage du moteur 
s'effectue sous tension variable obtenue par court- 
circuit mécanique des résistances montées en série 
dans le circuit inducteur de la génératrice. 

Sur la planche 12, figurent schématiquement 
toutes ces opérations, ainsi que les circuits annexes, 
lesquels sont : circuit de pompe à huile asservissant 
la machine ; circuit de pompe à palettes ; circuit de 
commande du moteur de déphasage avec stabilité 
automatique. 

D'autre part, des dispositifs de sécurité assurent 
l’arrêt de la machine en cas d’échauflement d’un 
des paliers ou d’accroissement d’intensité dû à une 
augmentation de puissance consommée corres- 
pondant à un défaut mécanique. 
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Pompe à huile 
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©) Excitation machine shunt 
—— Couplage mécanique 
_35- Contact de commande 
— fortes intensités 


Bac 


Sélecteur 


Géné. (Aux.) 


Groupe 
æxiliaire 


Circuit de com- 
mande automatique 


8u 3u 
= biphase 220V 220V 
Fic. 12. — Schéma des circuits électriques des machines assurant le fonctionnement du sélecteur. 


VI. Vibrations de la machine. — Les premiers 
essais avec la machine ont fait apparaître l’exis- 
tence d’un maximum accentué du niveau des vibra- 


Amplitude 
80 2 
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40 


tours /mn 
5000 6000 


4000 : 
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F1G. 13. — Amplitude des vibrations mesurées sur une des 
enceintes à vide (côté moteur). Les vérins sont en posi- 
tion extrême basse (sans huile). 


tions de la machine à une vitesse voisine de 
6 000 tours/mn, vitesse normale d'utilisation. La 
figure 13 montre la variation de l’amplitude des 
vibrations avec la vitesse de rotation. Ces vibra- 


…. 


tions sont détectées et mesurées avec un capteur 
Philips PR 9260 placé sur une des cloches à vide 
(côté moteur). Cette cloche jouant un rôle d’ampli- 
ficateur, les amplitudes ainsi mesurées sont plus 
importantes que celles des vibrations existant dans. 
les paliers. Diverses courbes ont été tracées en 
variant la position du capteur. Elles gardent lallure 
générale quelle que soit la position du capteur. 
Elles sont toutes caractérisées par trois pics prin- 
cipaux vers 1 500,3 000 et 6 000 tours/mn. Les fré- 
quences de ces vibrations correspondent essentiel- 
lement à la vitesse de rotation. ; 

Une étude systématique a été entreprise pour 
rechercher l’origine de ces vibrations. Les maxi- 
mums observés ne peuvent correspondre à des 
vitesses critiques des rotors, car ces vitesses cri- 
tiques ont des valeurs bien supérieures à celles 
correspondant aux vitesses de rotation. Ceci a été” 
confirmé par le fait qu’une amélioration de l’équi- 
librage dynamique des rotors et accouplements est 
sans action sur les vibrations. On a de même vérifiés 
que le moteur n’avait pas de vitesse critique dans 
le domaine de vitesse utile. Il semble donc bien 
que l’on soit en face des résonances du bâti de la 
machine. Ceci a été vérifié par la mesure de la 
réponse du bâti à des excitations artificielles de« 
fréquence variable. 


+ 
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L'élimination de ces résonances du bâti semble 
difficile car il est difficile de changer la rigidité du 
système. Deux choses ont alors été faites : a) calage 
soigné de la machine pour éviter toute oscillation 
d'ensemble du bâti ; b) amortissement des vibra- 
tions, en faisant reposer la machine sur des vérins 
hydrauliques. Une mince couche d’huile suffit à 
Pamortissement. Les résultats obtenus sont résu- 
més par la courbe de la figure 14. Lé niveau moyen 


2000 4000 5000 6000 


3000 


F1G. 14. — Amplitude des vibrations, les vérins étant sous 
pression. La différence entre les figures 12 et 13 fait 
apparaître l’effet de l’amortissement. 


de vibration ainsi obtenu est satisfaisant surtout 
si l’on tient compte du fait que le niveau effectif est 
beaucoup plus faible sur les paliers. 
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F1G. 145. — Vue d’ensemble de la machine. 
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Prochainement, des essais seront effectués pour 
changer la rigidité du bâti. 


Conclusions. — Le montage et les essais prin- 
cipaux ont été terminés en mars 1957. 

Certains éléments ont subi quelques mises au 
point, principalement : les garnitures tournantes 
d'étanchéité fournies par la Société PACIFIC ; les 
roulements à billes, oscillants à double rangée, du 
type SKF à cages en acier embouti, ont été rem- 
placés par des roulements SKF de fabrication 
suédoise, à cages massives en bronze. 

Les mesures des vibrations ont mis en évidence 
l'intérêt du montage sur vérins hydrauliques, mais 
l’expérience a montré que leur fidélité de réglage 
avec le temps n’est pas certaine ; actuellement, un 
système permettant d'obtenir les mêmes caracté- 
ristiques dynamiques que celles des vérins, mais 
conservant de manière rigoureuse la position dans 
l’espace, de la machine, est à l’étude. 

L'appareil installé à la sortie d’un canal de la 
pile EL 2 depuis mars 1957 a donné toute satis- 
faction, montrant qu’un tel système à deux rotors 
liés mécaniquement est parfaitement viable. Par 
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rapport à un spectromètre à cristal, ce système pré- 
sente une moins bonne résolution, mais une inten- 
sité plus forte. Le choix entre les deux techniques 
dépend du type de l’expérience envisagée. 

Nous tenons à remercier le C. E. M. H. pour toute 
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l’aide apportée dans les mesures (mesures des con- 
traintes, équilibrage et vibrations) et par les fruc- 
tueuses discussions que nous avons eues en parti- 
culier avec MM. Chaudron, Darricaud et Platzer. 
Manuscrit reçu le 23 janvier 1958, 
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INSTALLATION D'UN DISPOSITIF D'IRRADIATION 
POUR 1 KILOCURIE DE RADIOCOBALT 


Par Marc LEFORT et Xavier TARRAGO, 
avec la collaboration technique de M. VIAUX 


Laboratoire de Physique Nucléaire de la Faculté des Sciences de Paris à Orsay 


Résumé. — On décrit les dispositifs d'irradiation et de protection installés pour recevoir 1 kilo- 
curie de radiocobalt, destinés à des recherches sur les effets chimiques des rayonnements y et sur 
les réactions photonucléaires. Les échantillons de cobalt radioactifs sont placés dans six petits 
cylindres fixes disposés selon un anneau enrobé dans un tambour de plomb. Les diverses possibilités 
d'irradiation sont décrites et on indique les doses intégrales par heure en diverses positions. 


Abstract. — Design of a kilocurie $°Co gamma ray source is described. Irradiation facilities 
are given both for research on radiation chemistry and photonuclear reactions. Cobalt samples 
were placed in six small cylinders in a lead drum. Dose-rates are indicated for various exposure 
chambers. The volume of the central exposure chamber is of 4 000 cc and the integral dose rate 
is of 2.108 rads. kilog/hour. Tubes or electrical leads connected to external apparatus can easily 
be introduced with the sample. 


L'utilisation du radiocobalt pour l'étude des tions, gaz ou solides, et de permettre l’étude de cer- 
effets des rayonnements gamma est de plus en plus  taines réactions nucléaires. 
répandue [1]. Un dispositif pouvant recevoir 1 kilo- : 
, curie vient d’être installé au Laboratoire d'Orsay I. Sources de rayonnement et protection. — 
dans le double but de poursuivre des recherchessur L'activité totale est répartie dans six petits 
les effets chimiques des rayons y sur diverses solu- cylindres de laiton contenant des disques de radio- 


_ 


3 


Fic. 4. — Vue de la protection, le dispositif de levage des cylindres étant dégagé sur le côté. 
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cobalt de 120 à 135 curies fournis par le centre de 


l'énergie atomique britannique. Ces cylindres de 


19 mm de diamètre et de 25 mm de hauteur sont 
surmontés d’un petit champignôn de laiton facile à 
saisir au moyen de pinces à distance au moment du 
transfert depuis les récipients de transport à l’ins- 
tallation définitive (*). Les sources ont été intro- 
duites à l’intérieur d’une protection de plomb 
placée sur un socle de béton lourd dans un bâtiment 
léger spécialement étudié pour l'installation du 
radiocobalt. Le plomb antimonieux a été coulé 
dans une forme cylindrique en tôle de 5 mm 
d'épaisseur, de diamètre 90 cm et de 55 cm de 
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hauteur (fig. 1). Des canaux verticaux de diverses 
dimensions avaient été ménagés soit pour recevoir 
les sources en six positions réparties annulairement 
autour d’un cylindre de 8 em de rayon, soit pour 
faire coulisser des cylindres contenant des chambres 
d'irradiation. Les parois des tubes de laiton cons- 
tituant ces canaux ont été étamées de façon à 
obtenir une adhésion parfaite du plomb qui a été 
coulé avec des précautions spéciales pour éviter les 
retraits et la formation de bulles. Cet ensemble 
de 3 tonnes repose sur un socle de béton lourd. Le 
béton lourd a été coulé dans un montage en tôle 
comprenant des gaines d’accès en laiton (fig. 2). 


Fire. 2, — Montage noyé dans le béton lourd permettant de l’extérieur le guidage des sources. 


Ce montage ménageait au centre un puits central 
cylindrique de 10 cm de diamètre et de 25 em de 
profondeur et six tuyaux venant déboucher avec la 
précision la plus grande possible, sous l’empla- 
cement des six sources radioactives permettant de 
les atteindre de l’extérieur sans modification de la 
protection. Un demi boudin torique en plomb 
s'enfonce dans le béton et protège ainsi contre le 
rayonnement qui pourrait être diffusé horizon- 
talement entre le béton et le plomb. 


IT. Dispositif d'irradiation. — Les six sources 
sont placées dans six canaux verticaux fermés vers 
le bas par 10 cm de plomb et vers le haut par des 
carottes de 40 cm de long, amovibles et soigneu- 
sement ajustées pour limiter au maximum les 
fuites de rayonnement. Ces carottes peuvent être 
remontées à distance et les sources habituellement 


(:) Nous remercions M. Cottin qui nous a aidé pour cette 
installation délicate. 
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fixes peuvent être occasionnellement poussées vers 
le haut grâce à des flexibles placés dans les gaines 
qui traversent le béton lourd. 

On peut irradier en introduisant les échantillons. 
dans des cylindres de différents types qui coulissent 
dans les canaux ménagés dans le plomb (fig. 3). 
19 Au centre de l’anneau formé par les six sources 
coulisse un cylindre central dans lequel on peut 
irradier un volume d’environ 1 litre. 20 Entre les 
sources sur le même cercle et sur un cercle de plus 
grand diamètre peuvent être introduits des échan- 
tillons de quelques em°. 3° Enfin on peut dans 
certaines circonstances découvrir un faisceau verti- 
cal de rayons y irradiant en dehors de la protection: 

1° Le cylindre central de 83 cm de haut et de” 
9 cm de diamètre est un tube comportant unes 
niche de 20 cm de haut placée entre deux parties em 
plomb de 30 cm d’épaisseur. En position d’irras 
diation la chambre centrale est placée au niveau 
des sources et la partie inférieure en plomb pénètre 


% 


EL 
: 


No 


… dans le trou pratiqué dans le socle de béton. Pour 


placer les objets à irradier on soulève le cylindre 


« jusqu’à l’affleurement de la chambre d'irradiation 
… au-dessus du château de plomb, et la partie de 
… plomb inférieure arrête le rayonnement. On peut 


placer des volumes d’environ 1 litre dans cette 
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chambre d'irradiation. Le passage de canalisations 
ou de fils électriques à travers le bouchon de plomb 
supérieur a été prévu. On peut ainsi placer un 
thermostat, ou un vase Dewar et irradier à tempé- 
rature constante où à basse température avec une 
dose relativement uniforme. 


F1G. 3. — Coupes et vue de dessus de l'installation. 


20 Des cylindres analogues au cylindre central 
mais de diamètre très inférieur (16 mm intérieur) 


. permettent de descendre au voisinage des sources, 


des objets de volume plus faible où ils sont soumis 


» à des doses élevées. Six de ces cylindres ont des 
» positions intercalaires entre les tubes contenant les 
sources. Douze sont placés plus à l'extérieur où les 


doses reçues sont plus faibles. : 
30 Nous avons également pensé pouvoir utiliser 
le faisceau de rayons y provenant d’une source de 


. radiocobalt placée au fond de son puits dans lequel 


on a enlevé le bouchon de protection. Ceci permet 


de nombreuses expériences aussi bien physico- 
chimiques que physiques. Des secteurs de plomb 
sont placés autour du trou contenant la source 
avant relèvement du bouchon. On remonte ensuite 
à distance ce bouchon au moyen d’un câble passant 
sur une poulie et on fait glisser horizontalement 
l’objet à irradier jusqu’au-dessus du trou par une 
fenêtre ménagée dans le château. On peut ainsi 
placer dans le faisceau une cellule potentiométrique 
ou polarographique, un scintillateur ou tout autre 
appareil solidaire d’un ensemble que l’on ne veut 
ni détacher ni irradier. La position même de la 
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source de cobalt peut être ajustée à diverses dis- 
tances de la sortie du tube grâce aux gaines qui 
coulissent dans les six canaux coulés dans le béton. 


III. Dosimétrie et efficacité de la protection. — 
La dosimétrie chimique a été faite avec des solu- 
tions de sulfate ferreux 107% M dans l’acide sul- 
furique 0,8 N [2]. Les ions ferriques formés ont été 
dosés par spectrophotométrie à 304 mu en prenant 
des coefficients d'extinction compris entre 2 094 
et 2180 selon les variations de température 
de 20 0C à 25 0C. On a pris 15,6 équivalents oxydés 
par 100 eV [3]. 


a) Cylindre central : 


La dose est répartie de façon uniforme à 15 % 
près dans un rayon d'environ 4 em et sur 10 cm de 
hauteur. Le débit de dose intégrale pour 1 
litre est de l’ordre de 2.10$ rads.kg/heure 
(10,2.10?1eV/kg/heure). 


b) Positions entre les sources : 


Pour les six positions intercalaires entre les six 
sourées on a trouvé un débit de dose variant 
entre 1,40.105 et 1,6.105 rads/heure, soit 
8,7 + 0, 1 .1018 à 9,6 + 0, 1018 eV/g/heure. 


c) Cercle extérieur : 


Le débit de dose a été trouvé égal à 3,3.104 
rads/heure, ou à 1,4.10#rads/heure selon la position 
des canaux dans l’alignement des sources ou dans 
les intervalles. 
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Mesures de protection. — Des mesures concernant 
les doses reçues à l’extérieur du dispositif ont été 
faites par le service de protection du laboratoire. 
La dose varie entre 0,08 mr/heure et 12,5 mr/heure M 
suivant les différents points de la surface du 
château de plomb. Au moment de la mise en place « 
des objets à irradier dans un canal central elle 
est de 125 mr/heure juste au-dessus du canal. 
À 50 cm, la dose est de 0,08 à 0,16 mr/heure. Un 
séjour prolongé y est possible, étant donné que la 
dose permissible actuellement admise est deu 
7,5 mr/heure. Nous avons vérifié que des appareils. 
ordinaires de comptage 8 et y (tube Geiger- Muller) dl 
pouvaient être installés dans la pièce même du 
radiocobalt et que le mouvement propre est d’envi-} à 
ron 30 coups minute. Li 

Monsieur Joliot a autorisé l'installation de cette 
source dans un local spécial des nouveaux labora- 
toires de Physique Nucléaire de la Faculté des 
Sciences à Orsay. Nous le remercions vivement 
Nous sommes particulièrement redevables à 
M. Desgranges de tous les conseils techniques qu'il 
nous a prodigués pour l’exécution d’une autr 
installation de radiocobalt au laboratoire Curie 
conseils qui ont été très utiles pour élaborer les 
plans du dispositif actuel. Nous adressons égale 
ment nos remerciements au bureau d’étude et 
l'atelier du Laboratoire à Orsay où l’appareillag 
a été entièrement réalisé, et au Service de protec 
tion dirigé par M. Bonet-Maury qui a contrôlé @ 
surveillé avec soin le personnel chargé de PS 
lation du radiocobalt. 
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